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THE MAXIMUM PRINCIPLE OF L.S. PONTRYAGIN 
IN OPTIMAL-SYSTEM THEORY. PART III* 


L.I. Rozonoér 


(Moscow) 


The paper deals with questions related to the proof and employment of L.S. Pon- 
tryagin’s maximum principle in optimal-system theory. Certain new results are also 
contained in this work. The first part of the work deals with the optimization problem 
for the case of free trajectory right ends. In the second part, the maximum principle 
is formulated for a more general type of boundaryconditions. In the third part, the 
connection between the method of dynamic programing and the maximum principle is 
established, a method of solving the optimization problem in discrete linear systems is 
given, and a number of considerations are presented concerning the use of the maxi- 
mum principle for solving a definite class of problems which are related to the theory 
of dynamic accuracy of control systems, 


4, Relationship of L.S. Pontryagin’s Maximum Principle to R. Bellman's Dynamic 
Programing Method 








We show, in this section, that with certain requirements of a general character there exists a close relation- 
ship between the equations deriving from the maximum principle and the corresponding equations from the theory 
of dynamic prograr ing [1-3]. 


a) The Optimality Principle 





Initially, we give a brief presentation of the idea of the dynamic programing method, using as an example 
a problem analogous to that considered by us in section 2 of part L_ Let it be required to control the system de- 
scribed by the equations 


g, = f(z, 8,0, imi,...cme e@ed, (4.1) 


such that a given function, F (x (T)] of the state of the system at a fixed time T will be minimized, given that 
the system initially starts at point x’ at time ty, We denote the optimal control by u*(t) = u*(x°, to,t) and the cor- 
responding optimal trajectory by x*(t) = x*(x’,ts,t). We now consider the position on the optimal trajectory of 
the representative point x' = x*(x°,tg,t’) at some moment of time t = t’, ty < t' < T, and we pose for ourselves 

the same problem as was formulated above, but with the initial time and initial state now taken to be, respec- 
tively, t' and x’. It is now necessary to find the optimal control u**(t) = u**(x’,t’,t) which minimizes the func- 
tional F [x (T)] if the system is in state x’ at initial time t’. 


It is easily seen that the control u*(t), which is optimal on the time interval (ts, T) in the problem with in- 
itial state x®,t), must also be optimal on the interval (t',T) in the problem with initial state x’,t’, i.e., the values 
of the functional F [x(T)] for the controls u*(t) and u**(t) must be equal. 


—__ 


*Parts I and II of this work were published, respectively, in Automation and Remote Control 10 and 11 (1959). 
[See English translation]. : 
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Indeed, if this assertion were not correct, and the functional's value for control u**(t) were less than its 
value for control u*(t) with the system initially at x',t', then control u*(t) could always be “improved” in the 
problem with initial state x°,t»,by replacing it by the control 


u(t) texst<t’, 
v (7) _ * ’ . 

uv(t) sts, 
which would first take the system to point x‘ and thereafter to the corresponding final state. But since control 
u*(t) is, by definition, optimal in the problem with initial state x°,t, and cannot be improved, we arrive at a 
contradiction, and thus prove our assertion. This very simple fact is the essence of the so-called "optimality 
principle” which, following R. Bellman ({1], p. 83), we formulate in the following way. 


An optimal control has the property that, whatever the initial state and initial control are, the remaining 
control must constitute an optimal control with regard to the state resulting from the first stage of the control. 


Despite its simplicity, the optimality principle in many cases allows one to obtain an equation which ul- 
timately defines the optimal control. 


b) Basic Equation of the Dynamic Programing Method 





We now derive the proper equation for the problem we are considering. We denote by Sp (x®,t,) the value 
of the functional provided by optimal control u*(t) (t) < t = T) with the system initially at x°,te. In our case, 
Sy (x°,ty) = F (x*(x°,ty,T)]. By varying the initial data, x°,t», of the problem, we obtain a function S-7(x°,t,) on 
the points of the (n + 1)-dimensional space L(x}, . . . , x,t). When it leads to no ambiguities, we shall omit 
the superscript 0 and write the function S in the form S- (x,t). 


As before, we consider the state x',t' on the optimal trajectory x*(t) where x*(t,) = x°. When motion is 
along an optimal trajectory from state x°,t, the functional assumes the value S7(x°,t,). By virtue of the op- 
timality principle, controls u*(t) and u**(t) on the segment (t',T) correspond to one and the same optimal value 
of the functional, equal to S;(x’,t'). But since the functional is defined only by the terminal position of the 
representative point (at time t = T), we have that 


Sy (x, to) = Sr (z’, t’). (4.2) 


We consider the optima! control u*(t) on the segment [t,,t']. It is obvious that, for optimality of the con- 
trol u*(t), it is necessary that, to the point x’ at time t’, there correspond the least value of the function Sy (x’,t’) 
in comparison to all other values which could be obtained for all other states reachable from x’,t,,at time t’ by 
means of admissible controls u(t)@U. Taking (4.2) into account, we can write this requirement in the following 
way: 


Sz (x°, to) = min Sy (z’, t’), (4.3) 
u()eU 
i<t<t’ 


bearing in mind for this that x’ is a functional on u(t) and, moreover, depends on x’, ty. 


We now set t’ = t, + r, where r will be considered to be quite small. Then, bearing in mind that x*(t) = 
= x" and xf(te) = fi {x",ty,u*(t,)], we write 


z, = x; (to + t) = 29 + +f, [2°, lo, u* (to) +8, 


where ¢ is of a higher order of smallness than r. By assuming the existence and continuity of the partial deriva- 
tives of the function S-; (x,t), we can write 


OS 7 (2%, to) Sy (2°, to) 


Sz (2’, t') = Sr (x, to) + > ear awe (x, — 2?) + a oS (t’ — ty) + B, 
1 a | 
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where 8 is again small in comparison with r. By using the previous formula for Xj, and by substituting the ex - 
pression for S-(x',t") in (4.3), we get 


P 18 . 
Sr (x°, tp) = min [Sr (x°, to) +1 >; aS (=? te) 7 [z®, to, w° (to)} + 
ele  & 


OS > (x, to) 
+t —— + "| ‘ 
The only terms depending on control u(t) are those containing the f; and, in addition, the quantity y which is of 
a higher order of smallness than r, Therefore, the terms Sp(x*,t,) and _ Sy (=, to) can be "taken out from in 


at 
front of” the “min” sign. By simplifying and dividing by r, we get ‘ 
n 
aS (z®, to) 1OS_ (x®, to) 
et ——_—_— f,|z°, to, w (t t}. 
5, min, [a Al, toe OL + 4 
betel, +t 


This equation, just as (4.3), defines the choice of u*(t) on the segment [ts,t']. By now letting r approach 
zero, which also reduces the [tp, t, + r] segment to zero, we nee 7 an equation which determines 
choice of u*(t) at the single point t = ty. Therefore, by letting u*(t) = u*, and bearing in mind that lim + =0, 

+0 


we obtain, finally, 


OS (2%, to) 5 5r (x®, to) 


2°, to, u®). . 
Ato weU | az? fi( 0 u°) (4.4) 


This relationship is valid for any x’,tp. Therefore, the subscript (and superscript) 0 will henceforth be fre- 
quently omitted. Furthermore, by using the evident relationship max(—) = —min» (where » is an arbitrary 
function), we will use Eq. (4.4) in the form 


aS (2, t) . ae t) 
a = max \(—— ———) f,(@, t, u). (4.5) 


Equation (4.5) is essentially a specific partial differential equation which can be integrated if the corres~- 
ponding boundary conditions are given. The quantity u can be eliminated from (4.5) if one makes use of the re- 


sy (z, t) 


3, ) /,(z, ¢t, w) must be maximal for u for any values of 


quirement according to which the sum >) (— 
1 

x,t and 0S;/oxj. The elimination of u from (4.5) is carried out in exactly the same way as the elimination of 

u(t) from the function H in L.S, Pontryagin's maximum principle (cf. section 2, part I and section 3, part Il). 

For this, u is expressed in terms of x,t and 0S7/ dx;: 


aS p 
u=1(z, t, —s ). (4.6) 
For example, in the case of the equation 
aS asr . 
SE = max|— 5 (uz uz, +2)—SEu | 
as as: aS, 
with the assumption that 8S-/ xg > 0, the maximum of the expression— — xu — = wf — Sa, as a function 
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asp 


of u is attained when u = — _ e , Sothat after u is eliminated we obtain the equation 
Ox, 
aS 7,2 
rd) _ Sr, 
ot 1 8S > ox, ** 
Ot 


The boundary con~*tions for Eq. (4.5) are very simply found: for t = T, state x,t is simultaneously initial 
and final, and the functional's value does not depend on the control, and equals 


Sp (x, T) =F (x, ...,%n)- , (4.7) 


Now, the solution of Eq. (4.5) can be determined, Indeed, from (4.6) the derivatives dS-(x,T)/ dx; are de- 
termined for each value of x and, consequently, the right member of (4.5) is also determined. In the same way, 
we learn the value of the partial derivative 0S--(x,T)/ dt and, for sufficiently small r, we can find the value of 
Sy (x,T—T) in the form 


as, (z, T 
ceo Dod al 


Sr (xz, T —1t) = Sr (z, T) — 5 


Thus, moving step-by-step, we can determine the function S7 (x,t) for any yalue of t. If the function 
S77 (x,t) is known then, by means of Eq. (4.6), u is also determined as a function of x,t: 


u = 9(z, t). (4.8) 


Thus, we have uncovered the partial differential equation and the corresponding boundary conditions* 


= 


at ueU “ 
Sr (x, T) = F (zx), 


Sy (2, 1/  OSp(z, 
os me max >) (— aren) f, (x, w, t), (4.9) 


which allow us to determine the function S; (x,t) and the min-optima! control as functions of x,t. An analogous 
equation can be obtained for the case when one is seeking a maximum of the functional F[x(T)): 


Sy (x, t) arf Sp (, t) 
——eo = min 3}(——S— a, u, t), (4.10) 


Sy (x, T) =F (2). 


The arguments leading to Eqs. (4.9) and (4.10) can in no way be considered as proofs, They should be con- 
sidered as strictly heuristic, allowing us to surmise how the problem must be solved. However, as will be shown 
below, with certain assumptions the validity of the final results [i.e., an equation of the type of (4.9) or (4.10)] 
can be proved, 


c) Relationship of the Maximum Principle to the Dynamic Programing Method 
The structure of the right members of Eqs. (4.9) and (4.10) recalls the structure of the function 


*In the analogous equations [1], the partial derivatives were takenwith respect to T, rather than with respect 
to the initial time t, which leads to a change in sign in the formulas. 














n 


H = >; p,j, (x, u, 1) in the maximum principle. We convince ourselves that this similarity is not a casual one 
1 
by theorem 6, given below, which is formulated for the case of the problem which we considered in section 2, 


n 
part I [i.e., the problem of optimizing the functional S = Y cea (7)}. 
1 


We introduce the following notation: L is an n-dimensional space with coordinate axes (x, .. . , Xp,t)s 
(x,t) is a point of space L; u*(t) = u(x", to, t) is a min-optimal (max-optimal) control when the system's representa - 
tive point is at position x° at initial time t = tg; x*(t) = x*(x*,to,t) is the corresponding optimal trajectory (so that 


n 


x*(t,) = x"); S-y (x,t) is the value of the functional S = > c;z,(7T) on the min-optimal (max-optimal) control 
1 


if, at time t, the system is found at point x. 


n 
Theorem 6. In the problem of the min-optimization (max-optimization) of the functional 5 = > ea, (T) 
1 


let the function S7-(x,t) be continuous and continuously differentiable in the region I of space L. 
Then: 


1) For all t for which (x*(t),t)€ If, the min-optimal (max~-optimal) control u*(t) satisfies the maximum 
(minimum) condition with respect to p(t) = [p;(t), .... Pn (t)], where 


aS +r [x* (t), ¢] 
Oz, ’ 





P; (¢) =— 


to) eee * 


and where 


eS [2* (0), 


=; x" (1), p(), u°(), n= Da MAle (t), u* (0), t; 





2) The function S-(x,t) satisfies, in Tr, the partial differential equation 


Sy (x,t) 3 +e i 
“i = mex D(— TRE) my 0 


ae al 2 = mia )}(— o aE) A (a, u, >, 


ueU 4 
n 
where S; (xz, 7) = >) cai, if (x,T)¢€ Ir. 
1 


Theorem 6 has a graphic geometric interpretation. In section 2, part I, we recalled that the vector p(t), 
with respect to which optimal control u(t) satisfies the maximum (minimum) condition, defines the direction of 
"maximum acceleration” of the system in phase space X. As follows from theorem 6, this direction, at each 
fixed moment of time t, is determined in its turn by the gradient of the function S--(x,t) in phase space X. It 
thus turns out that, at each fixed moment of time, whatever the position of the representative point, the optimal 
control tends to “accelerate” the latter in the direction defined by the gradient of the function S-y (x,t). 


An analogous partial differential equation can also be obtained in the more general problem considered in 
section 3, part IL However, the boundary conditions here are significantly more complicated, and we shall not 




























































consider them here. In any event, the relationship of the function p(t) to the gradient of the function S is retained 
even in this more complicated case.* 


d) Analogy with the Equations of Analytical Mechanics 





The analogy of the equations deriving from L.S. Pontryagin’s maximum principle with the canonic equa - 
tions of analytical mechanics was already cited in section 2, part I and section 3, part IL The establishment of 
the relationship between the maximum principle and the dynamic programing method shows that this analogy is 
very profound. The equations deriving from the dynamic programing method are completely analogous with the 
Hamilton-Jacobi partial differential equations, whereby the function S (x,t) plays the role of the action function. 
It thus turns out that the connection between the equations of the maximum principle and the equations of the 
dynamic programing method is analogous to the relationship between the Hamilton canonical equations and the 
Hamilton-Jacobi equations. The sole difference is that, in optimal system theory, the “control” is explicitly in- 
troduced, whereas it is eliminated in analytical mechanics, and is expressed in terms of the coordinates, their 
derivatives and time. As physical analogies of the control we can use, in particular, the velocities of the system's 
mass points. 


e) The Role of the Requirement of Continuity 





In the present work we have always considered problems in which the set U, defining the limitations on 
the class of admissible controls, does not depend on the running value of the coordinates. Otherwise 
the maximum principle, in the form in which we have formulated it, turns out to be inapplicable, whereas the 
equations of the dynamic programing method remain valid in general. This circumstance is explained as fol- 
lows. The vector p(t) which enters into the formulation of the maximum principle has always been considered 
by us as a continuous function of time. But if the derivatives 0S (x,t)/ xj are not continuous in the entire space 
L and, moreover, if the optimal trajectory (x(t),t); to = t = T intersects a surface of discontinuity of the functions 
OS (x,t)/ Oxy in L, then the functions p; (t) = — dS-y (x (t),t)/ Ax; also experience a discontinuity. Therefore, it is 
necessary to generalize the maximum principle so as to allow us to consider discontinuous (in particular, piece- 
wise-continuous) pulse functions. Such a generalization was made by R.V. Gamkrelidze [4] for the case of op- 
timization with account taken of limitations placed on the system coordinates. The conditions for jumps in the 
functions pj (t) were formulated in [4]. 


Thus, the requirement connected with the continuity of the function Sy (x,t) and its derivatives has a non- 
formal character. If we follow the mechanical analogy, the requirement that p(t) be continuous corresponds to 
the requirement of no impulse skips, i.e., the absence of “impact” interaction. 


Whether or not the function Sy (x,t) itself is continuous is particularly important. The fact of the matter 
is that the corresponding partial differential equation is valid, generally speaking, only in the regions of contim- 
ity of the functions 8S-7/ dx; in the space L, The condition that Sy (x,t) be continuous allows one to “join” the 
solutions at the surfaces of discontinuity of the derivatives dS-;/ xj. The author knows of only one work [5] in- 
vestigating such questions in which continuity is proven for the optimal values of the functional (time control) 
in the linear problem of temporal optimization with limitations on the control of the type |u,,| <1. 


We note further that the existence of singular controls (cf. section 3, part Il) is also related to the presence 
of discontinuities of the derivatives 8S7;/dx;. A singular trajectory occurs precisely along a surface of discon- 
tinuity. 


f) The Synthesis Problem 


As has been mentioned, the knowledge of the function S-7(x,t) allows one to find the control u in the form 
of function (4.8), which rapidly gives a complete solution to the synthesis problem. However, solution of partial 
differential Eqs. (4.9) and (4.10) in closed form is possible only in the most elementary cases. In practice, there- 
fore, one should attempt to solve the synthesis problem by means of the dynamic programing method only if one 
has access to a computer. The development of convenient machine algorithms for solving these equations is, 
therefore, of importance. 








*In the general case, the surface S(x,t) = const is the geometric locus of the points which give the same value 
to the functional when the optima! trajectory begins with them. In particular, in the case of maximizing speed 
of response, the corresponding surfaces T(x) = const are the “isochrone” surfaces first considered by A. Ya. Lerner. 





con 











5. Concerning Optimal Processes in Discrete Systems 





In this section we shall consider several optimization problems in discrete automatic control systems, de- 
scribed by finite difference equations. For this, we assume time to be “discrete,” i.e., time takes the values 
0,T,...,™mrT,... where r is some constant having the dimensionality of time. The control processes in such 
systems can be described by the following system of finite difference equations: 


eet x thi (a, vee y TMU... Um, mt) (ie d,... 5m) (5.1) 
The symbols xf", uy uj, are the abbreviated notations for the variables x; and uz considered at time mr, i.e., 


= x;(mr), wi ujy,(mr). In the sequel we shall consider m as the time variable, assuming the integral 
Se 0, 1, 2, 


In system (5.1) we have, as before, denoted by x; (i= 1,...,m)andu,(k=1,..., 4) the generalized 
system coordinates and, respectively, the controlling actions. By introducing the vectors x = (xf?,... , |e 
and u™ - (uy, ras ur’), we can write system (5.1) briefly as 

emtl — 7™ = tf (x™, u™, m),. 


We shall consider as given the value of vector x™ at time m = 0, and shall denote it by x*. The space 
X (xq, . . . , Xn) in which the generalized system coordinates can vary is called the phase space. A “trajectory” 
of a discrete system in the phase space is a sequence of isolated points. 


As before, we shall call the vector u™, considered as a vector function of time m, the control. The limita - 
tions on the controlling actions are expressed by requirements in accordance with which, at each moment of time 
m, the control u'™ must lie in a closed set U of space R(uy,... , Ur): 


u™ ¢- U (ma = 0, 1, a on © (6-8) 


With respect to system (5.1), we can pose various optimization problems, completely analogous to those 
formulated for continuous systems in section 1, part I. The sole difference is that, in the proper places, the in- 
tegrals over the continuous time t are replaced by sums over the discrete time m, and derivatives with respect 
to t are replaced by the differences of variable values at neighboring instants of time m. A wide class of op- 
timization problems (in complete analogy with the problems for continuous systems) can be reduced : the fol- 
lowing basic problem. From the set of admissible controls u'™ €'U which take the system from point x* to a fixed 
closed set G of phase space X, it is required to so choose the control u'™ (m = 0, 1,..., M1) that the sum 


n 


S= Ds cx ata given moment of time m = M, assumes a minimum (maximum) value. 
1 


The control time here is assumed fixed. Problems in which the time is not given beforehand will not be 
considered further here. 


We remark that the probiem formulated above _ eae the — of minim Pe wae the 
function S with respect to rM independent variables u®, ... , ug; uj,..., uj; u uM-1 1 Indeed, system 
(6. 1) is a systern of recurrence relationships which allow one to daterusias mccevely the ries of the vectors 
xt x2. x™ in the form of functions of the coordinates of the vectors u°, ut, .. . “4 if just the initial 
value x® is given. By thus obtaining the quantities xM (i=1,..., Mm) and the function . itself as functions of the 
variables uf, ..., uM- 1 one can determine the extrema of s with the conditions that u™U (m = 0, » M~1) 
and x™(G, However, such a direct method of solving the problem in any complicated case turns out to be in- 
admissible in practice, due to the mass of computations required, even if a computer is available. Moreover, 
the necessity arises of providing a uniform method of solution valid for all problems of a given class and allowing 
the use of a uniform programing scheme. It is, therefore, necessary to seek methods which provide simpler (from 
the computational point of view) solutions of the problem posed, One of these methods is R. Bellman's dynamic 
programing method. Below, we present another method, related to an extension to discrete systems of L.S. Pont- 
ryagin’s maximum principle. 


In what follows, we shall consider systems which are linear in the variables x;, ... , Xn, i.e., systems of 
the form 


n 

ott — 2 = (2 amen + 9,(uy,-..,um)| (= 4,..-0m), (5.3) 
=1 

where the functions gj are assumed to be continuous, and the coefficients aj, can depend on the time m. 


The limitation of our consideration to linear systems is explained by the fact that the extension of the 
maximum principle to discrete systems is possible, generally speaking, only in the linear case.* 


We mm now to the formulation of the basic results. We consider n functions pf’, ..., py’ of time m 
which satisfy the relationships 


of, (z™, u™, m) 


m 
pr— pr? =— Dip, Bs, (i= 4,..., n), (5.4) 
1 





where, by virme of the linearity of the functions fs, 


of, (2™, u™, m) 





Oz, = a: (5.5) 
We call the vector p™ = (pj, . . . , pi) the impulse. 
By introducing the function 
H (x, p, u, m) =t D> Pole lis 2 5g Bue Bg, .., Bp; M) 
1 


and denoting the corresponding first differences by the symbol A (Azj = z!+! — ai, Ap! = p!+1— pl, 
t=0,1, 2,...), we can write relationships (5.3) and (5.4) in the form of the system 


0H (z™, p™, u™, m) 
? 
opt" 


0H (2™, p™, u™, m) 


™ 
Ax j 








Az = Apr = (i= 1,..., n). (5.6) 


We now turn our attention to the circumstance that the values of the variables x™, p™ and u™ at time m 
define the first differences of the coordinates at the same moment of time, but the first difference of the im- 
_ pulses at the previous moment of time. In other words, by knowing the values of the variables x™, p™ and u™ 
at time m we can determine from (5.6) the coordinates at the following moment, m + 1, and the impulse at the 
previous moment, m — 1, 


Let w™€U be some admissible control, x™ = a(m) the trajectory corresponding to it and p™ = 8(m) some 
one of the vector functions of time which satisfy (5.6). By substituting in the function H the values of x™ and 
p™, which are now definite functions of time, we obtain the quantity 


K (m, u,,...,U,) =H [a(m), B(m), u, m), 


which, for each moment of time m, is a function of the point u = (uy, . . . , u,) which lies in set U of space 
R(uy,..., Uy). We shall say that control u™ satisfies the maximum (minimum) condition with respect to vector 
function p@ = 6(m) if, at each fixed moment of time m (0 =m = M-~—1), the function K(m,u) attains an absolute 


maximum (minimum) on the set U for values of the variables equal to the values of the controls at the same mo- 
ment of time, i.e., for uz = uj 








*As r - 0, the solution of the discrete problem tends, in a definite sense, to the solution of the corresponding 
continuous problem (cf. [6], where this assertion has been proven for linear systems and limitations of the form 
|uzl| +61). Therefore, for sufficiently small r, the maximum principle turns out to be applicable in a sense. 
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Just as in the case of continuous systems, we first formulate the results for the problem with free trajectory 
right ends and then turn to the more general problems. 


In the problem with free trajectory right ends, the set G occupies the entire space, which we symbolize by 
G ~X. The following theorem establishes the discrete analogy to theorem 2 of section 2, part L 


Theorem 7, A necessary and sufficient condition for min-optimality (max -optimality) of control u™ (m = 
= 0, 1,..., M1) in system (5.3) for G ~ X is the holding, for u™, of the maximum (minimum) condition with 
respect to the vector functions p™ which satisfy the relationships 


pM- -~—c, (i=f,...,A). 


The proof of theorem 7 is completely analogous to the proof of theorem 2 for continuous systems, and is 
based on the following formula for the increment of value of function S with a variation of control u™: 











n M~1 
8S =>) cisai = — >) (A (a™, p™, w™-+ du™, m) — H (2™, p™, u™, m)|—1, (5.7) (5.7) 
1 m=0 
where 5u™ = (6uy', ..., uM) are the increments of control and n = 91 + ng are certain residual terms, where 
M—i n 
os Ny) [ af (y™, u™ + Bu™, m) aH (y™, u™, m) sym 
teen? a m os m Y, 
m=0 8=1 OY, OY, 
- ns I aH (y™ + oty™, u™ + §u™, ™) 
Paap, 5 J [ re? 
~ m=~0 8,q=1 ay, Avg 
ad (y™ + oP ty™, u™ + tu™, m) 





x $ ms) m. 
ayy auy" 4 Ma 

Here, 0 < 9y< 1, 0 < Sy < 1 and, for brevity, we introduced the vector y = (yy, . . - » Yen) (yi = xy Ynet = 
= py i=1,...,.m). For linear systems of the form of (5.3), the residual term n reduces to zero. The proof of 


formula (5.7) is carried out in complete analogy with the derivation of the corresponding formula for the incre- 
ment of the functional's value. 


We now turn to the more general problem wherein the “trajectory‘s” right end lies in some set GEX. It is 
further assumed that set G is convex and closed. 


We first formulate a sufficient condition for optimality. Just as in the analogous continuous problem (sec- 
tion 3f, part Il), we shall denote by bj (yy, . . . . Yn) (i= 1,... , m) the coefficients of the hyperplane bracketing 
set G at point y = (yy,...., Yn). If y is an interior point, we then set bj (yy, . . . , Yn) = 0. For each point 


n 
X (Xj, ...,» Xn) €G in the problem of minimizing function S we require that the inequality >) bi (Yis+ ++ sq) 
1 
(xj — yy) = 0 hold, whereas we require that the inequality hold with reversed sense in the problem where the 


function is to be maximized. 


Just as in the continuous case, we shall consider nondegenerate problems (section 3a, part II). In a nonde- 


n 
generate problem, the set G*€G, on which the function $(z) = > CjX, assumes its least (greatest) value in 
1 


comparison with the values for ail other points x ¢G, is not attainable from point x® within time M for the con- 
trols u™€ U. 


Now, if x®,x4,. . . .xM is an optimal “trajectory,” corresponding to control u®, u’,. . . , u™~4, the follow- 
ing theorem, establishing a sufficient condition for optimality, is valid. 


Theorem 8. If, in system (5.3), control u™ (m = 0, 1,..., M~1) satisfies the maximum (minimum) con- 
dition with respect to vector p™, whose coordinates at time m = M—1 assume the values 
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pM = — hg — abi (2M, coos aM) (i=1,...,A), 


where y is nonnegative and A is positive, then control u™ is min-optimal (max-optimal) with respect to S = 
n 


A 
2caé. 
! 


The proof of theorem 8 is based on theorem 7 and is carried out in complete analogy with the proof of 
theorem 4, 


In the formulation of necessary conditions for optimality, it is necessary that one further requirement hold, 


this requirement amounting to the following. We consider the mapping of r-dimensional space R(uy, . . . , Uy) 
on n-dimensional space X (vy, ... , Vn), given by the functions 
Vi = 9, (Uy,...,Ur) (i=1,...,%), 


which define the right members of Eq. (5.3). We denote by V the set of all points of space X which are images 
of points which lie in set U of space R. We shall require that set V be convex (which occurs, in particular, if 
vq = bgu and the variable u is of bounded modulus). 


The following theorem establishes a necessary condition for optima lity. 


Theorem 9. Let set V be convex. Then, if control u™¢ U (m = 0, 1,... , M~1) in the nondegenerate 
n 
problem is min-optimal (max-optimal) with respect to S = DSor™ , there then exist vector functions p™ with 
1 


respect to which control u™ satisfies the maximum (minimum) condition, where the coordinates of vectors p™ 
at time m = M~1 assume the values 


pM" = — hey — whi (eM, ..., aM) (i= 4,.-..m), 


where the constants A and » are nonnegative and do not vanish simultaneously. 


The proof of theorem 9 is obtained in the following way (for example, for the case of min-optimality). 
By the conditions of the theorem, the set Wy of the points of space X which are attainable in time M with the 
n 
control u™'-U is a convex set. The set G~, which is the set of those points x(-G for which dei (xj — a") <0, 
1 


is also convex, where the points common to G" and Wy (in particular, point xM) can only be boundary points of 
the two sets. The foregoing defines the hyperplane A which separates sets G’ and Wy and which is a bracket for 
set Wy. By so choosing the signs of the coefficients aj (i = 1, . . . , n) of hyperplane A that the inequality 
n 
ya (xi — ai" )>0, holds for any point x(-W7, we can convince ourselves that control u™ is min-optimal for 
i 

n 


sS= Yaai in the problem with free right ends and, thanks to theorem 7, satisfies the maximum condition 
1 


with respect to the vector function p™’ which, for m = M~—1, assumes the values po-1 =—aj(i=1,...,n). 
The relationships aj = Acj + bj (xy... . xM) are obtained exactly the same as in the continuous case (section 
3, part Il). 


Theorems 7-9 are used for solving concrete problems by methods completely analogous to the method of 
using the maximum principle in continuous problems. Specifically, the maximum (minimum) principle allows 
one to express the control at each moment of time in terms of the coordinates and impulses of the system, thus 
providing the capability of eliminating the control from Eq. (5.6). The system of 2n finite difference equations 
thus obtained, with 2n variables and with boundary conditions stated in the formulation of the theorems, defines 
the optimal trajectory and impulses and, consequently, also defines the optimal control. 
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6. Use of the Maximum Principle for Solving Certain Problems Connected with the 
Dynamic Accuracy of Automatic Control Systems 








The variational methods developed in optimal system theory have immediate application to the problems 
connected with the estimation of control processes. This is particularly the case when the problem of dynamic 
accuracy is so posed that the extemal stimuli on the system are assumed to be given only by a number of limita~- 
tions: limitations on modulus, known values of the integral of the squared disturbances, etc. Indeed, the basic 
problem in the investigation -of such cases amounts to choosing, from the class of admissible external stimuli, 
that one which is the most “dangerous” from the point of view of the accepted criterion of dynamic accuracy. 
Thus, for example, in B.V. Bulgakov’s problem of the accumulation of deviations in linear systems [8-10], one 
chooses, from the class of modulus-limited external stimuli, those which lead to the maximum deviation of the 
controlled quantity. In [11, 12], the class of external stimuli is assumed to be given by definite values of the 
integral of the moduli, or the integral of the squared disturbance, and one then computes the maximum yalue 
which the integral of the squared controlled quantity might assume. It is obvious that, in all these cases, a varia- 
tional problem is essentially solved. An explicit use of variational methods and, in particular, the methods de- 
veloped in optimal system theory, holds out broad possibilities, both in the sense of the investigation of nonlinear 


systems and in the sense that the class of admissible external stimuli might be given in significantly more diverse 
ways. 


B.V. Bulgakov's problem on the accumulation of deviations can be formulated in the following way, which 


is also meaningful for nonlinear systems. Let uj, ... , uy be the external stimuli acting on the control system 
described by the differential equations 


t= fi(t,..-, In; Usscecy Uy; t) (i =1f,, - - M). (6.1) 
All that is known about an external stimulus u = (uy, . . . , Uy) is that it satisfies some system of inequalities 
(in the particular case here, |u,,| < ay, k= 1,..., 1). It is required to find the maximum possible deviation 


X(T) which the controlled quantity x, can undergo at time t = T. 


The problem just formulated is identical with the problem considered by us in section 2, part I, whose 
theorem 1 gives the capability of solving the problem of accumulation of deviations for nonlinear systems with 
diverse limitations on the external stimuli. The results already obtained for linear systems are easily obtained as 
particular cases here.* It turns out that the “impulse” vector considered in the maximum principle is intimately 
related to the impulsive response of linear systems, 


For example, we now show that, by means of the maximum principle, the simplest problem of deviation 
accumulation is solved for a system with constant coefficients, described by an equation of the form 


ph) 


x’ + Qn— 2" © 4. seul aor =u, |u| <i. (6.2) 


The relationship between the system's impulsive response and the “impulse” vector is particularly clear 
here. For simplicity of exposition, we shall write the formula for the case when n = %, 


© +- a,x + aor = u, (6.3) 


bearing in mind that the analogous relationships subsist in systems of higher order. We now write Eq. (6.3) in the 
form of the system 


Le = 2, Ly = — Apo — A, 2, + U, 
where x = X, is the coordinate to be controlled. The problem consists of finding the maximum of the function 


S = X9(T) for the fixed time T. We must, therefore, solve the problem with free trajectory right ends, whereby 
Cy = 1 and cy= 0. In correspondence with the general methodology, we set up the function 





*In those cases when the derivatives of the external stimuli enter into the system equations, we should use the 
results obtained in [4]. 
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H = Poy — Pi (Ao%o + 423) + Py 


and the system of equations 


. Hi) . 0H 
Po = -- oP, Pi = — Ze Po + SP: (6.4) 
The boundary conditions are written in the form 
Po(T) = —eo =—1, p, (T) = — ec, = 2. (6.5) 


Since we seek to maximize the functional, we use the minimum condition. The minimum of function H 
with respect to u is obviously attained at the boundary points: for py > 0, u=—1 and, for py< 0, u= +1, ie., 


u (t) = — sign p, (t). (6.6) 


The alternation of the intervals on which the external stimulus assumes the values +1 and ~1 is thus de- 
termined by the variations in time of the function p;(t). We now replace system (6.4) by one equation in p = 
= pa(t) 


p—ap- agp = 0. (6.7) 


| In (6.7) we carry out the change of variables, W(t) = — p(T—t). Since W(t) = p(T—t) and W = p(T—1), 
(6.7) is written in the form of the following equation in the function W(t): 


W(t)}+aWt)+aW()=0 © (6.8) 
with boundary conditions [by virtue of (6.5)] 


W (0)=0, W(0) =1. (6.9) 





However, as is well known, the solution of Eq. (6.8) with initial data given by (6.9) gives the system's im- 
pulsive response. The left members of Eq. (6.8) and the sought-for Eq. (6.3) coincide. Thus, W(t) is the im- 
pulsive response of control system (6.3), and component p, of the “impulse” vector (pp,p,) equals 


P(t) = —W (T —1). (6.10) 
Since po = aop; and pp(T) = — 1, it easily follows that 


Po(t)=—1+h(T—2), (6.11) 


where h(t) = \ W (z)dz, i.e., the system's reaction to a unit step function input. The most “dangerous” stimu- 
0 


1) lus, by virtue of (6.6) and (6.10), is determined from the well-known formula 
u(t) = sign W (T — 12). (6.12) 


Also easily solved are problems in which, in addition to the limitation |u| =< 1, there are limitations of 
T v 
the type \ juldt< A or \ u* dt < A .* With such limitations, one may seek, not only the maximum de- 
0 0 
viation of the controlled quantity but also, for example, the maximum value of the integral of the square of this 





*We note that, when there are “integral” limitations, we have to deal, generally speaking, with a problem in 
which the trajectory right ends are no longer free (cf. section 1, part I). 


1528 





| 











quantity (analogously to [12] where, however, T = oo and the limitation |u| = 1 is lacking), In linear cases, the 
solutions of such problems can be expressed in terms of the system's impulsive response. 


The results of section 5 indicate how to expand the above-considered methods for investigating the dynamic 
accuracy for systems’ control of discrete linear systems. 


APPENDIX III 


” 


Proof of Theorem6. For definiteness, we consider the case where a minimum of the functions > cx, (T). 


1 
is sought. The theorem is proved in an analogous way when the case of maximizing the functional is concerned. 





As before, we shall assume that the right members of the system 


a, i; (x, u,*t) (i--1,...,n, u(t) GU) (IIL. 1) 


are continuous in the set of arguments (x,u,t) and have continuous first and second partial derivatives with respect 
to the arguments (xy,... , Xp). 


Let (x°,to) be some point in L, and let u* = u*(x°,t,t) and x*(t) = x*(x*,tg,t) be, respectively, the optimal 


n 


control and the optimal trajectory corresponding to a minimum of the functional p:: e,x,(T') for a fixed value of T. 
I 


Further, at time t = t' (t) = t" = T), let the system be in state x = x'. On the segment [t’,T] we shall use 
the control u*(t). With this, the function y cx,(7') assumes the value -(x',t") which depends exclusively on 
(x',t') since the control u*(t) is given. By siete the state of the system (x’,t'), we obtain the function 

O, (z’, t') == Oy [z’, t’, (u* (a, to, t)}), 


which has for its value at each fixed point (x',t')€L the value which the functional assumes on the control u*(t) 
(t' =t = T) for the system's initial state being (x’,t'). The symbol {u*(x*,to,t)} (which will be omitted in the 
sequel) indicates that the function 7 (x",t') is meaningful only if the control u*(t) is previously given. Let 
x(t) = x(x',t’,t) be the system's trajectory corresponding to control u*(t) with initial state (x',t’). By definition, 


n 
D, (z’, t’) = 2 C2; (x, 0’ 5 2). (II1.2) 
1 


Since the functional to be minimized is defined by the final value only (at time t = T) of the system's co- 
ordinate, it is then obvious that 


D, [x (t), t] = Dy (2’, t’) (<ts7), (IIL.3) 


i.e., the values of the functional on the control u*(t) for initial states lying on the trajectory x(t) = x(x',t’,t) are 
all equal to one another. 


It follows, from our assumptions as to the right member of Eqs, (IIL 1) and from the piecewise continuity 
of control u*(t). that the solution x = x(x',t’,t) of system (1) is continuous in the set of initial data (x',t') and has 
all continuous first-order partial derivatives and continuous second-order partial derivatives #x;/ 8x’, Ox'g, 
&x;/Ox,dt' (i,s,q=1,...,n) (cf. (7).* 





*In (7), the proof of the theorem on the continuity and differentiability of the solution with respect to the initial 
data is carried out on the assumption that the right members of the differential equations are continuous in t. 

By assuming that, at the points of discontimity, u*(t) = u*(t+ 0) and by using the so-called theorem on half-in- 
terval continuity, one can easily show that the solution of system (1) of equations is differentiable the proper num - 
ber of times with respect to the initial data even in the case considered, when control u*(t) is piecewise continu- 
ous and has a finite number of first-order discontinuities. oni 














From this follow the continuity and the existence of the corresponding partial derivatives of the function 


n 

0, (2, ') = > cz, (z’, t’, T). By taking this into account, and by then taking the total derivative of Eq. (IIL.3) 
1 

with respect to t, we get, by virtue of (IIL 1), 


n 


OD, [x (z’, t’, t), t) > OO, [x (x’, t’, t), t] 
ot od Oz, 





fy lz (z’, t’, 1). u* (4), t]. 





1 


This last relationship holds identically for any x',t' and t(t' st =T). Therefore, by setting t = t’, and 
keeping in mind that x(x',t’,t) = x' and then discarding the primes, we obtain the identity 


n 


OD, (z,t) AM (z, t) 
5 = - ae 


1 8 


f, (=, w* (t), #}. (III.4) 


We now consider the functions 


dD, [z* (t), ¢] 
7, ()=- On, (ioe i, ..., m). (IIL. 5) 





We first note that, by virtue of the well-known relationships 





dx, (2, T, T) {3 i=), 


dx? 0 is] 


(cf. for example, [7]), we have the equality 
n ” 7; Al 
Sr.) we ¥ «, oz, (2°, T, 7) ™ ¢; 3 ee oo 
aa} a29 





1 
and, consequently, 
A(T) = — &;. (IIL 6) 


We now derive the system of differential equations which are satisfied by the functions p,(t). We differen- 
tiate (IIL5) with respect to t (the corresponding derivatives exist and are continuous): 














° a OD, [a* (t), t| =f HPD, [x (1), t] +e 
ee 2 Fae tt). (111.7) 
From (III,1) 
z= f,[z* (t), u* (t), t] iT aaa (IIL.8) 
By differentiating (IIl.4) with respect to xj (i=1,..., n), we get 
PDp(z,t) ay O*Dy(z, t) ; “ 
ax,at — 2 ax,Ar, fale, w° (t), & 
- OD, (z,t) af, [z, u” (t), t] 
aT, ar; 2 (II. 9) 


1 
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Now, by substituting x = x*(t) in (IIL.9), by taking (IIL.5) into account and by substituting (III. 8) and (IL. 9) 
into (IIL.7), we obtain the system of differential equations sought: 


af, {z° (t), u* (t), #] 
Oz; o 





n 
P; (t) = — >} P, (t) 
1 


er) (IIL. 10) 


But, according to theorem 1 of part I, optimal control u*(t) satisfies the maximum condition with respect 
to vector p(t), defined by relationships (III.5) which are subject to system (III.10) of equations and conditions 


(IIL.6). Therefore, by substituting x = x*(t) in (III.4) and by noting that the right member of (IIL.4) is the expres- 
sion 


HI {x*(t), p(t), u(t), t}== >) py (t) f, [= (t), u* (0), 
1 


we obtain 








Or [rt] s) | a®, [x*(t), t] 


at ueU | \ az, 


) f,{2*(t), u, t]. (IIL 11) 
In particular, for t = ty, x*(tg) = x°, 


O, (x®, t Ms a®., (x0 t 
ere mas ¥ (- ere.) 1, (a, WF, ted. (IIL 12) 
Oto u’eU ax® 


We now establish the relationship between the functions $7 (x,t) and Sy (x,t). By definition, for any point 
(x,t) in space L the following equation holds 


®, (z, t) > Sp (x, t), (IIL, 13) 


for which at points x, t lying along the trajectory we always have: 
®, (x, t)= Sp (2, t) (x = r*(t)). (TIL. 14) 


It follows from this that, in the (n+2)-dimensional space (z,x,t), the surfaces z = ®7(x,t) and z = Sy(x,t) 
do not intersect, and have the common curve y, whose projection on the subspace L gives the optimal trajectory 
x = x*(t). In other words, the surfaces z = ®(x,t) and z = S-(x,t) are tangent in the curve y. Now, if the func- 
tion S- (x,t) is continuous and differentiable with respect to (x,t) in a neighborhood of the optimal wajectory x = 
= x*(t) (the continuity and differentiability of the function @7(x,t) were noted by us above), then the following 
relationships are valid everywhere on the optimal trajectory 


aD, (zx, t) OS 7 (x, t) O®,,, (x, t) OSy (x, t) 
at amy a Ox; Oz; 











(IIL. 15) 


By assuming, in accordance with the conditions of the theorem that (x*(t),t) I, and by thus guaranteeing 
that (III.15) holds, we obtain, in correspondence with definition (IIL. 5), 


OS, [x*(t), t] 


Ox ; 





A O=— (i 4,..., A). (111.16) 


By replacing 06-~/ dx; and 067/ dt in (IIL 11) by dS-¢/ ax; and dS-p/ dt, and by using (IIL 16), we obtain 
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AS» [x*(0), | 


wae II (a*(t), p(t), u(t), t). (IIL 17) 


Since we have proven that the optimal control satisfies the maximum condition with respect to vector p(t), 
defined by relationships (111.16), point 1 of theorem 6 is completely proven. 


To prove point 2 of the theorem, it suffices to replace the partial derivatives of function ®-7 in (III.12) by 
the analogous derivatives of the function Sy, in accordance with (III.15). If we note, with this, that the equa- 
tion obtained remains valid for arbitrary x°,t), and if we discard the superscript (and subscript) 0, we will obtain 


the identity* 
OS, (x, t) As x) OS, (x, t) 
<n oe _ (- aT | f(z, u, t). (11.18) 
The condition 
i S, (z, T) = >} 7, (IIL. 19) 


1 


is obvious in view of the definition of function S (x,t), if one realizes that, by virtue of continuity, 
ia Sy (x, t) = S, (x, 7). Formulas (III.18) and (IIL 19) prove point 2. Thus, theorem 6 is completely proved. 
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*In contradistinction to (II.18), one cannot discard the subscript 0 in (IIL 12) and consider the relationship thus 
obtained as a partial differential equation defining the function ®7 (x,t). Equation (III. 12) is not identically valid 
and, as is clear from (III.11), is valid only for points (x°,t») lying on the optimal trajectory x = x*(t) which cor- 
responds to some previously fixed optimal control u*(t). 
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A FREQUENCY METHOD FOR DETERMINING THE DYNAMIC 
CHARACTERISTICS OF OBJECTS OF AUTOMATIC CONTROL 
FROM DATA ON THEIR NORMAL USAGE 


V.V. Solodovnikoy and A.S. Uskov 


(Moscow) 


A frequency method is presented for determining the dynamic characteristics of 
objects of automatic control from data obtained in the course of their normal usage, 
specifically, correlation and cross-correlation functions. 


The paper considers linear objects with many inputs and outputs and with con- 
stant lags, and also multiloop linear systems in the presence of noise. 


It is assumed for the analysis that the random processes occurring in the systems 
are ergodic and stationary, and that the objects under investigation are stable. 


1. Present State of the Question and the Posing of the Problem 


As is well known, the ordinary methods of determining the dynamic characteristics of industrial objects by 
measuring their reactions to artificial stimuli of definite forms applied at their inputs [1, 2] are inapplicable in 
many cases, for the following reasons, 





1, It is undesirable or impossible to apply a disturbing stimulus of a special form at the object's input be- 
cause this leads to a disruption of the normal course of the processes in the object. 


2. Very frequently, random uncontrolled disturbances are superimposed on these stimuli. As a result, it 
turns out to be impossible to determine accurately the dynamic characteristics by means of standard input signals. 


In view of these circumstances, great importance has recently been attached to a statistical method which 
permits the dynamic characteristics — transfer function or impulsive response — to be determined from correlation 
functions which are obtained from data on the normal usage of the objects. The present work is devoted to the 
theory of this method and to the methodology of its use. Works [3-6] present methods for machine solution of 
the Fredholm integral equation of the first type which relates the correlation function and the impulsive response. 
This equation, which is the starting point in the obtaining, by the statistical method, of the characteristics of ob- 
jects with one input and output, is written in the following way: 


Ry (2) = Re (e— 9) k (9) a, (1) 


where 


T 
R, (t) = lim + { z(t 4-t)a(t)dt 


0 





























is the correlation (autocorrelation) function of the process x(t), x(t) is an ergodic and stationary random process 
at the object's input, 


T 


Rys () = lim =| y (t+) x(t)dt 
Toco by (3) 


is the cross-correlation function of the processes x(t) and y(t), y(t) is the ergodic, stationary random process at 
the object's output. 


The sought-for impulsive response k(t) must satisfy the condition of physical realizability 
k(t) = 0 fort < 0, (4) 


In [7], using a linear presentation of the problem, solutions are given for determining the dynamic charac- 
teristics of objects with several cross-correlated inputs and outputs. Two methods of solution are presented. The 
first method is based on the separation from each input of the component which is not correlated with the re- 
maining input stimuli, In the second method, the problem is reduced to the choice of stimuli which are corre- 
lated with only one or, in the extreme case, with a small number of stimuli at the object's input. 


In [8], this latter method is used for the analysis of single-loop and multiloop systems with noise applied 
inside the objects under investigation. Finally, in [4], there is considered an approximate method of determining 
the characteristics of objects closed by feedback paths and with internal noise present. 


In all these works (3-8) solutions are carried out principally in the time domain, i.e., the determination 
of the objects’ characteristics is reduced to the finding of the impulsive responses. This, in many cases, leads to 
horrendous computations, particularly when integral equations must be solved. It tums out that the problem can 
be significantly simplified if, by using the Fourier transform (cf., for example, [9]), one seeks the solution in the 
frequency domain, t.e., if one initially determines the objects’ frequency characteristics, from which one can 
always go to the impulsive responses. 


It is true that with this it is necessary to compute the spectral density from previously computed correla - 
tion functions, On the other hand, Eq. (1) for a stable object is very easily solved by dividing the spectral den- 
sities: 


: S x (w) 
@ (iw) = 5 (wy) (5) 
where (iw) is the object's transfer function, 
bes 
Syx(w) = \ Ryx(t)e—* dr, (6) 
rng 
Sz(w) = | Re (2)e-1* de. (7) 


Solution (5) of Eq. (1) is given in (8-10). In the first work it is given without proof, and in the two latter 
works it is obtained formally by means of the Fourier integra] without the condition of physical realizability (4) 
being verified. 


The task of the present paper is, first, to show that solution (5) is valid only for stable linear objects and, 
second, to apply the frequency method, based on the Fourier transform, to the investigation of a broad class of 
objects (with many inputs and outputs and with constant lags), as well as to the analysis of multiloop systems 
with internal noise present. 


The solution of the first problem, with the detailed computations, is given in the Appendix to this paper, 
where we consider a stable linear object with one input and output and with constant lags. We, therefore, turn 
immediately to the consideration of the general case. 





imi BS 


af 








2. The Frequency Method of Determining the Dynamic Characteristics of Linear 
Objects with Many Inputs and Outputs and Constant Lags 








As is well known, to obtain the dynamic characteristics of an object with n inputs and m outputs, it is neces- 
sary to find nm impulsive responses or transfer functions [7]. It is easily shown that, to solve this problem, it suf- 
fices to consider an object with n inputs and one output. 


Indeed, by considering each output independently of the others, and by determining m times the dynamic 


characteristics of an object with n inputs and one output, we can successively determine all nm impulsive re- 
$ponses. 


Thus, in the general case, the problem reduces to the determination of the dynamic characteristics of ob- 
jects with n inputs and one output. To solve this problem, we make the following assumptions: 
a) the random processes at the objects’ inputs and outputs are ergodic and stationary; 


b) the objects under investigation are stable; 


c) the autocorrelation and cross-correlation functions of the processes at the input and output of an object 
are known. 


In addition, for greater generality, we shall assume that the object possesses constant lags with respect to 
each input, and that the lags with respect to different inputs are not equal to each other. 


We consider two methods of solving this problem. 


A. Reduction of the Problem to the Solution of a System of Linear Algebraic Equations 





For an object with n inputs and one output, the output signal z(t), by virtue of the principle of superposition, 
may be written in the form . 


oo oO 














z(t) =\ x, (t —9) k, (9) dd +... 4 \ Tn (t —9) kn (9) d®, (8) 
0 0 
where x;(t), . . . , X,(t) are the disturbances at the inputs (Fig. 1), 
a(t) 4y) sR ae kp (t) are the impulsive responses. 
Sf) | nf ~~ 
mete. > z(t) If we multiply both sides of Eq. (8), successively, by 
asa X3(t-T)/T, .. . , X,(t—7T)/T and integrate with respect to t we 
Zed [Met ————=-~ obtain, by virme of relationships (2) and (3), a system of n linear 
at) [wt -- integral equations relating the n impulsive responses: 
Fig. 1. 
ry > 
R,., (t) = \ Rx, (: — 9) k, (98) dd +... + \ Rx,x, (= —) k,, (9) d%, 
‘di ‘dn 
ts Re ak Sh EE (9) 
Rix, (2) = \ Ras (< — 9) k, (9) d9 +... .4+ \ Ry, (t — 9) ky (9) a9, 
‘di “an 


where td etdye ss tan are the constant lags with respect to the first, second, etc.,inputs respectively. 


In the particular case when x,(t), . - 


- » Xp (t) are statistically independent or, in other words, are uncor- 
related, i.e., when 


Rxixi fr) =0 (10) 
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for all r and for i # j, we obtain from system (9) 


Ras, (2) =\ Re, (¢— 9) hy (0) a8, 


Pi ect a Cr ee ee ae (11) 
Rix, (t) = \ Re, (t—9) kn (9) dd. 


Fe—3g 


Equations (11) are solved separately by the method presented in the Appendix. We, therefore, consider only 


the case when x;(t),..., Xp (t) are cross-correlated. Thus, we will seek a solution of system (9), k,(t),.. -, kp (t), 


which satisfies the condition of physical realizability: 


k, (t) = 0 (t < ty), 
ke (t) = 0 (t < @), 
2 eae see (12) 
k,(t) = 0 (<ty) 

where the constant lags, tay td. veee tgn are not known beforehand. 


We first find the solution of system (9) in the frequency domain. For this, by carrying out the computations 
analogous to those in the Appendix, we present system (9) in the frequency domain in the following way: 


Szx, (@) = Sx, (@) D, (iw) +... + Sx, x, (@) Dy (iw), 
Sra Et ead oe ee I i (13) 
Szx,, (@) = Sx,x, (@) D, (iw) +.. +55, (@) D, (i), 


where ©,(iw), ... , ®, (iw) are the transfer functions corresponding to impulsive responses k;(t),..., Kp (t), 
©; (iw) = 9; (iw)e-iutg; (j= 1,2,...,m) and (iw) is the transfer function corresponding to the idealized link 
without lags, 


+00 
S2x;(@) = \ eter R,., (t) dr, (14) 


-—oo 


a 
Sxix, (@) = \ e* Rex, (t) de, (15) 


Hoo 
Sx, (@) = \ eter R(t) dt. (16) 


—0o 


By solving system (13) with respect to 4;(iw), . . . , ®,(iw), we find the object's transfer functions: 
n stale 
(= 1)"F9S,., (0) Ans 


0; (iw) = = (17) 





A 


Sx,(@) . . . Sx,x,(@) 
eS ee a ies ewlne ae . 4 (18) 
Sx,x,(@) . - - Sx, (@) 








is 
ar 














is the determinant of system (13) and Ay; ‘s the minor obtained from the determinant by striking out its k-th row 
and j-th column. 


By then applying the Fourier transform, we obtain the solution in the time domain as well: 


1 Y seen 
k;() = = \ D; (iw) e#'dw (j = 1, 2, ..., A). (19) 
—oo 
By repeating the argument given in the Appendix, we easily see that the solution found satisfies condition 
(12) for physical realizability. The method presented requires significantly less work than the solution in the 
time domain based on the separation of the component of each input which is uncorrelated with the remaining 
input stimull. 


We turn now to the consideration of still another method of obtaining the dynamic characteristics of an ob- 
ject with n inputs and one output. 


B. Determination of Objects’ Dynamic Characteristics by the Method of Selecting Uncorrelated Stimull 





As was shown earlier, for uncorrelated stimuli x;(t), . . . , X,(t), the dynamic characteristics of the object 
are determined independently of the remaining n—1 inputs. It can be shown, however, that in the general case 
when the stimuli x(t), ... , X,(t) are cross-correlated, the problem of determining the dynamic characteristics 
can be reduced to the solution of equations analogous to Eq. (1). 


Indeed, let us have succeeded in choosing stimulus e; which is correlated with the i-th input and uncorre- 
lated with the remaining n—1 inputs. Then, multiplying Eq. (8) by e;(t—1) and integrating over t, we get, by 
virtue of relationship (3), that 


Rie, (t) = \ Ryje, (t — 9) x (9) d9 (ims 1, 2, .s 0) Bh (20) 
0 


This equation differs from Eq. (1) only in that the cross-correlation function Ry sey (T) appears in the inte- 
grand. Therefore, the solution of this equation in the frequency domain can be obtained analogously to that of 
Eq. (1) (cf. the Appendix) by a simple division of spectral densities. However, despite this simplicity, the method 
presented here, in comparison to the previous one, possesses an essential disadvantage. This disadvantage lies in 
the difficulty of choosing stimuli e;(t) which are not correlated with the x;(t) (i # j) and with the lack of assur- 
ance that they are, in fact, uncorrelated. In the conclusion, 
we shall consider a case when, from the physical conditions, 
one can comparatively simply choose the stimuli e; (t) and, 
consequently, apply the method just described. 


3. Statistical Analysis of Multiloop Control 
Systems with Internal Noise 








With the system shown in Fig, 2 as an example, we 
shall demonstrate the application of the method of choosing 
uncorrelated stimuli to the analysis of multiloop control sys- 
tems with internal noise. By assuming that the three assump- 
tions of the previous section hold, we find, for example, the 
dynamic characteristics of link 3. As is easily seen, the pres~- 
ence of feedback in the system leads to the result that the 
stimulus at the link's input, x(t), is correlated with the inter- 
nal noise fs(t), and relationship (1) does not hold. Due to 
this, it is advantageous, in determining the dynamic charac- 
teristics of the object, to use the method of choosing stimuli which are uncorrelated with the internal noise. In- 
deed, let us, in addition to writing the process x (t) and y(t) at the object's input and output, write the process E(t), 

















Fig. 2. 












































which is correlated with the input and, consequently, with the output stimuli, but is uncorrelated with the dis- 
turbances inside the object. Then, the object's dynamic characteristics are easily determined from the equation 


co 


Ryz (2) =\ Res (t«—9) k(9) a8, (21) 


considered in the previous section. 
For the quantity E(t), which is uncorrelated with the noise f(t), we can choose, based on the hypothesis of 


the independence of the internal noise and the external stimuli, any one of the stimuli e;(t), eg(t), or eg(t) (Fig. 2), 


SUMMARY 


The paper presented a frequency method for determining the dynamic characteristics of objects of auto- 
matic control from data obtained in their normal processes of use — autocorrelation and cross-correlation func~- 
tions ~ applicable to a broad class of stable linear objects — those with many inputs and outputs and with constant 
lags. The application of the method presented was shown for the analysis of stable elements of linear multiloop 
automatic control systems with internal noise. 


The results presented for the method of statistical linearization of nonlinearities can be easily extended to 
quasi-linear systems, 


APPENDIX 


We consider a stable linear object with one input and one output. For greater generality, we choose the 
case when the object has a constant lag ty. Such an object can be presented in the form of two series-connected 
links: an ideal link without lags 1 and a lagged link 2 (Fig. 3). 
We understand by a lagged link an element of a control system in 


settles 4. age me pag. we a * which the output quantity reproduces, without distortion, the varia- 
; a(t t a(t-ts) tions in the input quantity but with some fixed delay t, [11]. 
' ' 
bee ee 4 Thus, if k(t) for the ideal link without lag has the form 
Fig. 3 shown in Fig. 4a, then the introduction of the additional lag link 


leads to a shift of k(t) to the right by the magnitude of the delay 
(Fig. 4b). 


By virtue of this, the impulsive response for an object with a constant delay k,(t) can be given in the form 


kg(t)=k(t—t@), t>td 
kg (t) = 0, t<td (22) 


If expression (22) is taken into account, Eq. (1) is written in the following form: 





R,y (t) = \ R,, (t — 9) k (® — tg) dd. (23) 
t, 
aft) alth-a(t-4) 
‘2 
a 4 b 
Fig. 4. 











Thus, the determination of the dynamic characteristics of objects with constant lags from known correla - 


tion functions R, (rT) and Ryx (T) leads to the finding of the solution of Eq. (23) which satisfies the condition of 
physical realizability: 


k(t tg)= ky (t) = 0 for t <t,y (24) 


With this, the constant lag ty is assumed to be unknown beforehand. We now wn to the solution of Eq. (23). 
By multiplying both sides of Eq: (23) by e"“7 and integrating over Tr from — to +00, we obtain 


\ Menge 7 (+) dt = \ e Td \ R, (t—O)k (& — ty) ds. (25) 


If, in the last integrals, we make the change of variable A = r — 9, we get 


+00 co +co 
\ i (z)dt = \ e too, (& — ts) dd \ R, (A) e~ io a). (26) 
—oo ts —— 


By introducing the new variable y = 9~ty, we transform Eq. (26) in the following way: 


+00 oo -++0o 
\ eR, (t) dt = tal e '°?k (9) dp \ R, (A) eo **'da. (27) 
—oo 0 —oo 
By using the spectral densities 
-+-0o 
S, (@) = \ e—i@tR, . (t) dr, (28) 
-—0O 
— 
S,. (@) = = \ e tot R, (t)d (29) 
—oo 


we rewrite relationship (27) in the form 


e (30) 
Sy (w) = eA (iw) S, (o), 
where 
co 
® (iw) =| e~ °° (9) do 
(31) 
is the transfer function of the ideal link without lags. 
By denoting by 
® (iw) = ® (iw) ed (32) 
the transfer function of the object with a constant lag, we have 
Syx (@) (33) 





® (iw) == S.. (@) 
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If, as a particular case, we sett, = 0 in expression (33), we obtain, for the ideal link without lags 


Sys (@) 


et hn 34 
5. (o) (34) 


® (iw) = 


Thus, the real and imaginary frequency characteristics of a stable linear object are found by a simple di- 
vision of the mutual spectral density Syx (w) by the spectral density S, (w). 


If we again use the Fourier transform, we obtain the solution of Eq. (23) in the time domain: 


P +00 . ‘ +20 Sy. (@) , 
k(t)= ~y \ ® (iw) "dw = 35 \ rm all dow. (35) 
—o0 — * 


The integral in (35) can be comparatively simply computed, either by the theory of residues, for which it 
is necessary to go from the frequency characteristic given in some frequency range to the analytical expression 
for the transfer function, or by a graphicoanalytic method [2]. 








awe yg In the second case, it is not necessary to find the analytical expression for the 
transfer function beforehand. It suffices to limit oneself to the graphic presentation of 
one real frequency characteristic. 
o|_ —* We now show that the function k(t) just found satisfies condition (24). For this, 
we turn to the computation of k(t) for t < tg, keeping in mind that, by virme of the ob- 
ject's stability, all poles of the function lie in the left half-plane, i.e., that the function 
(s) is analytic and bounded in the right half-plane. 
le By presenting the expression found for k(t) in the form 
Fig. 5. ‘ “+00 
k(t) = >= \ ® (iw) e'—-Bdw, (36) 
—oo 
we consider the integral 
a | @(s) e'-Ha 37 
2ni (s)e das, (37) 
Cr 


along semicircle C* with radius r, and lying in the right half-plane (Fig. 5). 


For t—tg< 0 and for |r| + a, integral (20) tends to zero [12]. We may, therefore, write 


4 +ico ‘ +-ico ‘Se 
k (t) = 3 \ ®(s) e°""@ds = 5 \ ® (s) Ws + ==> \ ® (s) e®—'dds (38) 
or 
4 
k(t)= xr? @(she"@ds = for t— 1d <0, (39) 


where the integration is carried out along the closed contour consisting of the imaginary axis and the semicircle 
of infinite radius lying in the right half-plane. 


According to the theory of residues, the value of this integral equals the sum of the residues of the inte- 
grand for all poles lying within the contour, multiplied by 27i. However, by virtue of the system's stability, 
function @(s) does not contain any poles in the right half-plane. Therefore, 
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-+ ico 


1 
k (t) = aT \ @ (s) et ds = 0 


Be for t—tg<0 (40) 
—to 
or 
i 
k(t) = > \ Dio) edw=0 for t—W<0. (41) 
—oo 


Thus, for stable linear objects, the solution of Eq. (23), which satisfies condition (24) for physical realiza - 
bility, is expressed by formulas (33) and (35). The constant lag ty is determined after the impulsive response is 
computed (Fig. 4b) from the graph of k,(t). 


In the proof just given, we assumed that S, (w) and Syx (w) were exact spectral densities. Then, by virtue 
of the relationship 


1 . 
S_., (w) ia Y,7 (iw) Xp (iw) Y, (io) 
xy . T 
AY 1 ° 





where 
T 


T 
X_7 (iw) = \ x (t)e**at, Y,7 (i) = \ y (t) e—tot ge, 
0 0 


the function Syx (w)/S, (w) for a stable object will be analytic and bounded in the lower w half-plane, since 
7 (lw) is factored out, and the corresponding singularities drop out. 


In actuality, S,(w) and Syx (w) are computed with some error, and the zeros of S,(w) may not exactly co- 
incide with the zeros of Syx (w) in the lower half-plane, i.e., Syx(w)/Sx(w) may have poles in the lower w half- 
plane. This latter circumstance can be reflected in the results of computing (iw) by formula (33). 


Indeed, in the case when the object under investigation has a yery small margin of stability, at the points 
on the w axis close to the poles of w)/S,(w), the computed frequency characteristic may have sharp jumps. 
Such values of (iw) should be either rendered more exact by varying the path of integration in expression (25), 
or should be ignored completely. 
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ON THE SYNTHESIS OF LINEAR AUTOMATIC CONTROL 
SYSTEMS WITH VARIABLE PARAMETERS 


S.V. Mal'chikoy 


(Moscow) 


A method is presented for determining, from the optimal weight function of an 
entire system and the weight functions of the known links in the system, the optimal 
weight function of the system's correcting link, which may be placed either onthe for- 
ward path or in a feedback path. The method presented here allows one to circum- 
vent the difficulties associated with the necessity of solving Volterra integral equations 
of the first type. 


In designing an automatic control system which will operate under noisy conditions, one must know its op- 
timal dynamic characteristics, i.e., those which will provide the greatest possible accuracy of system operation 
under the given conditions. The recently developed theory of optimal systems allows one to determine the op- 
timal dynamic characteristics of an automatic control system if one knows the statistical characteristics of the 
input signal and noise. Once one knows the system's optimal dynamic characteristics, one can try so to choose 
the design schematic and the parameters that the dynamic characteristics thus obtained will be close to the op- 
timal. There is still no general solution to this problem of automatic control system synthesis. In the present 
paper, the attempt is made to cite one of the possible ways of solving the problem of synthesizing a linear auto- 
matic control system if, as the system's optimal dynamic characteristic, its weight function W (t,r) is given. 


A designer is ordinarily given part of the design schematic and some of the parameters of the system to be 
designed, variations of which are either impossible or, for some reason or other, undesirable. In order to provide 
the required dynamic characteristics of the system in such a case, one can include in it some correcting links. 

It is possible to determine the optimal weight functions of the connecting links from the given optimal weight 
function of the system, found by means of optimal system theory, and from the weight functions of the given links 
of the system. Knowing the weight functions of the correcting links, one can find the differential equations by 
which they must be described [1,3]. The differential equations which describe the uperation of the correcting 
links completely determine their structure and parameters, which can be found for each concrete form of differ- 
ential equation. The questions involved in designing the schematic for the 
correcting links from their differential equations are not considered in this 
paper. 


In order to determine the optimal weight functions of the correcting 
links, one must know the relationships between the weight functions of the 
individual links and their interconnections. Any automatic control system 
consists of three types of link interconnections; parallel (Fig. la), series 
(Fig. 1b), and with a feedback path (Fig. 1c). If we know, therefore, for 
each of these three types of connection, how to find the weight function of 
one link from the weight function of the entire connection and the weight 
function of the other link in the connection, then we can always find the 
weight function of a correcting link which enters into a system of any de- 
gree of complexity. 
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The relationships between the links’ weight functions and the weight functions of their standard intercon - 
nections can be derived in the following way. For a parallel connection of the links we have, as an obvious con- 
sequence of the system's linearity, 


Wie (t, t) 7 W, (t, t) + Ww, (t, t), (1) 


where the subscripts 1 and 2 denote the link to which the corresponding weight function cosine and the sub- 
script 12 denotes the connection of links 1 and 2 (Fig. 1). 


By applying a delta -function signal to the input of series-connected links (Fig. 1b), we must obtain, at the 
output of link 1, the weight function W,(t,r) of this link and, at the output of link 2, the weight function W»(t,r) 
of the connection. If we now assume that W,(t,1r) is the input signal to link 2, and take into account that, from 
the condition of physical realizability of the system, 


W, (t,t) =040dW, (t,t) =0 for t <x, 
we obtain 


Waa (t,t) =) Wy (21,2) Walt, 1) de. (2) 


Finally, if we apply a delta-function signal to the input of a connection with feedback (Fig. 1c), we obtain 
the weight function W»(t,r) at the output. Then, in accordance with (2), at the output of link 2 we get 


™%: 


{ W 49 (ts, t) Wa (t,t) ty, at the input of link 1 we have a signal of the form 


t 
t 


8 (ty — )—{ W 43 (te, t) Wo (11, t2) dt2, and, at the output of link 1, we again have the weight function W 4(t.r), Le, 


t %; 


W 32 (t, 7) = \ [3 (t; —*t)— \ W432 (t2, t) We (t1, te) des] W, (é, 1) dry. 


‘ t 
If we keep in mind that 3 (t; —t) W, (t, t1) dt, = W, (t,t), we obtain, finally, 


Tv! 


t 
Wis (t,t) = W, (t,t) m- W;, (t, 1) dt, j W 32 (t2, t) Wo (tr, ta) ate. (3) 
If there is an amplifying link with unit gain in the feedback path, then W,(r4,T:) = 5(7,—T 2), and formula 
(3) takes the form: 


Was (t,t) = W(t, 2) —) Wi (t, 1) Wie (41, 2) de. (4) 


It follows, from a consideration of (1), (2), amd (3) that only in the case of a parallel connection of the 
links will the weight function of any link in this connection be easily determined from the weight function of the 
entire connection and the weight function of the other link. In all other cases, the desired weight function will 
enter into an integrand and, consequently, it is necessary to solve an integral equation in order to determine it. 
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If one is considering a connection with feedback, and it is necessary to find the weight function of link 1 
(Fig. 1c) from the given weight function of the entire connection and the weight function of link 2 then, for its 
determination, it is necessary to solve, for W;(t,7), integral Eq. (3), which is a Volterra integral equation of the 
second type with kemel 


as | 


\ W 2 (t2, t) Wo (ts, ta) ate. 


t 


The solution of a Volterra integral equation of the second type presents no serious difficulties, and may al- 
ways be carried out by the method of successive approximations [2]. 


If the weight function of link 1 is given, and it is required to determine the weight function of link 2, which 
is in the feedback path, then it is necessary first to solve Eq. (3) for the integral \ W 12 (t2,t)Wo(t1,t2)dt2, 


and then, knowing the value of this integral as a function of r and ry, to determine Wg(rT 4,7). For this it is 
necessary both times to solve a Volterra integral equation of the first type. The problem of finding the weight 
function of either link in a series connection from the weight function of the other link and the weight function 
of the connection also leads to the solving of a Volterra integral equation of the first type. The solution of Vol- 
terra integral equations of the first type is a very arduous task which can be solved with sufficient speed and ac- 
curacy only by computers. In addition, Volterra integral equations of the first type may, in general, not have 
solutions. 





+47} Las. arated 
» rer ee 


Se 











Fig. 2. Fig. 3. 


Let us attempt to avoid the difficulties connected with the solution of Volterra integral equations of the 
first type by using the concept of inverse links. Two links will be called inverses of each other if, when con- 
nected in series, they effect the identity transformation on any signal applied to the input of this connection, L.e., 
an inverse link carries out a transformation of its input signal which is inverse to the transformation effected by 
the direct link on its input signal (Fig. 2). A link which is inverse to a given link will be denoted by the same 
symbol, but with a superscript minus sign. It follows from Fig. 2 that it is possible to interchange inverse links 
which are series-connected. Since the weight function of the series connection of inverse links equals 6(t~Tr), 
the following relationships hold between the weight functions of mutually inverse links: 


t 
|W (x, t) W~ (t, =) dt, = 6(¢ —*), (5) 


t 


t 
\w- (t,, t) W (t, = )dt, = 8(¢ —*). (6) 


t 


Suppose now that we have a series connection of links (Fig. 1b) and we know the weight function W»(t,T) 
of this connection and the weight function W,(t,r) of link 2. It is required to determine the weight function 
W,(t.7) of link 1. We add, to the given series connection, a link inverse to link 2 (Fig. 3a). The weight func- 
tion of the connection thus obtained equals the weight function of link 1. On the basis of formula (2), we get 
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t 
Wi (t,t) =) Wis (1, 2) Wa (tad. (7) 


tT 


Thus, the solution of a Volterra integral equation of the first type was reduced to a simple quadrature. 
The weight function Wg(t,r) of the inverse link is easily determined from the differential equation which de- 
scribes the operation of link 2. The differential equations of known links are ordinarily given. 


If link 1 is known, and it is necessary to determine the weight function of link 2, then the given series con- 
nection should be connected to the output of a link inverse to link 1 (Fig. 3b). We then obtain 


t 


W(t, 2) =\ Wi (1,2) Walt, 1) dey. (8) 


t 


Let it now be necessary for us to determine the weight function of link 2, located in the feedback path 
(Fig. 1c), for the given weight function of the entire connection and the weight function of link 1. We show now 
that, by using the concept of inverse links, we can reduce this problem to the solution of a Volterra integral 
equation of the second type. For this, we find the weight function of the link which is inverse to the entire con- 
figuration shown in Fig. Ic. We rewrite formula (3), which defines the weight function of a feedback connection, 
in the form 


t T%: 


W, (t,t) — Wy, (t,t) = \ W(t, 1) dey \ W 40 (to, t) Wo (41, te) dtp. (9) 


t t 


To analyze this formula, we introduce the following notation: 


™1 


Wa(t,t)= \ W 12 (2, t) We (41, te) dts, (10) 
W, (t,t) = W, (t,t) — Wy, (tt, *). (11) 


With this notation, formula (9) takes the following form: 


We (t,t) =| Wa (2) Wilt.) dey. (12) 


t 


It follows from formula (12) that Wp(t,r) is the weight function of a series connection of links with weight 
functions Wa(t,r) and W,(t,7), respectively. In its turn, W,(t,r) is, according to (10), the weight function of a 
series connection of links with weight functions Wy»(t,r) and W,(t,r). Consequently, formula (9) expresses the 
fact that the weight function of a series connection of three links with weight functions W»(t,r), W,(t,7r), and 
W,(t,r), respectively equals W,(t.r) = W,(t,T) — W(t, 7) (Fig. 4). 


The weight function of link A, enclosed in dotted lines on Fig. 4, is defimed by formula (10). On the other 
hand, it can be determined by the method of inverse links from the weight function of the entire configuration 
and the weight function of link 1: 


‘ t 
Wa (t,t) = \Ws (t, ,t)Wy(t, 1) dt, ={ [W, (t,t) — Was (11, ) |W; (t,t) dt, = 


t 3 (13) 
= 8(t—1)—) Waa (n1, 2) Wilt, ty) dry. 


7 


By comparing (10) and (13) we find, after some elementary transformations, that 
























W 32 (71, 7) |W. (t,t) + W, (t, *:)| dt, = 8(t —t), 


from whence, on the basis of (5), it follows that 
Wo (t,t) = We (t,t) + W, (t, 2. (15) 


From this formula, Wg(t,r) can be easily determined if the inverse weight function Wg(t,7r) is known or 
can be easily determined. In practice, it can be determined only in the rare cases when the differential equation 
describing the operation of the whole configuration is given or can be easily found. We now show how to deter- 
mine weight function W2(t,T) if Wg(t,7) is not known. 
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Fig. 5. Fig. 6. 


In accordance with (15), the block schematic of the circuit (link) which is inverse to the link shunted by a 
feedback path is represented as is shown in Fig. 5. 


Since, in the block schematic obtained, links 2 and 1” are completely equivalent, the block schematic of 
the direct link can be represented as shown in Fig. 6. 


In the forward circuit of the feedback system just obtained, which is equivalent to the system shown in 
Fig. 1c, there is a link inverse to link 2. The weight function of this link is found from an equation of the form 


t 


Wie (t, t) - W, (Z, t) 4 \ WwW; (t, t,)dt, \ Wis (Te, t) W, (t, T?) dt, (16) 


which is a Volterra integral equation of the second type for Wg(t.T). 


The method just given allows one to determine the weight function of a link which is inverse to the link 
in the feedback path. For further solution of the synthesis problem, i.e., for finding the optimal weight functions 
of the simplest links, and also for finding the differential equation which describes the operation of the link, it is 
completely indifferent with which weight function one operates, whether with the direct or with the inverse. 


It should also be mentioned that the method given here for finding the weight functions of links from given 
weight functions of their interconnections proves the existence of the solutions of Volterra integral equations of 
the first type of the form (2) or (3) for finding the weight function of a link in a feedback path. 


EXAMPLES 
1. Let there be given a series connection of links 1 and 2 whose weight function is 
Wy, (t, t) = 1(t — t) A(t) eC), 


where link 1 is described by the differential equation 
Tx +a2=K(t)y. 














We now determine the weight function of link 2. 
The weight function of the link inverse to link 1 equals 


i T | 
Wi = RT 1) + RG Pt — 2». 


From formula (8) we obtain 


t 
+ es 4 5 

Wate VE Eh d+ RE_ra—D]L = AGE Vae. = 
t 


t 














t 
a \ 1 (t — 7) eae oe? (lt) 8! (ty — t) dt + 1 (¢ —t) mc et — Wg x, — 
= [100 APA] su vf oe. 
=i 2) Eee e* (=) _4 (t — 2) [r= Fact e* U~*) 5. ho he—* «-9] + 
+4(t—7) Sen eA Ul—*) 
=a! (¢— 2) + 1(0¢—2) oe ae a rh? (z2)| ght), 


2. We now consider a feedback connection (Fig. 1c) whose weight function is 
Wis (t, t) = 1 (t —t) A(t) eC), 
where link 1 is described by the differential equation 
Tx, + 2, = K (t) y. 


We determine the weight function of link 2. 


From formula (15) we get 
Ws (t, t) = W,, (t, t) — Wy (t, =), 


where 
- T we 1 
W, = Ki a’ (i —t) + Ky &(i—7). 


In the given case, the differential equation which describes the operation of the connection as a whole is 
easily determined: 


he A(t) 
Ld dy, tess 5 
Whence, 


- 1 dh 
Wi.) = ZH (¢ — 7) + AG ° ¢—*). 


Consequently, 





; i v1. i 1 
wi. =[Am -xa]* sy [am * Ray |b. 








3. We now consider the feedback connection (Fig. 1c) whose weight function is 
t 
Wis (t, t) == { (t— t) A(t) exp (— \ B (t) at), 


and where link 1 is described by the differential equation 


a, (t)zy + dy (t) 2, = y. 
| 


We determine weight function W,(t,r) without first determining Wy(t,T). It is obvious that 


Wy (t, T) = ay (t) 8° (4 — t) + ay (t) 8 (t — 7). 
From Eq. (16) we obtain 


we (t, t) = 1(i — t) A(t) exp ie B (t) dt ) a 


4A" oO = 


™%; t 


Writ, 1) dt \ 1 (t; — 1) A (+) exp (— \ B (t) at) fay (41) 8° (+; — 49) + 


t 


4 


re 


+ ay (7) t (T — T)] dt. 


After some elementary transformations we get 


A(t) 4 
Ait) ai(t) exp(—| Bear) + 





W, (t,t) =1 (¢—*) >= 


A(t) __¢ 
+ = (=) ay (t) \ Wy (t, 71) [4 (t1) == Gp (71) B (t:)] exp (— \B (t) dt) dt. 


t 





This equation is solved by the method of successive approximations. As the first approximation we take 





- A(t . 
Ww, (t,t) = 1(¢— 7) veto me (=) exp(—| B at) 


By substiniting the first approximation W;(t,7) under the integral sign in the original integral equation we 
obtain the second approximation. The process of successive approximations continues until the last approxima - 
tion does not differ appreciably from the previous one. 
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SYNTHESIZING THE ELEMENTS OF LINEAR AUTOMATIC 
CONTROL SYSTEMS 


I.A. Orurk 
(Leningrad) 


A method is proposed for synthesizing the elements of linear automatic control 
systems. The method is based on the use of time characteristics and relationships ob- 
| tained from the system's (and links’) integral equations with kernels in the form of 

1 polynomials, in conjunction with the D-decomposition of parameter space. 


The method can be used, in particular, for programing the synthesis problem for 
computer solution. 


For the construction of an automatic control system with given dynamic characteristics, it is required to 
determine the individual parameters of various links, and to find completely certain links of the system, 


A number of engineering methods have been developed for synthesizing the elements of control systems 
[1-6], but it is impossible to consider the problem as having been completely solved. 


The method proposed in the present work for synthesizing linear systems with constant parameters, this 
method being based on the use of time characteristics (transient responses) and Volterra integral equations with 


polynomial kernels in conjunction with the D-decomposition of parameter space, allows one to make broad use 
of computers, 


1. Basic Relationships 


One can, from the coordinate representation of a system, obtain an integral equation in this coordinate 
. where the kernel is in the form of a polynomial (6). Let the transient response in the system under investigation 
be known: 


Dan + Oy —yP +++ byp™— 1+ b0p™ 


2) XP) ag yh bed pe aap 





(1) 


where m =n. 


We divide both sides of the equation by p": 


7 (Qn + Qn—1Pp +2°* -+ agp) =P = (bmn + bm ap + ° + + bop”). (2) 


We then introduce the notation 





2) _y (p), @ 
P 




















which, when we transfer back to the time domain, gives 


t 


++ (a(dt...dt 


0 


— 7 a hd 


z()=y"@, 9) =\ 





(4) 


For the auxiliary function y(t), as follows from expression (4), there are zero initial conditions, y) = yy = 


De ioe yo) = 0, since for t= 0, the function y(t) and its n—1 derivatives vanish. 


By substituting expression (3) in (2) and returning to the time domain, we obtain 


any (t) 4 Any’ (t) ile 1 ay” (t) _ 


1” {1 {*- m 


= bm A Din fail soot O 


9 (n—m)! * 
We then write the Maclaurin expansion for function y (t) 


t 
Pe Teter = (¢—s)*" 
y(t) = yy tut +--+) ym (pds =| 2 (ods. 


0 0 


(5) 


(6) 


If we differentiate Eq. (6) n times and substitute the expressions for y(t), y‘(t),..., etc., thus obtained in 


Eq. (5), we obtain a Volterra integral equation of the second type in x(t): 


a,x (t) — \ K (t, s)(s)ds = N (0). 


0 


Here, the kernel is 


and the free term is 


“ap pr—mtk 


N (t) = >) be amare. 


k=-0 


(7) 


(8) 


(9) 


The integral equation with the kernel given in (8) occurs in the literature [7]. There, however, the equa- 


tion is in a higher derivative, and not in the coordinate x(t) itself, as is the case in the present paper. 


The kernel in (8) and the free term in (9) of Eq. (7) are linear functions of the representation coefficients 


aj, and by. 


The definite integral in expression (7) can be computed approximately by using one of the quadramre for- 
mulas. If we take d as the integration increment and then apply the trapezoid formula for calculating the in- 


tegral, we obtain the well-known expression for the series of successive ordinates of curve x(t) [4]. 


For the first ordinate, 


, % .- 
Nitd> Ks 


d 
rus 





(10) 














































For any moment of time ty = vd, 


N,. 4 a(5 K, +2:K,_,+.-.-+2, Ax) 
rp — ———- 2 —______,-_—____ (11) 
a + 45 
where 
, : ~ ee 
Ly = 2 (ly), Ny = Nite) nw Ky= —> Qk Te—Apt (12) 
k=1 


The error of this method is determined by the error introduced by the trapezoid formula in expression (11), 
In order to decrease this error one can, instead of the trapezoid formula, use Laplace's quadrature formula. 


Equations (7) and (11) establish linear relationships between the coefficients of the numerator and denomi- 
nator of representation (1) and the particular values of the coordinate x(t), and are used below for the purpose of 
synthesis. For given numerical values of ay, and by, the same expressions (and somewhat transformed expressions) 
can be used for finding successively particular values of coordinate x(t) (cf. [6)). 


2. Determining the System's Parameters 


Initially, we consider the synthesis problem as thus posed: the structural scheme of the system {s given; 
defined are several parameters of one or several of the system's links (gains, time constants, etc.). Let the repre- 
sentation of coordinate x(t) be given by expression (1), where a), and by, depend linearly on the system parameters 

to be determined. 


of + fio The parameters must be determined so as to provide: 
" 1) a given degree of stability and oscillation of the system; 


' 

! 

3 ‘af 2) an accepted form of the transient response curve for the coordinate 
‘ x(t), a minimum and maximum value of the transient response, a given 





e p=brjur velocity and duration of the transient response. 


a é Point 1 is attained by introducing into the calculations the conditions 
ie ; , which determine the placement of the roots of the characteristic equation 
» all Bed inside the given region of the complex p plane acnpbu'n'a (Fig. 1), called 
s* \ +3 the domain of given quality in the sequel, for which one uses, as suggested 
/ by Yu.L. Neimark, the D-decomposition of the plane of one or two param- 
A eters [5]. 








Point 2 is attained by moving the points of the actual transient re- 
' sponse curve of x(t) close to the given curve x"(t), for which one uses ex- 
lai ;, pressions (9), (11), and (12). 


a 








Several parameters being sought may be determined or eliminated 

/ from the aforementioned equations by using the well-known relationships 

4 of the operational method, providing coincidence of the actual process x (t) 
with the given x(t) for t = 0 and t = o: 





Fig. 1. 


b 
| x (0) = x° (0) =X (coy =, x (00) = 29 (00) = X (0) =—. (13) 
In accordance with the methodology being considered here, the remaining q parameters of the system are 
determined in the following way. 


1, We set up q—2 (or q-1) equations from formula (11) for q—-2 (q—1) points lying on the curve of the 
given process (Fig. 2): 
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d \ 3 ra : 
#y (My 1 ay +) N, -!- d(- Ry, -] PR. teow | Ks 1 oy (14) 


These equations are solved for the q-2 (q—1) parameters, i.e., we obtain equations of the form 
31 = fy (Sq, Sq), Fe = fo (Sq, Sq -a)y---, Fq-2 = fq (Sy, 74-1) (15) 


The expressions obtained are substituted in the system's characteristic equation D(p) in order to eliminate 
part of the unknowns. As a result, we obtain an equation with two (or one) parameters Oq Fq-1 


D(p, Sq, Sq—1) = An + Qn—1P +... + agp". (16) 


2. If we substitute p = yell in Eq. (16) and then so vary x and y that the end of vector p in the p-plane 
slides along the boundary of the domain of given quantity (Fig. 1), we obtain a representation of this domain on 
the plane of parameters 0g, Oq-1. 


3. Inside the plane of given quality (on the parameter plane) we choose a point, and from these values of 
Og and Gq-, we determine the remaining parameters by formulas (15). 


The determination of the parameters by the representation method guarantees system stability as well as 
the other requirements enumerated above on control quality. 


The most advantageous form of the region of given quality is the six-sided one (Fig. 1) which gives the 
limiting value of the real roots of the characteristic equation and of the system's frequency and oscillation. 
Thanks to the existence of well-defined relationships between the parameters, introduced by Eqs, (14), the cutting 
off of region » ng (Fig. 1), as a rule, is insignificant. Finally, the parameters 
of the given curve x*(t) must be matched with the parameters of the Ry, Re, and 
a region. 


If high accuracy in reproducing the given process x*(t) is not required, 
one can simplify the computations by representing, not the domain of given 
quality, but the boundary of the region of stability (the imaginary axis of the 
p-plane). An even greater simplification occurs in the case when the system 
is stable for all practically possible values of the parameters to be determined. 
There is then no necessity of using the D-decomposition: all q parameters are 
found as a result of solving the system (14) of g equations, followed by a test for stability of the system. 





It should be emphasized that the method allows one to determine a large number of parameters. In the 
Appendices, we give examples of synthesis with four and with three parameters to be varied. 


3. On the Synthesis of Individual Links 


- By the synthesis of a link we understand the determination of its transfer function, its form (circuit) and its 
parameters. Without carrying out a detailed consideration here, we state that the method presented can also be 
used for the synthesis of individual links, including correcting links. The controller as a whole can be considered 
as an individual link; in this case, one determines the control law, the gain, etc. Links can be determined which 
are connected in complicated ways (not necessarily in series or in paraliel) with the remaining portions of the 
system. 


We first find the desirable transfer function and form of the link. Let a dynamic system be described by u 
equations with u generalized coordinates and definite disturbing stimuli. The last equation of the system, the 
equation of the link desired, has the form: 


‘ Tr 
Q(p)ar = Rp) tq X (p—) = = FB, (17) 


where Xp and x; are, respectively, the link's output and input coordinates and X (p) is its transfer function. 


If we are given the desired process for one of the coordinates x), we solve, in operator form, the first u-1 











~ rR engine 
















equations for x, and x and, after substituting in Eq. (17), we find the desired link’s transfer function x%p). 


An exact reproduction of X%p) is generally disadvantageous, due to its complexity; we therefore deal with 
a compromise solution which permits the quality condition to hold approximately for the simplified link circuit, 
If the degree of the numerator of X°%(p) is greater than the degree of the denominator then, by separating the in- 
tegral part, we obtain terms of the form Ag and Aypy which define the law of transformation of the input signal, 
and which are reproduced, with a known error, by differentiating links. If the representation of X°(p) is a proper 
fraction, then we first determine the characteristic in the time domain, x°(t). From this, using the condition of 
approximate coincidence of the time domain characteristics x*(t) and Xapp(t), we determine the order (form) of 
the transfer function of the simplified link Xapp(P) = Xapp(t). To alleviate the choice of Xapp(p), one can set 
up 4 “dictionary” of standard links and their transfer functions. 


The coefficients of the numerator and denominator of Xapp(p) are determined by using relationships (13) 
and (14) and the condition that one obtain the best point approximation of the actual process x,,(t) to the given 
one when link Xapp(p) is connected in the circuit. For this, the condition of stability of the entire system must 
be met, where the system includes the link being determined. This latter requirement is met by constructing the 
region of stability (or of given quality) as was shown in section 2. 


To simplify the calculations, Eqs. (14) can also be set up from the condition that Xapp(t) approximates to 
x(t). In this latter case, the error of reproducing Xy(t) will exceed the error in the approximation of Xq pp{t) to 
x(t). For the first error to lie within the limits of 10 to 12%, the second error, as a rule, must not exceed 3-4% 
with this, of course, stability of the whole system must be guaranteed. 


The circuit implementation of the simplified link with transfer function Xa pp(P) is carried out either im- 
mediately, or by application of the ordinary methods of circuit theory. 


Finally, we draw the following conclusions as to the basic advantages of the method of synthesis considered 
here. 


1, In using the method, one obtains linear relationships between the representation coefficients a} and by 
(in contradistinction to the method of integral estimates, for example) which makes it possible (in all cases when 
there is a linear relationship between the a, and by; and the parameters being sought) to determine any number 
of parameters. More than that, the accuracy of reproduction of a desirable process in the system increases as the 
number of parameters to be varied increases, with a relatively small increase in the amount of computational 
work required, Existing methods ordinarily allow one to determine no more than two parameters of the different 
links in a system. 


2. For the parameters found by the synthesis method presented here, it is guaranteed that one obtain a stable 
system and the given quality indicators, defined both by the representation’s denominator (domain of given qual- 
ity) and by its numerator (point approximation to the desired transient response curve). 


3. The parameters sought in the synthesis method considered here can appertain to different links of the 
system, One can find the form (circuit) and parameters of an individual link whose switching in determines the 
course of the process with given indications of quality. A link can be connected in a complex way (not neces- 
sarily in series or in parallel) with the rest of the system. 


4, The computations in the method presented here are carried out in accordance with a series of rules (cf. 
section 2) which are based on the execution of simple mathematical operations: calculating the coefficients of 
polynomials, solving systerns of linear algebraic equations and calculating the values of polynomials for certain 
values of a variable w. Thanks to its lack of complexity, the method does not require special preparation for 
the computer. Since the course of the calculations is formulated as a series of algorithms, it can conveniently 
be turned over to a computing machine. The method can also be used for programing the synthesis problem for 
an electronic computer. 


APPENDIX 1 


We now determine the four gains in the regulation of the excitation of a synchronous generator by the 
derivatives of current and yoltage. We determine the coefficients from a consideration of the transient response 
and the stability of small oscillations when the generator operates on a load via a long transmission line. The 
representation of the angle 6 between the vectors of the generating station emf and the voltage of the receiving 


system is given by expression (1), where: 
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n == 4, m == 2, 

by =: 0,287 — 0,389 By — 0,632 yo, 

by = 4,86 — 0,389 B, — 0.632 y,, 

bs = 111 — 0.389 B, — 0.632 yo, 

ag == (1,21—1,63 By — 2.65 2) 10-2, 

ay — (20,5—-1, 63 By — 0.124 Bg -+ 2.65 7, -+ 0 201 ye) 10°, wr 
ay =~ (18.7—1. 63 By — 0.124 By -+- 29.3 By — 2,65 yo — 1201 y,— 39.9 yo) 10-2, 

aq -= 2.10—0.124 - 10-* By + 0,293 B,; — 0.201 - 10-% y, — 0,399 y34, 

a, == 0.298 4 0.293 By — 0.399 yp. 


Here, 9 and By are the gains when control is by deviations of the voltage and the current from their nom- 
inal values. By calculating them from the given static characteristics, we find that By = 1.0 and 7, = — 15.0. 


To determine the gains when control is by the first and second derivatives of current and voltage, denoted, 
respectively, by By, Bg, y, and yg, we start out from the following requirements: the system must be stable and 
have a degree of stability of Ry = 2; the transient response must be monotonic (or with small overshoot) and prac- 
tically without oscillations, the duration of the transient response is T ~ 2.5 sec, the maximum velocity of the 


transient response is y ~ 2,5 sec™* 


To satisfy these requirements, it is necessary that the following conditions hold. 


1, The roots of the characteristic equation must lie inside the domain of given quality (Fig. 1) for which 
Ry = 2, a = 120°, R, = 12 and H= 10, The values of R, and H limit the real parts of the roots of the characteris- 
tic equation and the frequency of oscillation, and are arbitrarily taken here; a is determined by the oscillation 
of the system. 


2. The curve of the transient response must be close to the exponential curve with Eq. (1.55) (1—e*%y, 
giving the values X» = 0, xy = 0.62, xg = 1.1 and x(oo) = 1.55 for d= 0.25 sec. This also guarantees that the 
maximum velocity of the transient response is v ~ (x;—X,)/d = 2.5 sec™', The duration of the transient response 
is determined both by the given curve and by the magnitude of Ry. The point t, = 2d of the approximation to 
the desired curve is adopted from the condition that t, = Tmax, Where Tmax = 1/R, is the maximum time con- 
stant component of the transient response, With this, a small overshoot is possible for t > Trpyax- To avoid such 
a possibility, one can increase the number of points of the approximation. 




















Fig. 3. 


The calculations are carried out in the following order. 


The expressions for Ny and K, for t = d and the expression for N, for t = 2d are found from formulas (9) and 
(12); 


Ny = (2.44—0.101 By — 1.22 B, — 0.165 y, — 0.198 ys) 10-*, 
K, = (43.1—0.751 Ba -} 7.49 Be _ 3.95 vi —_— 10.17 2) 10- 3 
Nz = (17.69—0.811 B, — 4.86 By — 1.32 y; — 7.90 yq) 10-*. 
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For moments of time t, and t,, one sets up system (14) of equations and, after substituting the expressions 
for Ny, Nz and Ky in them, one solves the system for y; and y2: 


11 = — 0.949 3; — 2,67 Be + 7.76, v2 = — 0.0435 8, — 0.964 3. + 1.26. (19) 


The expressions obtained for y; and yz, as well as the values of y» and 8», are substituted in formulas (18) 
for the coefficients a), and these latter are then substituted in the characteristic polynomial given in (16). 


After this, one maps the domain of given quality (a six-sided figure) in the p-plane on the plane of param- 
eters 6,,8, by the method of D-decomposition, as shown on Fig. 3. Points a and b of the hexagon on the param- 
eter plane correspond to lines a and b, while the remaining mapped and image points have the subscripts, re- 
spectively, of c,d, .. .,k. 


In Fig. 8 is constructed also the image of a trapezoid in the plane p (Fig. 1, points f, g» h). As is evident, re- 
gion ngy is ignorable since in the plane 6;, B; it is not a working region. 


Since the contender, in the domain of given quality in the 8;,8, plane, has marks on the four roots which 
are greater than the minimum marks (zero root) (Fig. 3) and since the characteristic equation is of fourth degree, 
the marked domain 4 (Fig. 3) is the region of given quality. Approximately in its center we choose a point 

(marked by the cross) and then obtain the values of the 
Bi. parameters 8; = 5.8 and 6, = 2.8. Then, by Eqs. (19), we 















































ae , arr 

x(t) 2%(t) z(t) find that y; = — 5.23 and y, = — 1.69. 
16 => 

“Fevt) | With this, the synthesis of the system, i.e., the de- 
aad a eae termination of the gains, terminates. 
i , 
# 7 Beas To verify what we have done, we substitute the values 
f 25 7) 5 s 7 — found for 6;, Bs, yy and yz in the coefficients a; of the 
* sec 


characteristic equation and then compute the roots. We 
Fig. 4. obtain two pairs of complex conjugate roots: pj, = 
= — 2.014 1.18 j and ha? = 8.98 + 5.0 j. 


These results show that the roots lie within the given region, i.e., the hexagon. From the roots which we 
now know we find the exact value of the time-domain function x(t) = 5(*), and we then construct the curve of 
the transient response. 


Figure 4 shows the constructed curve for x(t) and, by the dashed line, shows the given curve x°(t), on which 
the points of approximation, x, and x,, are marked by crosses. In the interval 0 < t < 2d the two curves prac- 
tically coincide. For t > 2d, curve x(t) has a small overshoot, one within the 4% limits. Thus, the actual curve 
differs little from the given, optimal one. 


APPENDIX 2 


We now consider the same system as in the previous example, but with control of excitation by deviations 
from nominal values of the current and voltage and by angular derivatives, with gains o}, which are to be deter- 
mined (k = 1, 2, 3 is the order of the derivative by which control is carried out). 


The representation 5(t) = x(t) in the transient mode for small oscillations is given by formula (1), where 
(for the nominal mode) 


n == 4, m = 2, by == 0.287, by = 4.86, by == 10.2, (20) 
oz {2 . lo-?, Q> 0.205 a 0.74 Sg, 2 = 0.516 - 0.74 Sa, Wg = 2.13 -+-) 74 Si. 4 = 1 572. 


Since, in the plane of parameters oy, the system has a broad open region of stability, we determine the 
parameters sought from the condition that 4(t) approximate to the given process §t) = 1.55 (1-e*%) at five 
points. In addition to the points t, = 0 and ts = o, in which both time-domain functions coincide for any values 
of the o;,, we choose as the approximation points t, = d, ts = 2d, and | = 3d, where d = 0.25 sec. Using formula 
(14), we set up for them the system of three equations whose solution gives the following values for the gains: 


6,=7.13, 62°=1.5, 63 =6.5- 10-3, 


1556 
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Fig. 5. 


In view of the smallness of the product 035, we take og = 0. [The derivative § has order 7, which is clear 
from the expression for 5°(t).] To verify the accuracy of reproduction of the given process we find, from formula 
(1), the time-domain function 6(t) with the values of 04, 0, and og obtained as the result of our computations, 
Figure 5 shows the given curve &%t) and the actual curve 6(t). As is obvious, the maximum relative error is 
within 6%, 
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FREQUENCY CHARACTERISTICS OF RELAY SERVO SYSTEMS 


Ya.Z. Tsypkin 
(Moscow) 


A method is presented for constructing exact amplitude-frequency and phase- 
frequency characteristics of relay servo systems with arbitrary forms of external peri- 
odic stimuli. At the base of this method is the concept of the generalized character - 
istic of a relay system, a particular case of which was used in [4] for the investigation 
of periodic modes. 


Examples are given which illustrate the method described. 


In those cases when a relay servo system operates in the conditions of a periodic external stimulus, im- 
portant quantities which characterize the quality of its operation are the maximum value of its output quantity 
and the phase shift between the output quantity and the external stimulus. 


The dependence of these quantities on the frequency of the external periodic stimulus for a constant "amp- 
limde” of the latter determines the frequency properties of the relay system. 


The known approximate methods for constructing frequency characteristics, based on the harmonic balance 
method or its offshoots [1-3], are connected with a host of essential limitations. 


Among these limitations are, for example, the assumption of harmonic form of the output quantity and ex- 
ternal stimulus. Since these assumptions are not always fulfilled or justified, the approximate method can lead 
to unwarranted conclusions. 


In the present paper, we introduce the concept of the frequency characteristic of a relay servo system, and 
present a method of constructing it which is not related to any limiting assumptions. This method uses the con- 
cept of the characteristic of a relay system [4]. By the frequency characteristic of a relay system we understand 
the dependence of the maximum yalue of the output quantity of the continuous part (or a segment of it) on the 
frequency of an external periodic stimulus of arbitrary form, and the dependence of the phase shift between the 
output quantity and the external stimulus on the frequency of the latter. With this, the maximum deviation of 
the external periodic excitation is assumed to be constant. 


Preliminary Remarks 


We consider a system consisting of a relay element and a continuous portion (Fig. la), The relay element's 
characteristic has the form shown in Fig. 2a, and the continuous portion is characterized by the transfer function 


P( 
K(p) =5e. (1) 


or by the frequency characteristic 





K (jw) = U (@) + jV (@) = Ky (w) e®, 



























or, finally, by the time-domain characteristic 


n 
“10 | > Cw erv! . (3) 
vert 


h(t) ( 


In expressions (1) and (3), Q(p) is a polynomial in p of degree n and P(p) is either a polynomial or a trans- 
cendental function 


P (p,) 


P (0 
a Co = po, ” (4) 


Coo Q (0) K (C), 


and p,, are the roots of the equation Q(p) = 0, assumed to be simple, 
, stow | 1: 
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Fig. 1. Pig. 2. 









































As in [4], we denote by z(t) the steady-state quantity at the output of an open-loop relay system (Fig. 1b) 
when the input quantity varies periodically with frequency w».* As was shown in [4], the expression for Z(t) can 
be given in various forms: 


- 1 
5 {/) k. jo | s Cw a Say ; (5) 
’ @, 
if the transfer function of the continuous part is given; 


_ My Sy Ko [am — Neal 5 (2m — A) tag | O((2m -— 1) ml, (6) 
rT p | 2m — 1 


m=-1 


z (2) 
if the frequency characteristic of the continuous part is given; 


> (1) — NI ‘t-4- .. =\\ 
z (t)=-hr jh (0) +- 2 tan’ j- (ke ma be (7) 


if the time-domain characteristic of the continuous part is given. 
The conditions for the existence of forced oscillations, engendered by the external stimulus 


F(t) = Afy (ort — @) (8) 


are defined by the relay system's characteristic (Fig. 3) 








*We note that z(t), for we = const, does not depend on the form of the input quantity. 













J (w) = _?@) oe Jz (=) = ReJ (m) +- j ImJ (@) (9) 


@ 


by means of the simple constructions described in [4]. 
The critical value of the amplitude is 


A cr = |Im J () |. (10) 


It depends both on the frequency of the external stimulus and on the form of 
J (w). 


If the amplitude of the external stimulus is greater than A,,, there then exists 
a mode of forced oscillation, with frequency of the external stimulus, in the relay 
system. If the amplitude of the external stimulus is less than Ag; (A < Agr) there 
are then no forced oscillations and the system may have either a pulsed mode (if it 
is capable of autooscillation), motion with interruptions or the so-called subharmonic 
oscillations. 


J 
i 








For A > Agr, the phase shift ¢ between the external stimulus and the error 
X(t) is determined by the expression 


~ * A 
fy (x i ?) | ss lfo (—?) | (11) 


or by 


‘toad & 12 
p= fF}. (12) 


where },* denotes the function inverse (reciprocal) to }y. 


Construction of the Frequency Characteristics 


The frequency characteristic is comprised of the amplitude-frequency characteristic, defined as the de- 
pendence of the maximum of Z(t) on the frequency of the external periodic stimulus of arbitrary form (Fig. 4): 


R (@) = max | z(t) A (13) 


and the phase~frequency characteristic, defined as the phase shift 
between z(t) and f(t) (Fig. 4): 

















On (w) = arg] (t) — arg z(t) = —e + > (14) 








for A = const. To determine the amplitude-frequency and phase- 
Fig. 4. frequency characteristics, we introduce the concept of the gen- 
eralized characteristic of a relay system 





J (a, m) == — -—--+~—- jz (a *) , (15) 


where a varies from 0 to 1. Fora = 1, the generalized characteristic coincides with relay system characteristic 
(9). 


The real and imaginary parts of the generalized characteristic of a relay system 








ReJ (a, o) = — 








ImJ (x, ) = —zio 


(2 —) (16) 


bed J 


may, on the basis of formulas (4), (6), and (7), also be presented in various forms, which are given in the table. 


As a varies from 0 to 1, the imaginary part ImJ(a,w) defines the form of the periodic quantity 2 (t), 


— »9 
at) 


and the real part ReJ (a,w), in the scale of 1/w, determines the form of its derivative. It can be shown that 


J (a,w) possesses the following properties: 


ReJ (a, w) = — 


TABLE 


1 dIlmJ (a, o) 
oa 


, 


J (0, w) = —J (1,). (17) 
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for different and nonzero fields py; 


u--l an 
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for the presence of single, zero field P,, = 0. 





Note to table. In the last formulas, c'y = P(0)/Q*(0). 


To determine the amplitude-frequency and the phase-frequency characteristics, we use the relay system's 
characteristic J(w) = J(1,w). From each point w = we of this characteristic we construct, from the expression for 
the generalized characteristic J (a,w), segments of curves by varying a until these curve segments intersect the 
axis of ordinates for a = ay, and the axis of abscissas for a = ay (on Fig. 5, the values of ap, and ag are denoted 
by Qjm and ajo, respectively). Fora = am 
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Im J (a,, w) = —z (a, =| = 0, 


z(a,=)| e(e=)| 


'ReJ (a,, @)| = Fs)! == max 


w (19) 
In accordance with (14) and (18), the amplitude-frequency characteristic R(w) will equal 
R() = |ImJ (am, ¢)|. (20) 


Thus, the segments of the axis of ordinates cut off by the curves J(a,w) for w = const and a varied equal to 
the ordinates of the amplitude-frequency characteristic for w = const. Obviously, R(w) exists only for those fre- 
quencies for which A = Ac, [4]. We note that R(w) does not depend on the form of the external stimulus. 


The phase-frequency characteristic Op(w) can be determined as 
the algebraic sum of the phase shifts of z(t) and F(t) with respect to 
the error x(t) (Fig. 4), i-e., 


4, 
a Im 


J 
\ 4 How) 
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9n (wo) = $>—¢ =x (a,—1)—¢ 
or, according to (12), 


- A 
Op (0) =n (1 —a)-+ 7 (Set), (21) 


Thus, the phase-frequency characteristic is determined by the 
values of a = a» for which J (a,w) intersects the axis of abscissas; it 
depends on the amplitude and the form of the external periodic stimulus, 





If the forms of F(t) and z(t) are almost symmetric about their 
maximum values, the phase shift may then be determined approximately 
from the distance between their maxima [5] (Fig. 4). In this case, by 
setting Y,, = 9 + 1/2, we can write 


Fig. 5. 


T 4 \ 
Ope (00) Om = Rm — Gm = Nm — F—-- =K( am — sz} + (22) 


+i *(= ) 


Formulas (20) and (22) are basic for the graphicoanalytic construction of the amplitude-frequency and 
phase-frequency characteristics, i.e., the frequency characteristics of relay systems. We note that, for the simp- 
lest relay systems, the frequency characteristics can be given in analytic form. 


Example. We consider a relay servo system whose continuous part has a transfer function of the form: 


bin ky 3 
KP) = “Pp +) (Tp = 1) 





. 
From the table we determine J(a,w). In this case, the poles of K(p) are known, so it is convenient to ex- 
press J(a,w) in terms of the transfer function in the form 


2h hs FH o aT’ : P oT; 
w (7; — 72) 2 a Be. _E (23) 
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We assume that the relay system has the following parameters: kek, = 75, Ty = 0.05, and T, = 0.03. With 
these values, J(w) = J(1,w) will have the form shown on Fig. 6 [here the representation J(0,w) = —J(w) is used]. 
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Fig. 6. 


By setting w = const and letting a vary from 0 to 1, we determine the curves J (a,w) which join points with 
equal w on the characteristics J(1,w) and J(0,w). The values of a for which the J (a ,w) intersect the coordinate 
axes define ay, and ag, and the values of the ordinates equal to R(w). 


By being given, for example, the value A = 0.2, and by using the dependence of A,, on w (Fig. 7), we 
find that 


and then, from formula (20), we find @p(w). The functions Rg(w) and @(w) are given on Fig. 7 for triangular and 
rectangular forms of external stimuli. The frequency characteristics exist for those values of w for which A = A,,.. 


We now assume that T, = 0. The transfer function of the relay system will then have the simpler form: 


ky 5 


In this case, R(w) and @p(w) can be expressed analytically. Indeed, we find from (23) and (24), for T, = 0, 
that 
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By letting ReJ(a,w) vanish, we finda = ay: 


wT; - 2 


6 gt fee 28) 


By substituting this value in ImJ (a,w), we get 


= T 2 
no tb (- sig 2g) 


To determine @p we use approximate expression (22); for ay), defined by expression (28), we will have 


¢ 2 Loe 
9, () ~ 0, oT, In aes a = » }- fe | x cr | ° @9) 


1+e oT, 
Computation of the exact values of @p(w) is connected with the solution of the transcendental equation 


Im J (Qo, @) = 0) (31) 


ina». Equations (29) and (30), in a somewhat different form, were obtained in [5}. 


Figure 8 gives the amplitude-frequency and phase -frequency characteristics of the relay system under con- 
sideration for sinusoidal, triangular and rectangular forms of external stimulus Fo, constructed from formulas 
(29) and (30) for k-ky = 75 and Ty = 0.05. 


The author wishes to thank N.A. Korolev for discussing the results and for aid with the computations. 
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ANALYSIS OF TRACKING FAILURE 
IN AUTOMATIC CONTROL SYSTEMS 
IN THE PRESENCE OF FLUCTUATING NOISE 


ILA. Bol’shakoy 


(Moscow) 


The failure of tracking in automatic control systems under the action of intense 
fluctuating noise is investigated by means of the Fokker-Planck equation. The Ritz- 
Galerkin method is used to solve the boundary problem where, by a failure, we under - 
stand an increase of the error in the system above some definite quantity. For a sys- 
tem with a smoothing circuit in the form of su integrator, formulas are found for the 
probabilities of absence of failures in the system as functions of the system parameters, 
noise leyel and time. The results are used in the analysis of the noise stability of the 
AFC (automatic frequency control) system of a receiver of continuous signals. 


as 4a ass «2.2 


In the investigation of the processes in automatic control systems, one frequently bypasses the fact, due to 
the mathematical difficulties involved, that the element for isolating the error signal (the discriminator) has a 
narrow range of linearity in the majority of systems. This is true, in particular, of the servo systems of radio en- 


gineering. y 


Thus, in heterodyne automatic frequency control (AFC) systems [1], the parameter in the input signal which 
is tracked (followed) is its average frequency. Due to the narrowness of the signal spectrum as compared to the 
average frequency, the linear portion of the frequency discriminator's characteristic is also narrow, its width be- 
ing approximately equal to the width of this spectrum. In automatic range-finding devices [2,3], one follows the 
time position of a pulse whose duration is significantly less than the repetition period. With this, the linear por- 
tion of the time discriminator’s characteristic is roughly equal to the pulse duration. Similar cases occur in AFC 
systems of television scanners, in various servo systems with phase comparisons of signals, etc. 


The number of such examples can be greatly multiplied. In all these cases, the nonlinearity is unavoidable 
in principle and, for an error exceeding a certain limit, the feedback path is opened and tracking failure ensues. \ 
One of the causes of failure may be the presence of intensive fluctuating noise (for example, internal noise of the 
receiver) which, strictly speaking, cannot be assimilated to an equivalent variation of the parameter being 
tracked [3]. 





The aim of this paper is to analyze random processes in servo systems with nonlinear discriminators in the 
presence of intense fluctuating noise. It turns out to be convenient, in view of the nonlinearity of the problem, 
to use the Fokker -Planck equation for the probability density of the system errors in order to find the most interest- 
ing characteristics. 


1. Mathematical Posing of the Problem 


The block schematic of the servo systems to be analyzed is given in Fig. 1, and is roughly divided into two 
parts, namely, a discriminator and, in the other part, an amplifying-smoothing circuit together with the driven 
devices, Here, uy(t) and ug(t) are the voltages applied to the discriminator and ug(t) is its output voltage; x(t) 


el 


te 


1566 


st 





and z(t) are the input and output values, respectively, of the parameter being tracked and y = x~z is the current 
value of the error. Thus, in an AFC system we understand by x, y, and z the values of the corresponding frequen- 
cies and, in a radar range finder, the time intervals. The dual system of notation is introduced in order to em~- 
phasize the fact that, between the voltages u,(t), u,(t) and their parameters x 
and y, there is not a single-valued correspondence. Thus, u,(t) is proportional 
































oh 4 m,(t) ? to y only in the absence of noise and in a narrow range of values of y. 
(J t 
rs. u,(t) We now consider the case, of practical interest, when part 2 of the scheme 
2(t) in Fig. 1, which ordinarily contains some fundamental inertia, is an integrator 


with transfer function K,/p. The generalization of the following analysis to the 
case of a system with a higher order of astatism is possible, but leads to a sig- 
nificant increase in computational difficulty. The servo system's integrator is 
described by the differential equation 


Fig. 1. 
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Here, u,(t) is obtained from u,(t) and ug(t) by means of various linear and nonlinear operations intended 
to determine the sign and the magnitude of y. Thus, in an AFC system, these operations might include the 
parallel amplification of the signal of two tuned circuits, amplitude detection and subtraction of one result from 
the other; in a range-finder, they might include passage of the signal through two offset gates, integration within 
the limits of each of these and, again, subtraction of the results, etc. 


We shall assume that the processes in discriminator 1 (Fig. 1) occur significantly more rapidly than those 
in the closed circuit as a whole. Then, to investigate the properties of u,(t), one need not consider the transient 
response in the discriminator induced by variations in y(t), but may characterize u,(t) only by the mean yalue 
g(t) = a{y} and the spectral density in the low-frequency range 
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Sty) = | [ar@as@-+ ) — a] dr, 
which are functions of the current value of y; a{y} is called the discrimination characteristic; a typical form of 
it is shown in Fig. 2. A failure of tracking is observed if the error moves off the “slope” of this characteristic 
(for example, beyond the point —c/2, +c/2). S{y} characterizes the noise component in u,(t); in view of the 
assumptions made above, this may be considered as white noise where spectral density varies jointly with y. A 
typical form of S{y} for an AFC system is given in Fig. 3. Finding these two quantities for a discriminator of 
each concrete form is frequently a very difficult task [2]. If we consider these functions as being known, we may 
write u,(t) in the form 


ws (t) =a {y} + VS {y) ee), (2) 
where €(t) is white noise with unit spectral density. By substituting (2) in (1) we get 
d | dx > > —s 
Fc VO) = | Fe — Koa ty) |+ KV SHE. (3) 


It is known from the theory of random processes [4, 5] that, if the function sought y(t) is subject to a sto- 
chastic equation of the form* 
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*In the mathematical literature on such cases, one frequently uses the concept of a Wiener process, i.e., an in- 
tegral of white noise. 






































then the probability density W(t,y) of this function satisfies the differential equation (the Fokker-Planck equation) 
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By comparing (3) with (4) and (5), and by denoting 
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we have the following equation for the probability density of the error 
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The mathematical and physical literature on diffusion 
equations is very extensive (cf., for example, [6,4,7]). In 
the theory of oscillations, work on this question goes back 
-cf2 y to (8). This apparatus was applied to problems of servo sys- 
0 oz tems in [9], where the questions as to failures were first con- 

sidered rigorously. 
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As follows from physical considerations, and also from 
a more rigorous analysis [10], the noise component of us, (t) 
Fig, 2. can be considered as white noise only when the effective time 
constant of the dynamic system is significantly greater than 
the correlation interval of the noise voltage. When the noise 
Sty level is low, this condition is sometimes poorly met in servo 
systems, but as the noise increases, the transmission factor of 














the discriminator K, a {y}|y—o generally falls, and 
<= Z the effective time vielen is increased. In the study of the 
tL questions of tracking failures, the case of high noise levels 
7 is the more interesting one. Thus, the conditions for appli- 
Fig. 3. Output spectral density of an AFC cability of the diffusion equations hold precisely when, in 
system as a function of the error: 1) low view of the nonlinearity of the problem, the use of this ap- 
level of fluctuations; 2) high level of fluc - paratus is the most advantageous. 


— In the linear theory, the stationary random process at 


the servo system's output is ordinarily characterized by its 
mean square error. In the case of a nonlinear discrimination characteristic, of the form shown in Fig. 2, this 
process always remains nonstationary. This is explained by the probability of tracking failure which, by urning 
the system into a nonclosed integrator, leads to arbitrary wanderings of the value of y(t) with ever-increasing 
dispersion. To describe such a process, one must introduce other characteristics. Frequently interesting are the 
probability of failure from the moment of switching in of the noise up to some succeeding moment, and also the 
average (over the ensemble) of the time until failure. By the failure of the system at time t we shall understand 
that the magnitude of y(t) has exceeded some fixed level, y; = — c/2 or yp, = + c/2, with the condition that the 
magnitude was found in the interval [—c/2, + c/2] and had never left this interval until time t. This leads [5] 
to the boundary problem for the function W(t,y) in which, at the boundaries, are given the boundary conditions 


W (t,—4e)=0, W(t,4+--4c)=0 


_ (the “absorbing boundary* conditions). We choose the boundaries at the boundaries of the selected domain (pass~ 

“ band of the basic amplifier in an AFC system or the gate cutting off the receiver in a range-finder) right at the 
"slopes" of the discrimination characteristic. In contradistinction to [9], this characteristic may be of any form. 

The probability P(t) of no failure up to time t is expressed in terms of the solution of the boundary problem as 
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P(t) = ( W(t, y)dy, (8) 
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and the distribution law for the times of first failure has the form: 
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With this, the mean time until failure equals 
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The functions W(t,y), P(t), p(t) and T quite sufficiently characterize the process of servo system failure 
under the action of noise. 


2. Setting Up the Equations for the Eigenvalues 


Let the quantity to be tracked be constant (x = const) and the diffusion coefficient B be independent of y 
(B(y) = Bg). It can be shown that, with increasing noise, the dependence of B on y vanishes (Fig. 3). In section 4 
we shall discuss these assumptions. We consider the error at the initial moment to be zero. It is then required 
to solve the equation 
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with the auxiliary conditions 


W (t, —3c)=0, W(t, +) =0, WO, y) = d(y), 
where 5(y) is the Dirac delta-function. We seek a solution in the form 
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W (t, y) = >) Cie Yi (y), (12) 
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where the CT; are defined by the initial conditions. We have arrived at a problem in eigenfunctions(e.f.) and 
eigenvalues (e,v,), described by the equation 
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Since there is no exact method of solving Eq. (13) for A(y) of arbitrary form, we choose, in accordance 
with the Ritz-Galerkin [11] approximation method, a system of coordinate functions which form a complete set 
on the given segment and which satisfy the boundary conditions, In such a system, we can conveniently choose 
the functions 
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which are the ef for the equation 
BR 2 a 
> ¥"(y) = —+Y (y), (15) 


to which (13) degenerates with increasing noise, i.e., when A(y) + 0. If we discard the odd ¢, from considera - 
tions of symmetry, we get Y(y) in the form 


: n 2k — 14 2 
Y (y) = > Cx—1 co y. (16) 
k=1 


By expanding A(y) in a Fourier series in the complete system (for odd functions) %g,%4 --- 2k» ---» by 
then substituting A(y) and Y(y) in series form in (13) and then by equating coefficients of ¢g,-; (k= 1,2,..., n), 
we obtain a system of equations in the C,y-;(k = 1,2, . . .,n): 
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For (17) to have a nontrivial solution, it is necessary that its determinant equal zero: 
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We have thus obtained an equation for the n first e.v. We now choose the number of coordinate functions. 
Since the fundamental result of the analysis must be the probability of no failures up to time t, expressed by (8), 
we estimate the accuracy with which a finite number of coordinate functions can satisfy the initial conditions in 
the sense of this relationship. The function W(0,y) is expanded in a series of terms of the chosen coordinate func- 
tions in the form 
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W(0,y) =8(y)==>' cos CAF y, (20) 
k=1 
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By integrating this equation over the entire interval [—c/2, +c/2}, we get 
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By choosing two coordinate functions we replace the series in the right member by the quantity (4~ #)[1 — 
— (1/3)] * 0.85, Le., we get a 15% deficiency, while by choosing three functions we obtain a 10% surplus. From 
considerations of continuity, it is clear that the approximate values pe2)(r) and p(8)(rp, obtained with two and three 
coordinate functions, respectively, include the exact expression between themselves at the initial moment of time. 
Therefore, by carrying out the computations in parallel for n = 2 and n = 3, and by choosing as P(t) the arithmetic 
average of *)(t) and P@)t), we obtain a 2% accuracy in satisfying the initial conditions and a satisfactory accu- 
racy in the approximation of P(t) for any t, A small number of coordinate functions does not make it possible to 
determine W(t,y), even for t+ 0. However, the desired probability density is easily found by solving the linear 
problem (A(y) = K,Kuy) on the infinite line since, in the first moments, the probability of the error going beyond 
the Linear portion of the discriminator, not to mention reaching the boundary, is very small. It is easily shown 
that 
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3. Finding the Eigenvalues and Setting Up the Solution 





By taking (19) for n = 2 and by seeking the solution in the form of a series in powers of 1/8%, we obtain ap- 
: proximations for the first two e.v. 
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(the superscript indicates the number of coordinate functions used). 


For n = 3 we get 
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We remark that the terms for the partial series (22) and (23) which include the A; have one and the same 
order for different A),-;, and that the first terms increase as (2k—1)*. Consequently, in comparisca with the de- 
generate case of (15), only the first few e,v. are subject to “perturbation,” while the higher ones remain virtually 
invariant, This is explained by the fact that the higher e.v., which form a sharp peak in probability density func- 
tion (21) in the first moments of time, are damped simultaneously with the dissolving of this peak until the mo- 
ment when the inertia of the feedback, expressed by the coefficients A;, comes into play. Therefore, the varia - 
tion of the feedback parameters has very little effect on the speed of damping of the higher e.v. 


Having the e.v., we can set up the solution. Let ro sate .4 be the coefficient of the (2i—1)-st coordinate 

















function for the (2j/—1)-st e.v. (mn = 2,3, i,j = 1,...,m). Then, for n = 2, we obtain from (16) and (22), with in- 
itial conditions (20) taken into account, 
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Here we have discarded terms of the order of 1 We obtain the coefficients cf) 1,2j-1 analogously. 
Finally, taking (12), (24) and the expressions for the C3/-3 9j-; into account, we obtain, for W(t,y), 





3 3 
c® . fe Cc) ' om 
W (t, v) = r> >) 2i—1,2j—1 5 2i—1,2)—1 42»; 11 cos (2i =S)s y= 
=[(2 7 — Fear ) 0085 vt E08 * 08 ' 


= — —___—_ os — ry a i 
c 4n BR 24nB? 4n B? As 


5 (A + As) St n\2 Be nm Ay 
ser Wiheady = y | exp{— [ (3) : ey 7 | Y 


a A; — Ag oon S54. (2 — 244A) _ a A) 
4nB? e an B? 16x B2 








3r 
cos —. y — 


5 (A; + Ag) 3n\2Bi 30 As 
feo i - > y | exp{— (=) Saat] th + 


3 
os —y + yt 


A, — As TT 3 (A; — Ay) 
24nBR 16n B? ¢ 





i 5 (Ag + As) __ 5 (A, + Ay) or 
+ (=— 24n B? 16nB% ) eos e v] 4 


x exp{— (=)* 0 t} ; (25) 


where, in correspondence with section 2, we have discarded the last term in the expressions for the e.v. The 
probability of no failures up to time t is approximately equal to 
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In the course of time, the terms in (25) and (26) with higher e.v. are damped and P(t), for example, tends 
to the asymptote 


im P(e) —| 4 264i— A) Fe _ = At) ,) 
tin P«) = | 4 neem |2{—[(F) 2-3] 4. 
shown by the dashed line on Fig. 4. 


From (26) and (9) we easily obtain the distribution law for the times of failure. The average time until 
failure is approximately 














Pe 3) 

T =— + —_ ; (27) 
4B, n° BS 

where the first term equals the average time until failure with no feedback, and the second term shows the in- 

crease of this time attributable to the feedback. Relationship (27) can be verified by solving approximately the 

equation for T as a function of the initial value y(0) = yo [8] 
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for V = 0, cA(y)/B% « 1 and boundary conditions T(+c/2) = 0, and by substituting y, = 0 in the solution. 


By substituting the values of B and the A, in the form of func- 
\ tions of the circuit parameters and the input signal-to-noise ratio 
\ (S/N) in (25), (26), and (27), we obtain, for concrete servo systems, 
the functions W(t,y), P(t), and T expressed in terms of the quantities 
with which one ordinarily works in analyzing noise stability. By fix- 
ing the tracking time, we obtain the dependence of absence of fail- 
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Sas t ures on the S/N ratio. The threshold value of the S/N ratio, deter- 
mined in some way or another, obviously depends on the time inter- 
Fig. 4. Probability of no failures, as val during which faultless operation is required of the system. 


a function of time. Relationships (25)-(27) become inapplicable only for a sig- 


nificant lowering of the noise level, when failures become improb- 
able. Thus, the S/N ratio threshold is the boundary of the region in which the formulas given here are applicable. 


4, More Complicated Cases 





The method of analyzing tracking failures given in sections 2 and 3 can be extended to more complicated 
cases when the diffusion coefficient depends on the error and when the quantity being tracked varies with a con- 
stant velocity x = Vt if the basic condition, namely, a large effective passband of the discriminator, is met. 
With a constant velocity V of the measured quantity, the equation for determining the e.f. and the e.v. takes the 
form: 
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where L{ } is a linear operator. For the Ritz-Galerkin method [11], it is necessary that the left member of (29), 
after the linear combinations of the k coordinate functions are substituted in it, be orthogonal to all these func- 
tions. The presence of the velocity introduces an asymmetry in the coordinates, so that the entire set in (14) 
must be taken. As a result, for example, we have for three coordinate functions, instead of (19), the following 
equation for the e.y.: 
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By means of this equation, we can again obtain the probability of no failures. We give the expression now 
for the average time until failure, obtained from Eg. (28): 
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A comparison of formulas (30) and (27) shows that the mean time until failure decreases when a significant 
velocity V exists. This is explained by the occurrence of a velocity error which moves the “representative point" 
to one of the boundaries, For some critical value of velocity, the increase in T due to the feedback is completely 
compensated by the presence of this velocity, as a result of which any prolonged tracking becomes impossible. 


Taking the variability of the diffusion coefficient into account may be done by any of several methods. 
One of these is a transition to new functions and new coordinates [7]. We introduce the functions 
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We can then show that 
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We thus arrive at a previous equation, whose method of solution on the interval [n(—c/2),»(+c/2)] is known. 


Another method is the direct application of the Ritz-Galerkin method. The difference from the case of 
section 3 will consist of the appearance. of the harmonic function B*(y), expanded in a Fourier series on the inter- 
val [—c/2, +c/2]. 


5. Tracking Failures in AFC Systems 





An AFC system is an example of a servo system with a nonlinear discriminator for which the analysis given 
above is valid. The block schematic of an AFC is given in Fig. 5. In this case, the mixer, the band amplifier 
and the frequency discriminator properly so-called correspond to the discriminator in the circuit of Fig. 1, while, 
for the second part, there correspond the smoothing circuit, 
the de amplifier, the frequency control devices, and the 
heterodyne, u;(t) and ug(t) are the voltages applied to the 
mixer, u,(t) is the output voltage of the frequency discrim- 
inator, z(t) is the heterodyne frequency minus (plus) the 
mean frequency of the band amplifier, y(t) is the current 
difference (error) of the average frequencies of the hetero- 
dyne signal and the amplifier. The operations used to ob- 
tain ug(t) from u,(t) and us(t) include heterodyning, amp- 
lification, limiting, new (separate) amplification by the 
Fig. 5. Block schematic of an AFC: 1) mixer; two tuned circuits, amplitude detection, and subtraction of 
2) band amplifier; 3) limiter; 4,5) tuned cir- the results of detection. 
cuits; 6,7) amplitude detectors; 8) subtrac- 
tion circuit; 9) smoothing circuit; 10) fre- 
quency control device; 11) heterodyne. 











For the spectral densities of the noise and the ran- 
dom signal at the input we take, respectively, Sp = const 
and S, = exp {—(w—ws)*/A2}. With the dimensionless fre- 
quency characteristic of the band amplifier equal to 
exp {—(w—w)*/28*}, the impulsive response of the tuned circuits equal to wee * cos (wy + wy)t, with a square- 
law detector and for a, wy, As, B «K wo,ws, We determine the discrimination characteristic a {y} and the output 
spectral density S{y}. For S/S, = x > 1, we get the expressions 
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We now determine the coordinates of the “absorbing barrier” in terms of the effective width of the selected 
domain. From (31) we have that 
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Then, by using (6), (18), and (31), we can obtain the following approximate values for the Fourier coeffi- 
cients: 


A = —"s_-_ 0 
‘ape Vi + 


—_ —_ 


2V x Kk w ni? 
“aoa 











Substitution of (32) and (33) in (25), (24), or (27) leads to the sought-for probability functions. For con~- 
crete values of the parameters, all necessary characteristics of tracking failures can be obtained from these for- 
mulas, 


SUMMARY 


Random processes in nonlinear servo systems in which the element for isolating the error signal has suffi- 
ciently low inertia in comparison with the inertia of the closed-loop system, can be conveniently investigated 
by means of diffusion equations. With a significantly high noise level, the conditions for the applicability of this 
apparatus are generally met, due to the decrease in the system's gain. In this work, we obtained formulas for 
the probability of failure of tracking under the effect of fluctuating noise for any form of nonlinearity of the dis- 
crimination characteristic. 


The development of these methods, as applied to problems of automatic control, begun in [9], is apparently 
fruitful, and can be carried out in several directions. It is of interest to consider systems with high astatism. The 
solution of such problems would, in many important practical cases, deliver one from the necessity of using nu- 
merical methods for the statistical calculation of servo systems containing nonlinear elements. 

LITERATURE CITED 
{1] M.R. Kaplanoy and V.A. Levin, Automatic Frequency Control [in Russian) (Gosenergoizdat, 1956). 


[2] LN. Amiantoy and V.L Tikhonov, “The effect of fluctuations on range-finder operation,” Automation 
and Remote Control (USSR) 19, 4 (1958).* 


[3] G.P. Tartakovskii, "Radar range finders as affected by noise and by fluctuations of the reflected signal,” 
Inzhenerno-fizicheskii Zhurnal 1, 6 (1958). 


[4] J.L. Doob, Probabilistic Processes [Russian translation] (IL, , 1956). 


[5] S. Chandrasekhar, “Stochastic problems in physics and astronomy" [Russian translation of article from 
Reviews of Modern Physics 15, 1 (1943)] (IL , 1947). 


{6} A.N. Kolmogorov, "On analytical methods in probability theory," Uspekhi Matem. Nauk 5 (1938). 





*See English translation, 


[7] W. Feller, *Parabolic differential equations and the transformation sem{-groups corresponding to them” 
Coll.: Mathematics 1, 4 [Russian translation) (IL , 1957). 


[8] L. Pontryagin, A. Andronoy and A. Vitt, “On the statistical consideration of dynamic systems,” J. Exptl.- 
Theoret. Phys, (USSR) 3, 3 (1933). 


[9] A.M. Vasil’ev, “Application of the theory of Brownian motion to the investigation of noise stability of 
servo devices in radio engineering,” Nauchn. Dokl. Vysshei Shkoly. Radiotekhnika i Elektronika 2, 1 (1959). 


(10) P.L Kuznetsoy, R.L. Stratonovich and V.L Tikhonov, “Correlation functions in the theory of Brownian 
motion,” J. Exptl.-Theoret. Phys. (USSR) 26, 2 (1954). 


{11] S.G. Mikhlin, Variation Methods in Mathematical Physics [in Russian] (Gostekhizdat, 1957). 


Received February 16, 1959 

















ON THE COMPUTATION OF THE MEAN SQUARE ERROR 
OF PROCESSING A STATIONARY RANDOM SIGNAL 
BY A LINEAR AUTOMATIC CONTROL SYSTEM 


N.I. Sokolov 


(Moscow) 


A formula is presented for the approximate computation of a convolution integ- 
ral(a Duhamel integral) which has a high degree of accuracy for a large integration 
step. 


A method is given for the approximate solution (evaluation) of a double-convolu - 
tion integral which permits the rapid computation of the magnitude of the mean square 
error of processing a stationary random signal by an automatic control system. 


The relationship between the correlation functions of the input stimulus and the error of a linear automatic 
control system can be presented in the form of some double-convolution integrals [1,2]. Consequently, the prob- 
lem of computing the correlation function of the error of a linear system leads to the computation of the double- 
convolution integral obtained, An exact solution of a convolution integral is an arduous task, and is not always 
possible in practice. 


In automatic control theory, a great deal of attention has been given, and is being given, to the develop- 
ment of methods for the approximate computation of convolution integrals. In this direction, particularly great 
interest inheres in the results obtained by A.A. Krasovskii, G.S. Pospelov, Ya.Z. Tsypkin, and O.B. Lur’e. 


The method of approximate computation of convolution integrals presented in this paper is based, funda - 
mentally, on the results obtained in [3, 4]. 


The high accuracy in computing convolution integrals by the method given here, using a large integration 
step, allows one to use this method in computating the aforementioned double-convolution integrals, 


The ordinary graphicoanalytic method of computation [1,2], based on the Fourier transformation of the 
double-convolution integral, is more arduous and less accurate. 


Approximate Method for Computing a Convolution Integral (Duhamel Integral) 





If f(t) is the stimulus acting on the system and k(t) is the impulsive-response function, then the reaction 
X(t) of the system to stimulus f(t) is expressed by the integral 


t 


x (t) — V/s) k(t —+)dr, (1) 


0 


where f(r) = 0 for r = 0. 


Or, going to the Laplace transforms, we obtain 








¥ (p) = F (pyK (p), ! (2) 


where F(p) is the Laplace transform of the stimulus f(t), K(p) is the Laplace transform of the impulsive response 
k(t) or, for zero initial conditions, the system's transfer function and X(p) is the Laplace transform of the output 
function x(t). Obviously, the stimulus function f(t) can be the output of a link or a system of automatic control, 


We replace the graph of the continuous function f(t) (Fig. 1) by a 
broken-line curve which coincides with the continuous curve at points 1, 2, 
3,... [3,4]. This broken-line curve can be given in the form of isosceles 
triangles with bases equal to 2r and with heights equal to the ordinates of 
the continuous function at the points 1, 2, 3,... , shifted with respect to 
= the origin of coordinates by r, 2r, 37, ... We call pulses of this form A- 
functions. 


A(t) 





We denote the representation of an elementary pulse with unit amp- 
litude by h{p). The quantity 


(i?) 


h(p) Se TP 


tp? 


we call a unit A-function. 


According to the theorem on lags, the representation of the elementary pulse of height f (ir) at point t = 
= ir will equal f(ir)h(p)e!7P, 


Consequently, the representation for the piecewise linear approximation to function f(t) can be written in 
the following way: 


ww. 


F (p) ~h(p) >; f (iz) 


i-o 


From whence expression (2) can be rewritten in the form 


I 
X (p) =h(p) K (p) >) 1 (iz) e?* (3) 


i=o 
for 
X (p) = f(t) h(p) K (p) e-* +f (2t) hh (p) K (p)e*P 4-... 
oe tf (it) h(p) K (p)e-*? +... -+ J (ih (p) K (pe. 


By turning to the inverse Laplace transforms, we obtain 


u(t) > J (t) kag (t — t) + J (2t) kag (t-- 2t) 4-.-. 
weet f (ith kg (t —it) +... 4+ J (le) ka (t —ls), (4) 


where k,(t) is the transient response in the system when a unit A-function acts on its input and f (ir )k,(t—ir) & 
the transient response in the system when a A-function of height f (ir), shifted by ir with respect to the origin 
of coordinates, acts on its input. The values of the ordinates at the points which are multiples of r are defined 


by the expressions 


o (=) = f(z) ka (9) + f(t) ka (2t — =) + f(t) ka (Bt —1) +... 
.-+ J (2t) kg (0) -+- f (22) hag (2 — 1) 4- J (24) eg (8t — 1) +... (5) 


With a different grouping of terms in expressions (5) we obtain the following formulas: 

















x (It) = f (it) ka (72) (7) 
where f(ir) and k Agr) are numerical sequences which are multiplied by the rules governing polynomials.* 
Formulas (6) and (7) give approximate values of the integral in(1). In these formulas, the ordinates kaGir) 


are unknown. These ordinates can be found from the formulas k Ad) = L~*h(p)K(p) or 


an *P a K a 
ka (he Le — or UO) 4. [AED eee, (8) 





where r is the integration step. 


The size of the integration step is determined from the condition that the approximation of the curve for 
f(t) by the A-functions be sufficiently accurate. For practical computations, it is sometimes more convenient 
to use the formula 





Ka (yt) = ky (t + 1) (0.5; 0; — 0.5) — ky (0 + 1) (0.5; —0.5), (9) 


where kg(t) is the system's transient response (to a unit function at its input). 


If kg(0) = 0 for t = 0, then formula (9) takes the form: 





ka (yt) = ky (jt + t) (0.5; 0; te 0.5). (10) 


These formulas are the inverse Laplace transforms of the following expression [3]: 





_ o~atP ke (0 —¢e7tP) 
Ka (9) = Ko (a) = — =, (11) 


where Ky(q) is the discrete Laplace transform of the transient response of an automatic control system when a 
unit function is applied to its input and K AC is the discrete Laplace transform of the transient response in the 
system when a unit A-function is applied to its input. 


Estimate of the Error in Computing a Convolution Integral by the Approximate 
Formulas 





Formally, the approximate solution obtained by the use of the formulas given above coincides with the or- 
dinary approximate computation of a convolution integral (or Duhamel integral). In practice, the ordinary ap- 
proximate expressions for the convolution integral are almost never used, which is explained by the low accuracy 
of this method when a finite interval r is chosen. If the interval r is decreased, then the accuracy of the com- 
putation of the convolution integral increases, but there is a corresponding increase in the amount of computa- 
tional work. 


It can be shown that the computation of a convolution integral by the ordinary formulas gives a significantly 
greater error than its computation by our method. This is due to two factors. 


The effect of the first factor consists of the following. If the stimulus function is presented in the form of 
an arbitrarily chosen transient response (Fig. 2) then, to solve the convolution integral by the ordinary approxi- 
mation formulas, the ordinates are taken at points which are multiples of r. By multiplying these ordinates by 
T we obtain the hatched area on Fig. 2 which is bounded by the broken lines. Both the area bounded by this 
curve and the form of the curve are approximated very roughly. With the other method, the one used by us (Fig. 3), 
the approximation is effected by means of the triangles very accurately for the same interval r. 





*A wavy line above indicates a numerical sequence. 











The essence of the second methodical error is explained as follows. In computing the transient response 
by the ordinary approximate expression for the convolution integral, the system's transient response is taken for 
the stimulus by a pulse function, but the approximation of the stimulus is done by means of rectangles. With 
very sinall intervals r, the transient response for a pulse stimulus and the transient response for a stimulus of rec- 
tangular shape and finite duration, with a unit area, will coincide for all practical purposes. 


The greater is r, the more significant will be the difference between the transient response obtained from 
a pulse stimulus and that obtained from a rectangular-shaped stimulus and, consequently, the greater will be the 
error in computing a convolution integral by the ordinary approximate method. 
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Fig. 2. Fig. 3. 


In the method herein presented, we choose for our ordinates, not a pulse transient response of the system, 
but the system's transient response to unit A-functions at its input, as a result of which the analogous error does 
not occur. 


In computing the convolution integral by the approximate formulas given in [3], the first methodical error is 
just about equal to the first methodical error obtained with the method presented here. In contradistinction to the 
method presented here, the second methodical error appears in full force when the earlier formulas are used. 


In computing convolution integrals by the approximate formula given in [5], the second methodical error 
is eliminated, but the first remains in its entirety. 


It is clear from what has been said that the accuracy of computation of convolution integrals by the ap- 
proximation formulas given here will, for large intervals r, be much higher than the computational accuracy of 


the approximation formulas cited. This permits the approximation formulas given here to be used in practical 
computations for automatic control systems. 


Computation of Mean Square Errors 





We assume that a stationary random function rm (t) is applied to the input of a linear automatic control 
system (Fig. 4), where k(t) is the system's impulsive error response. The 
system's task is to process this random function with some mean square 


mt) | Ht) ! L(t) error. 


The mean square value of the function m(t) will equal R69). The corre- 








We denote the correlation function of the input signal by R,,(T). 


— lation function of the system's error will be defined by the expression [1] 
co co 
R(x) = \ af k(n) k (4) Rm (t + 4 — 4) de. (12) 


—o —co 
We make the change of variables: 
t+rk—y=—& yo=tr+A—F. (13) 


We then find the new limits of integration: 



























N=—oo, F=t+-+4+ 00= 00, 
N= co, —E=—t+A— co = — 0.7 (14) 


Expression (12) is transformed to the form 
R.() = | k() dr \ Rin () k (« + }—8) dk. (15) 


By changing the order of integration we get 


Re (s) = | Rm dt | kA)E(E+A—Had. (18) 
We introduce the new variables: 
tc+h—E=0, A= O—1+£, (17) 


We again find the new limits of integration: 
hk =— oo, 6 = +t—E— co = — on, 
k =o, 6 = +t—£E-+ co = oo, 


7 With the introduction of the new variable 6, expression (16) is transformed to the form 


R.(*) = \ Rim (t) dé \ k (0) k (— (« — 0 —8)] dd. (18) 
We now consider the first integral: 
ks (t«—®) = \ k (0) k{—(« —0— 8] 0. 


The function k(t — @—€) differs from the function k[-(r — @—-€)] only by the sign of the argument. 
It is well known that the impulsive response of any stable system must satisfy the condition 


co 


[1k (2) de < 00. 
0 


Moreover, from the condition of physical realizability, we get 





k (0) =0 for 0< 0, 
k{[— (« —9—§)] = 0 for 9 < t—&. 


With the conditions of physical realizability taken into account, the first integral of expression (18) assumes 
the form: 


ke(s—)= | k@kI—~@—E—O) a. 
—E 


t 
















The approximate expression for the first integral will be 


co co 
\ k (9) k{—(e«—0 —8)] d9 = (&(0)k[—(9—0)} dd= 
TE 8 
i=l (19) 
a >} k (i9;) ka (— (1%, — 8,)), 
i=j 
af-(18-10)] where 9 = r—£, 91 = 64, 6; is the time interval chosen, k,[—(16, — i6)] 


is the ordinate of the system's transient response when a unit A-function 
is applied to its input [this transient response coincides with the transient 
response k,(1 6; — 164) with the difference that the argument is a negative 
quantity (Fig. 5)], k(i@,) is the ordinate of the A-function by which the 
Fig. 5. pulse transient response k(t) is approximated. 





The second integral equals 


co 


Re (2) = | Rm (hx (te —8) db (20) 


—0o 


In formula (20), Rp)(€) and k,{7—€) are even functions (Figs. 6 and 7). We present the function k,(T —€) 
as the sum of two functions k, (¢ — &) = kx, (t —§) + kx, [— («— 8], 


where ky, (t—§&) = 0 for t—& <0, kx,[—(t-—§)] =0 for t+—E>0. 











Fig. 6. 


Then, integral (20) can be given by the expression* 


R(t) = \ Rm (6) (Ax, (t —&) + kx, [— (t —§)]} dé = 
= Rn @kx(e—H de + (Rm @) hs, [— (2 — Dd - an 


—oo 


co 
| Fn (Be (F—B4E=0, since ky (7-3) =0 for <é, 


t 


\ Ry, (5), [—(t — &)] dE =0, since =k [—(t —£)] = Ofor cy. 


—co 

















By using the approximation formulas, we can present expression (21) in the form 


R, (71) _ >) Rm (iE) ks, (75, . 13) + 7, Rm (i) ks, [— (JE: ae : i,)] ’ (22) (22) 


where R,,,(1&3) are the ordinates of the A-functions by which the correlation function is approximated. 


To determine the Ky AGEs  {€;) and the ky al-Gés ~1€))), it is necessary first to compute the ky AGEs _ 
~ 1&3) by formula (19) where, instead of k(i6,), it ts necessary to substitute k Aion, and thereafter to present 
ky (Es — 163) as the sum of two functions, as was stated above. 


The values of the ordinates k,(i6;) of the transient response 
when unit A-functions act on the input of the system are found from 


hz (t~§) 
formulas (9), (10), or from formula (8). 


If the interval €; is so chosen that the continuous function of 
— > time kg(t) can, with sufficient accuracy, be approximated by a broken 
e e line, then formulas (9), (10) will be sufficiently accurate. 


Fig. 7. 
8 , The ordinates of the functions ky AGEs — 1&3) and Ky al-Ges - 
— 1&3)] will equal, respectively, 








k... (0) * 
S— > Kea(€:), Kxa(28s), Aea(26:),---. Aeali— A), abs 


k, (0 
ate, kesl— (| heal—(— 90 eal IE os 


ksal—(—j7 +1)8&), kxal— (—7é))- 


As was stated earlier, the function k, AGEs — 1£4) is even, and the ordinates 


ka (E:) andkzal— (— E,)), ksa[— (2€,)] andkx4[— (— 26;)). 
Kkxa (3&) andtxal— (— 3€,)] 


are, respectively, equal among themselves. 
Therefore, to compute the correlation function R,(1T), one can compute one of the two integral sums given 
j 
in formula (22). For example, R, (jt) = >) Rm(i€s) kxa(/&i— i1)- The second integral sum is obtained 
i=—l 


from the computed one by replacing £; by ~€;. 

By adding the ordinates of the function thus computed, we obtain the ordinates of the error correlation func~- 
tion of the automatic control system, 

By using formula (22), we compute the value of the correlation function's ordinate for if, = 0: 


R.,(0) = >) Rm (is) Kea (Ob: — &,) = 


i=c—/ 
= Rm (— 6) kxa (16) + Rem (— Her + €1) Rea (ls — bi) + 
0 
+ Rm (— 1E&, + 28) Kea (l&y — 261) + +. + Rin (— §&) kxa(€s) + Rm (0) a, 


i=l 


Re, (0) = >} Rm (is) keal— Ob — i) = 


i=j<0 














= Rin (/€,) kzal[— (— l€,)] + Rm (7, Rix, €,) ksa(— (— le, + E,)) + 
+ Rn (lE — 2€,) ksa(— (— le, + 26,)] + eee 
k., (0 
. + Ron (G1) hxal— (—-8:)) + Rm (0) 28%. 
Since Rm (/€1) = Ry (— l€,) and kza (/E,) = kzal— (— l:)}, then 
R, (0) = Re, (0) + Be, (0) = 2m (— 1B) ka (16) + 
+ Rm (— 16, + Es) Kea (Es — &1) + 20m (— LE, + 28) ea (1G, — 261) + --- (23) 
se + 2m (— Es) Kea (61) + Rn (0) eee (0) 





i=/ 


Rz (0) = Rm (0) ka (0) + 2 >} Rin (i) kxa (i€;). (24) 


ij 


Thus, to obtain the null ordinate of the error correlation function of the automatic control system, it suf- 
fices to multiply the ordinates of the functions R,,,(ifs) and ky a(i€s) with identical subscripts for £; > 0, to double 
the product thus obtained (except for the product of the null ordinates) and to add them. In certain cases, formula 
(24) can be replaced by a simpler expression (cf. the Appendix). 


The mean square value of the error is found from the expression « =y R, (0). 


APPENDIX 


Computation of the Mean-Square Error 





Using our approximate method, we now compute the mean-square error of an automatic control system for 
a radar set with a random stationary noise function at its input. In the given case, « =/R_ (0) = Y/R,(0). The data 
which characterize the properties of this system are given in [2]. 


The system‘s open-loop transfer function has the form: 


80 (0.36p + 1) 
p (16p + 1) 





W (p) = 


The closed-loop transfer function is 


1.8p+5 
K (p)= 377 tep ts (25) 





The generalized transfer function of the system when a unit function is applied to its input is defined by 
the expression 


18p+5 


Kelp) = Tip +1 86p +5) 





and the transient response can be presented as 
ko (t) = 1 — 1.084 e—°-" cos (2.04t + 0,4). 


The curve of the transient response is given in Fig. 8. This curve can be approximated sufficiently well 
by straight-line segments chosen with an interval of € = 0.4. 


The ordinates of the transient response, taken at intervals of €; = 0.4, are given in Table 1. 
TABLE 1 





t, sec.| 0 | o4 | os 1.2 1.6 | 2.0 2.4 | 2.8 | 3.2 | 3.6 | 4.0 | 4.4 


el 0.966 
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function. 


We observe in formula 10 the ordinates of the transient-response process for reaction of a single input A 





0.5k9(jt+-7) 
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By adding the columns of this table, we obtain the values of the ordinates of k(t) (cf. Table 2), 






































TABLE 2 
t,sec | 0 0.4; 0.8) 4.2 1.6 2.0 2.4 | 2.8 3.2 | 3.6 | 4.0 
ky (t) | 0.385} 0.608) 0.270) —0.013) —0.172} —0.122) —0.016| 0.010} 0.032 0.024 0.007 


We now find the ordinates of the “nominal” system when a unit A-function acts on its input, i.e., we com- 
pute the ordinates by the formula 


I 
SY) ka (184) ky (— (1E, — 18) = hg HE) kg (HE, — ET 


ie=j 


By multiplying the numerical sequences obtained in accordance with the rules governing polynomials, we 
obtain the even functions k, AtiED- By virtue of the symmetry of these functions, it suffices to compute only the 
ordinates lying to the left (or only to the right) of the axis of ordinates. 


The multiplication process is given in Table 3. 









































TABLE 3 
iE, 0 0.4 | 08 | 4.2 | 1.6 | 20 | 24 | 28 | 82 | 86 | 40 
—(1e, —48,) |—& |-8.6 |-3.2 |-28 |-24 |-2.0 |-1.6 |-1.2 | -0.8}-04 | 0 
ka (te, 0.385 | 0.608} 0.270 |—0.013 |—0.172 |—0.122 |—0.016 | 0.010 | 0.032] 0.024] 0.007 
ka [—(1&:—48,)} | 0.007 | 0.024 | 0.082 | 0.010 |—0.016 |—0.122 ein —0.013 | 0.270} 0.608} 0.385 
0.0048) 0 ona ~0.0001 0.013}. 0010] 0.00 0.0001] 0. 0.0001 
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By summing the columns of Table 3, we obtain the values of the ordinates k, ,(i€,): 





if, —4 -3.6|-s.2| -2.8|-2.4|-2.0 |-1.s |-1.2 |-0.8|—0.4 0 











ka (i&1)| 0.003 0.014 000 0.059 0007-0. 00] -0.153| 0.19 0.05 04nd 0.637 











: 





0 f 2 7 4 ésec 


Fig. 8. 


Figure 9 shows the graphs of k, ,(i€s) (the right side), 


In order that formula (24) might be used for computing R,(0), it is necessary that the correlation function 
of the noise at the system's input be computed. 


The normalized autocorrelation function of the input signal (noise) is expressed approximately in the form 
Ry (t) =e ** !*! cos 40. 
The normalized autocorrelation function of the commutated signal will have the form: 


Ri (t) = 2 cos ,tR, (7), 


where w, is the scanning frequency, 
R,,, (+) = 2e—** '!*! cos 188 + cos 40. (26) 


The ordinates of the autocorrelation function, computed by formula (26), are given in Table 4. 


With a sufficient degree of accuracy, the function R,,,(T) can be approximated by line segments with £; = 
= 0.005, 


We note that the duration of correlation function R,,,(7) is much less than the duration of correlation func- 
tion ky a(iéy). Correlation function ky ,(i€s) is approximated sufficiently well by a A-function of width 2€, = 0.8 


1586 
































































































sec, but Rj,(T) = 0 for r > 0.4 and, consequently, correlation function R,,,(T) can be replaced by a A-function 
of width 2€,; = 0.8 sec with an area equal to the area bounded by the R,,(r) curve and the axis of abscissas. 


This area is SRin(T) = 2(1+ 1.03 - 0.43 -1.03 +... + 
+ 0.22 + 0.08 — 0.12 — 0.12 — 0,02) 0.005 = 0.0027. 


Ky (8) 


The height of the A-function of the same area and with 
width 2E = 0.8 will be equal to h = 0.00675. 


The value of the ordinate R,(0) is found from the expression 
R, (0) = hk, (0) = 0.0043. 


The actual value is Ryy(0) = 212- 107° radians*, Consequently, 





R, (0) = 0.91 « 10-*, 








The mean square value of this error will equal 


e= Vi, (0) =0.975 - 10* rad. 
The experimental value of the mean square error is 


“exp, = 1.04 - 10-8 rad. 
As is clear from these results, the accuracy in computing the mean square error is high while the difficulty 
is small (if it is considered that the correlation function of the input signal is known). 


If the duration of the correlation functions Ry (7) and ky A(iés) are commensurable, then formula (24) 
should be used to compute the mean square value R,,,(0), which does not increase the difficulty of the computa - 
tion very much. 


If the duration of correlation function ky ,(i€y) is much less than the duration of correlation function R,,(T), 
then correlation function k, ,(i€3) should be replaced by a A-function in the same way as described above. In 
this case, the correlation function at the system's output will approximately coincide with the input correlation 
function. 


SUMMARY 


1, The approximate formula for computing a convolution integral has a high degree of accuracy for a 
large integration step, and permits the transient response of a linear automatic control system to be determined 
when an arbitrary function of time, given either analytically or graphically, is impressed on its input. 


2. The computation of the mean square error of processing a stationary random signal by a linear automatic 
control system leads to the solution (evaluation) of some double-convolution integrals, 


The approximate formula given here for computing convolution integrals allows this given double convolu- 
tion integral to be solved with a high degree of accuracy with a small amount of computational work. 
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JET POWER EFFECT IN NOZZLE-FLAPPER HYDRAULIC AMPLIFIERS 


I.M. Krassov and B.G. Turbin 


(Moscow) 


The experimental results are given for an investigation of a hydraulic amplifier 
of the nozzle-flapper type for various combinations of its parameters, Basic attention 
was paid to the power interaction of the jet and the flapper. A short description is 
given of the objects tested, the program, the setup, and the methods of carrying out 
the tests. 


Hydraulic amplifiers of the nozzle-flapper type have found wide usage in the devices of automation [1]. 
In these amplifiers, the jet of working fluid flowing from the nozzle has a power effect on the flapper, which is 
one of the specific idiosyncrasies of the amplifiers’ operation. For choosing the controlling elements for nozzle~- 
flapper amplifiers and for carrying out static computations and analyses of the dynamics of systems which contain 
one or several such amplifier stages, it is necessary to have numerical data on the magnitude of the forces on the 
flappers around which the fluid jets flow. 
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Fig. 1. Types of nozzles investigated. 
The information to be found in the literature on the jet's power effect on the flapper is insufficient in 


many cases [2] and somewhat inaccurate, since it contains data for a limited range of the parameters which 
characterize the amplifier (valve diameter, working fluid pressure, gap between valve and flapper, etc.). 


In this paper we present the results of the experimental investigation of jet power effect on nozzles, ob- 
tained with increased pressures of the working fluid, positive and negative working fluid temperatures and for a 
broad range of variation of the gaps between nozzle face and flapper. 
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The three types of nozzles shown in Fig. 1 were investi- 



























gated. 


The program for testing to determine the jet power effect 
on the flapper was divided into two series. In the first of these, 
testing was carried out with constant pressure drops on a valve 
with variable cross section, i.e., on a nozzle with a flapper. Noz- 
zles of the first type (Fig. 1a) were tested, with d = 0.6, 0.8, and 
1 mm and D = 0.95, 1.15, and 1.35 mm, respectively. The pur- 
pose of this series of tests was to determine the discharge Q of 
working fluid through the variable-cross-section valve and the 
force N on the flapper as functions of the gap h between the nozzle face and the flapper for various constant pres~- 
sure drops Ap across the valve. 





Fig. 2, Membrane a and collar b of the 
laminar valve. 


The second part of the testing was carried out for sets of nozzle-flappers with valves of constant cross sec- 
tion. For the latter, we set up a laminar valve (valve block) composed of a series of membranes a clamped be- 
tween collars b (Fig, 2). The membranes were collected into blocks in which their apertures were diametrically 
opposed. The number of membranes was so chosen that there was a flow through the valve of 1425 cm*/min of 
the working fluid at a temperature of from 30 to 35°C and a pressure drop of 25 kg/cm across the valve. 


The second part of the testing program was made for various constant pressures p, at the input of laminar 
valves for nozzles of the third type, with diameters of 0.8 and 1 mm. We determined by this the discharge of 
working fluid through variable valves, the force on the flapper and the pressure p, in the intervalve chamber as 
functions of the gap h for various pressures p, and for constant atmospheric pressure pp. 
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Fig. 3. Functional schematic of the experimental setup. 


The working fluid for the tests was AMG~10 oil. In the first series of tests, the temperature of the working 
fluid was between 50 and 55°C, while three temperature intervals were chosen for the second series; —40 to —45°C, 
45-50°C, and 80-90°C, 


All the experiments were carried out on the setup whose functional schematic is shown on Fig. 3. This 
setup is analogous in many ways to that described in [2]. However, in contradistinction to the latter, the distance 
1, from the axis of rotation of arm 1 to the point of application of the load 3, which balances the force on flapper 
7, is always maintained constant. Moreover, the setup used by us was supplied with damper 5, which was a disc 
with a 100 mm diameter immersed in oil. The basic arm (shaft) 1 and the weight pan 4 were suspended on kuife- 


edge bearings. 
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Fig. 7. Results of testing a nozzle of the third type with a diameter 
of 0.8 mm and working fluid temperatures between 80 and 90°C. 


The measurement of pressure on the setup was effected by the class 2 manometers 11, the discharge of 
working fluid was measured by calibrated cylinder 12 with an accuracy of within 10 cm*/min, the gap between 
the nozzle face and the flapper was measured by micrometer 6 to within 0.005 mm and the force on the flapper 
was measured by the weight pan to within 5 grams. To prevent splattering of the working fluid, the nozzle and 
the flapper were surrounded by protective device 10. 


The tests of the first part of the program were carried out without laminar valve 9 and damper 5. 


The experimental method amounts to the following. With zero gap between the face of nozzle 8 and 
flapper 7, lever arm 1 is balanced by bob 2 in a horizontal position. Then, by means of micrometer 6, the given 
gap between nozzle and flapper is established and, by manometers 11, either the constant pressure drop, Ap = 
= Py ~ Pg, on the nozzle-flapper (first part of the program) or a constant pressure at the input of laminar valve 9 
(tests of the second part of the program) is established. Thereafter, the discharge of working fluid and the force 
on the flapper are measured. 


The force N acting on the flapper due to the jet from the nozzle engenders a torque M = NI relative to 
the axis of rotation of lever arm 1, where / is the distance from the axis of rotation of lever arm 1 to the center 
of flapper 7. This torque is balanced by the torque created by weight G on balance pan 4, which is equal to M, = 
= Gl. Since ly = 2/2, then N = G/2. For the magnitude of G we took the weight of the load in balance pan 4 for 
which the flapper broke away from the detaining device of micrometer 6. The moment of break-away was fixed 
visually by the motion of flapper 7 itself, and also by the arising of oscillations of the working fluid pressure at 
the input to the nozzle. 


The results of the experiments of the first part of the program are shown in the graphs of Figs. 4-6. The 
character of the curves shows that, for a constant pressure drop on the nozzle-flapper, as the gap h increases the 
discharge of working fluid through the nozzle-flapper first increases, and then remains constant. Simultaneously, 
with the increase of the discharge, the force on the flapper also increases, due to the increase in both the velocity 
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Fig. 8. Results of testing a nozzle of the third type with a diameter 
of 0.8 mm and with working fluid temperatures between —40 and 
—45°C. 
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Fig. 9. Results of testing a nozzle of the third type with a diameter of 1 
mm and with working fluid temperatures between 45 and 50°C. 


and the mass of the fluid coming in contact with the flapper. On Fig. 4, for example, it is clear that when the 
increase in the discharge of the working fluid through the variable valve ceases, the increase of force on the 
flapper also ceases. If the gap between the nozzle face and the flapper remains constant, then an increase of 
the pressure drop on the variable valve and an increase in the nozzle diameter give rise to an increase of force 
on the flapper. The special features of the designs of the nozzles investigated did not show any noticeable in- 
fluence on the magnitude of the discharge or of the force. 


The results of the tests in the second part of the program for nozzles of the third type and with laminar 
valves, for various temperatures of the working fluid, are shown in the graphs of Figs. 7-9. It follows, from a 
consideration of these figures, that with a constant pressure at the input of the laminar valves there is, with an 
increase in the gap h, an increase in the discharge of working fluid but a drop in the force on the flapper and in 












the intervalve chamber pressure. As the discharge of working fluid becomes more stabilized, the pressure in the 
intervalve chamber and the force on the flapper (in the Limits posed by the location of gap h) approximate to 
constant values. With a constant gap, an increase of pressure pp, leads to increases in the discharge, the force and 
the pressure in the intervalve chamber, but an increase in the nozzle diameter has a more telling effect on the 
pressure drop in the intervalve chamber. 


A reduction of the working fluid's tempcrature leads (for the same values of supply pressure p», nozzle di- 
ameter d and gap h) to a decrease in the discharge Q, the force N and to an increase in the pressure p, in the 
intervalve chamber (Figs. 7 and 8). 


A comparison of the results obtained with those given in [2] shows that the functions N = f (h) for pp = const 
and also N = f (ps) for h = const have an identical character. However, the functions N = f (h) for p; = const are 
essentially different. This, in our opinion, is due to the fact that the data of the experiments given in [2] com- 
prise a range of variation of the force N (for constant pressure drop across a variable valve) within the limits of 
which the variations of the gap h do not lead to an increase in the discharge of working fluid through the nozzle- 
flapper. 

The results obtained from our experimental investigation allow one to design hydraulic amplifiers of the 
nozzle-flapper type with increased working fluid pressures, larger nozzle diameters and larger gaps between noz- 
zle and flapper, and so are useful for practical work. 


SUMMARY 


1, An experimentai investigation was made of nozzle-flapper type valves both with and without constant 
cross-section valves, for working fluid pressures up to 90 kg/cm, nozzle diameters up to 1 mm and gaps between 
nozzle and flapper of up to 0.5 mm. 


2. New data were obtained characterizing the jet farce effect as a function of a number of system param- 
eters, 
LITERATURE CITED 


{1] Fundamentals of Automatic Control [in Russian] edited by V.V. Solodoynikov. Vol. 2, part 1 (Mashgiz, 
1959). 


(2) V.N. Dmitriev and A.G. Shashkov, “Jet action on flappers in pneumatic and hydraulic controlling or- 
gans of the ‘nozzle-flapper* type,” Automation and Remote Control (USSR) 17, 6 (1956).* 


Received April 21, 1959 


*See English translation. 











st 


e- 


59 











INVERSION AS A WAY TO RATIONALIZE 
CLASS H RELAY CIRCUITS 


I.L. Oifa 


(Moscow) 


Gra phic inversion of plane circuits is carried out by the method of transfigura - 
tion, and the equality of this way with the inversion of spatial circuits is proved. The 
method cited opens up a formally reasoned way to transform plane circuits into spatial 
ones. 


Posing of the Problem 





The graphical method of transforming class H circuits to their inverses is the most widely used. The known 
method of inversion amounts to this, that each closed contour of the original circuit is replaced by a node. Simul- 
taneously with this, the elements in the lines of the circuit are replaced by their inverses, and series connections 
of elements are replaced by parallel ones, and vice versa. Every class H circuit structurally consists of nodes and 
contours connected continuously one to another. Since contours cannot be formed without nodes, it is clear that 
an inverse circuit can be obtained by carrying out one of some characteristic operations: either all contours are 
replaced by nodes, or all nodes are replaced by contours, as is done in transfiguration. 


In fact, if the nodes of the original circuit are expanded to contours in the new circuit, then the contours 
of the original circuit must be transformed to nodes since, first of all, in the new circuit as in the original one, 
nothing but nodes and contours is to be found and, secondly, in the transformation of the original circuit, only the 
nodes are expanded. In other words, replacement of the contours of the original circuit by nodes and replacement 
of nodes by contours are equivalent transformations, which lead to the inversion of plane circuits. 


Inversion by Transfiguration 





The standard graphical way of inverting consists of determining, from the original circuit, the locus of the 
nodes of the inverse circuit, so that the successive replacement of each loop (contour) by a node is carried out on 
the same drawing. This method makes it possible to orient graphically the dual circuit. In practice, it is more 
convenient to carry out the inversion by the successive replacement of each node by a contour, since the graphic 
construction in this case is carried out on a new drawing. 


As an illustration, we carry out the inversion of the circuit, borrowed from [1] (Fig. 1), by individual steps 
from node to node in the original circuit M. 


We first expand node 1. For the circuit of this node we determine that, in the first closed loop of the in- 
verse circuit, there must be the following elements: a, c and #. We construct this loop in the form of a triangle 
(Fig. 2). From the circuit of node 2 we determine that, in the second loop, the elements f, e + d, and k must 
be contained. 


We construct the second loop on the first so that element f is common to both loops. We carry out the ex- 
pansion of the other nodes of the circuit in the same order. Naturally, the number of sides in a closed loop must 
equal the number of rays of the node in the original circuit. 





The expansion of each node should be carried out in one, previously 
; established, direction (clockwise or counterclockwise) and the expansion 
| 5 should be begun each time with the element in common to the two neigh- 
e ; boring nodes. If these requirements are met, the construction of the in- 
1 . e ai verse circuit is carried out automatically. These requirements are conse- 
4 | quences of the topological-geometric theory of cizcuits by which relay 


? ead circuits are considered as subsets of the set of connected graphs [2]. 





The beginning and end of the circuit are determined in the follow - 
2 ing way. 


r ay The shortest paths from X to Y in the initial circuit were afk and 
i 





b(g+i). Consequently, the beginning of the inverse circuit must be the 
3 point at which all the elements of the first short circuit are connected in 
y parallel (a + f + k), and the end must be the point at which all the ele- 
ments of the second circuit are connected in parallel (b + gi). Finally, 
Fig. 1. the concepts of “beginning” and “end” of a circuit are purely conventional. 
These are the points at which the circuit is supplied. 


We now place the circuit obtained over the original one (Fig. 3). Having convinced ourselves that the 
number of elements in both circuits is the same and that all elements are replaced by their inverses, we compare 
the structural admittances of the two circuits (Gx-yi Gy-w). If Gy-y = 1, then Gy.y = 0 and conversely. Con- 
sequently, the circuits are inverses, q.e.d. 


Scheme N 

















Fig. 2. Fig. 3. 


Special Features of Graphic Inversion 





Obtaining inverse circuits in the class of the simplest pi-section circuits presents no difficulty, and is car- 
ried out either analytically or graphically with the same ease. For example, we represent graphically the circuit 
described by the equation F = n(ab+cd) (Fig. 4a). 


We now invert this circuit graphically by replacing the sole node in it by a loop, while simultaneously re- 
placing all the elements making up the branches of the original node by their inverses (Fig. 4b). 


We then write the equation satisfied by the inverse circuit thus obtained: 


F = n+ (a+ d)(c +d). 
If we write the inverse equation 


F = n(ab +- cd), 





we can convince ourselves that the expressions are isomorphic. 








al. 


re 


bey 
cuit 


re~ 














Two conclusions follow from this: 
1) The inversion of the original circuit by graphic means was correct. 


2) The search for the supply points of the inverse circuit obtained should be carried out by comparing both 
circuits and, if the current path in the normal circuit was anb then, in the inverse circuit, the supply point will 
be the point at which the inverses of these same elements are connected in parallel (a +n+b), f.e., all the in- 
verses of these elements are connected in parallel. 
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Fig. 4. Fig. 5. 








We now carry out the graphic inversion of the simplest H-circuit (Fig. 5). The analytic expression which 
describes the original circuit (Fig. 5a) has the form: 


F = a(b -+- kd) + c(d + bk). 


This form of the equation is not isomorphic to the circuit structurally, but is isomorphic to it dynamically, 
in the sense of the processes occurring in it, since the increase in the number of terms in the equation over the 
number of circuit elements is fictitious, 


We obtain the inverse of this expression analytically: 


F = [a + b(k + d)\fe + d(b+ k)| = ac +-ad (6 + k) + cb (k 4- d) + 
+ bd (k-+- d)(b+ k) = ac +- abd-- adk +- chk + chd + bd = 
= ac + adk +- chk +- bd = c(a + bk) +- d(b + ak). 


If we write the analytical expression for the graphically obtained inverse circuit, starting with elements c 
and d, we can easily convince ourselves that the inversion was carried out properly. 


This example allows us to draw the following conclusion. 


When the nodes of an original circuit are expanded into loops of the inverse circuit, the ends of the branches 
should be considered either as the following nodes of the original circuit or as supply points of the circuit, de- 
pending on where the branch under consideration ends. Thus, in the first branch of the first node one element a 
(and not the two elements ac) is contained; in the second branch of this node there is again just one element k 
contained, etc. 


We now carry out the inversion of a somewhat more complicated circuit (Fig. 6a), using the graphical 
method, 


By placing the circuit obtained (Fig. 6b) on the original circuit in such wise that the nodes of the resulting 
circuit are inside the loops of the original circuit, we convince ourselves that the inversion was properly carried 
out. This example allows still another important conclusion to be drawn. 


If to the poles of the circuit there are connected, not two, but three or more branches, then the poles of 
the new circuit should be sought in the boundary circuits of the original circuit which contain the least numbers 
of elements, i.e., in the shortest circuits. 
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Domain of Applicability of the Graphical Method 





The graphical method of inversion is applicable to all known types of relay circuits, including those which 
are called “nonplanar” in the theory, or spatial circuits, An example of a nonplanar circuit is given in Fig. 7. 
A circuit of this type is called nonplanar because it contains arms which do not 
intersect in the plane of the drawing. On Fig. 7, such arms are those which 








T | contain elements b and k and the one containing g. 
1 me However, it is known from relay circuit theory that any H-circuit can be 
t- } expressed as the sum of several simple pi-circuits. This expansion is carried out 
e ° by the method of determining the circuits which act on the given elements [3]. 

d 
f—e— + Thus, as applied to relay circuits, there are sufficient grounds for asserting 
x 1 that, in the set of these circuits, there cannot exist a circuit the dynamics of which 





would be impossible to represent without intersections in the plane of the drawing. 
Fig. 7. In other words, one can assert, in correspondence with the method of determining 
the circuits which act on a given element, that circuits constructed for the same 
set of conditions and differing from each other by the presence in one of them of 
some greater number of single-valued elements, do not differ in theory, but only structurally, i.e., H-circuits are 
simplified sums of series-parallel circuits. 


Inversion of Nonplanar Circuits 





We now invert the nonplanar circuit of Fig. 7. A trivial way to invert such a circuit is to first decompose 
it and then to invert it by parts [4]. 


We decompose the given circuit into two circuits by determining the circuits on which element a and ele- 
ment g act (Fig. 8). Both circuits, the sum of which is equivalent to the original circuit, are laid on a plane 
without mutual intersections, and may be easily inverted, as is done in Fig. 9. 
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The inverted terms of the original circuit are multiplied, i.e., are joined in series. We then carry out a 
graphic simplification of the circuit, by which we obtain an inverse circuit with the same number of elements as 
in the original circuit. Since this cannot always be done we then, to convince ourselves that the circuit obtained 
is actually the inverse of the original one, carry out the operation of inversion on the inverted circuit and, since 
f = F, the result of this operation must be the original circuit. 
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Fig. 9. 


The direct result of this repeated inversion is the circuit of Fig. 10. This circuit is equivalent to the original 
one, since the closed loops ebked are completely identical and, in any case, their mutual connections can be 
considered as one contour. The points where the initial and final elements are connected to such a loop are also 
known. Certainly, the inversion is correctly executed. 
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The circuit of Fig. 9, however, is inverse to the original circuit only by the existence of the transformation 
carried out on it. It is clear from what has been done that not every spatial circuit can have a dual, in the full 
sense of this definition. 


SUMMARY 


The present work provides a foundation for the following conclusions. 


First, although the inversion of a nonplanar circuit by the application of our method did not lead to thé 
obtaining of a dual, it is clear from the proof of inverseness given above that, by applying inversion as a way to 
synthesize sequential circuits, one can succeed in simplifying them until a spatial circuit is obtained. 


Second, in spite of the fact that the method of transforming bridge circuits to series-parallel circuits by 
isolating the circuits which act on given elements makes it possible to describe any complex scheme by homo- 
morphic analytical expressions without the use of special symbols (3 and ‘|| ) and, seemingly, opens a path for 
applying, without Limitations, the algebra of relay circuits to such expressions, the graphical method, in fact, 
turns out to be the most effective in practice. One can, with its help, obtain the desired results much more rap- 
idly and with less expenditure of effort. 


For example, the analytic expression which describes the circuit of Fig. 7 is 


F =a{(e + bked) m + (bk + ecd) n} + g{(c +dkbe) m + (d + cebk) nj. 


The inversion of this expression is 


= {a + [e(b +k +d +c) + myo + ke +c +d) +n} x 
x (et+icd+k+5+e) + miyde+e+b +h +n)}. 





Simplifying this expression analytically is significantly more difficult than doing it graphically, as was 
done on Fig. 9. 


Third, graphic inversion can be implemented either by compressing contours into nodes or by the trans- 
figuration method. 
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MAGNETIC LOGICAL ELEMENTS FOR AUTOMATIC CONTROL CIRCUITS 


N.P. Vasil'eva and N.L. Prokhoroy 


(Moscow) 


The functional schematics of basic systems of logical elements based on mag- 
netic cores and semiconducting diodes are considered. 


Comparative evaluations of these systems are made. 


INTRODUCTION 


In many systems for the automatic control of productive processes, relays frequently play the part of the 
elements lying between the sensor elements (transducers) and the executive elements. 


The basic purpose of the relays in such schemes is to execute various elementary operations which permit 
the requisite character of executive element control to be obtained as a function of the form of the signals ap- 
plied to the relay inputs. The logical operations in this case as, to be sure, in the majority of discrete (digital) 
devices, are based on the binary code, for which only two possible output states are used, For a relay, these are 
the two different positions of the contacts. 


Magnetic logical elements are advantageously employed to increase reliability, simplify usage and increase 
the length of service of such systems. Such elements have no moving contacts and are also characterized by two 
output states, namely, the presence or absence of an output voltage. 


The basic simplest logical functions, on whose basis all possible control operations can, with sufficient 
flexibility, be implemented are the functions of “and,” “or,” “not,” “memory,” and time delay. We now briefly 
characterize these functions. 


The "and" element has several input circuits, and a signal appears at its output only in case it has signals 
simultaneously at all its inputs. By its action, this element is analogous to a circuit consisting of several series- 
connected relay contacts. Such a circuit is closed only when all the contacts are closed. 


The “or” element also has several input circuits, but a signal appears at its output whenever there is a sig- 
nal applied to at least one of its inputs. This circuit is thus analogous to a circuit consisting of several relay 
contacts connected in parallel. 


The “not” element has one input circuit, and a signal appears at its output only when there is no signal at 
its input and, conversely, where there is a signal at its input, its output signal disappears. This element is anal- 
ogous to a relay with normally closed (break) contacts, The "memory"* element has two inputs and is charac- 
terized by the fact that, when a signal is applied to the first input, a signal appears at the output and remains 
there even after the signal at the input is discontinued. In order for the output signal to disappear, it is necessary 
that a signal, perhaps only a momentary one, be applied to the second input. Thus, a “memory” element is 
analogous, to a polarized relay, or to a lock-up relay. 


The time delay element can implement various functions, analogous to those characteristic of relays with 


*This element can also be called a “dynamic memory” or a “dynamic trigger.” 








time delays. For example, a “switch-in delay” can be implemented, when the signal at the output appears with 
a definite lag after the appearance of an input signal, but disappears simultaneously with the removal of the in- 
put signal, or a “switch-out delay” can be implemented whereby the signal appears without any lag (i.e., at the 
output) but disappears with a definite lag. 


The basic logical functions considered above can be implemented by the most diverse circuits and ele- 
ments, With only diode circuits all the logical functions can be obtained, but the signals in such circuits will be 
damped, and intermediate amplifiers have to be connected in the circuits. The use of magnetic cores in com- 
bination with diode circuits allows one to obtain the logical functions together with the amplification necessary 
for operation of circuits with large numbers of individual elements. Here we shall consider only magnetic logical 
elements consisting of rectifiers and of cores with rectangular hysteresis loops. 


Comparative Survey of Circuits of Magnetic Logical Elements 





Circuits of magnetic logical elements, based on the use of cores with rectangular hysteresis loops, can be 
classified in various ways, depending on what is adopted as the basis for the classification. Thus, for example, 
as a function of the form of the output voltage, logical elements can be divided into pulse elements and elements 
with de outputs. Another essential basis of classification by which different formis of logical element circuits can 
be grouped is the manner of connecting the load with respect to the elements’ working windings. All circuits 
may be divided, on the basis of the load connection, into circuits with elements in series with loads (or circuits 
with voltage sources) and circuits with elements in parallel with loads (circuits with current sources), 


To the class of circuits with elements in series with loads are the circuits of pulse elements which underlie 
the Ramey [1] fast-acting magnetic amplifier and the ordinary one- and two-core magnetic amplifiers in the 
amplifier and relay (switching) modes, 


To the class of circuits with elements in parallel with loads are the pulse one~-core transformer circuits, 
which are widely used in domestic computing and remote control devices [2],.and circuits of pulse one-core ele- 
ments with parallel connection of stages. 


A, Pulse Element Circuits of Elements in Series with Loads 





There are two types of logical elements whose operation is based on the Ramey circuit. 


Figure 1 shows the functional schematic of the first of these, developed by the Westinghouse Corp.* in the 
United States [3]. With no signal at the input, the operation of this circuit is analogous to the operation of the 
Ramey circuit with shorted input. This is attained by having a flow of current in the circuits of diodes 1 and 2, 
given an auxiliary supply, the magnitude of these currents being somewhat greater than the magnetizing currents 
of windings w, and wz, respectively. Thus, diodes 1 and 2 are open for the magnetization polarity-reversing cur- 
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Fig. 1. Functional schematic of a repeater. Fig. 2. Circuit for an “or” element. 


rent of both windings. The switching in of the signal U,, which is of sufficiently great magnitude, cuts off diode 
1 and the diode circuit of the control winding, and output voltage U, appears at the circuit's output, since the 
core will still not have its magnetic polarity reversed. Thus, the basic cell of elements of this type is a repeater 
with a half-period shift. Such a cell can form the basis of shift registers and pulse counters. 


The "or" function is obtained very simply with elements of this type by connecting the signals through an 
auxiliary circuit which contains a number of diodes equal to the number of inputs (Fig. 2). 


*This same principle of operation underlies the elements of the UNIVAC AF/CRC computer [4]. 
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circuit shown in Fig. 3. In theory, the circuit permits 
the reliable logical multiplication of any number of 
inputs, but, with this, the input circuits must contain 
three times as many diodes as there are inputs. 


j j j The “and” function is obtained by means of the 
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The “not” element, or inverter, is obtained by 
means of a simplification of the circuit of Fig. 1. 
The circuit for such an element is shown in Fig. 4, 
UL, VL2 where the control voltage U, is connected in the con- 
Fig. 3. Circuit for an “and” element. trol circuit instead of the supply voltage. For the crea - 
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Fig. 4. Circuit for a “not” element. Fig. 5. Circuit for a memory 
element. 


tion of a memory element in this system of logical elements, at least two “not” elements and two two-input “or” 
elements are necessary. The memory element has the form of a trigger (flip-flop), whose block schematic is 
given in Fig. 5. 


The advantages of the logical elements just described are, first, the independence of reliable system opera - 
tion on fluctuations of the voltage supply within quite wide limits (430%), which is attained by virtue of the bias 
of the diode circuits and, second, the capability of obtaining the basic logical functions “and,” “or,” and “not* 
with a theoretically unlimited number of inputs on one and the same repeater cell without controlling or varying 
the parameters. This latter capability makes systems of these elements enormously flexible, and decreases the 
total number of elements necessary for implementing any given control problem. 
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Fig. 7. Circuit for an “and” element. Fig. 8. “Circuit for a memory element. 














The disadvantage of this type of logical element circuitry is the large number of diodes, particularly in 


the logical elements which implement the “and” function. 
: 


The second type of logical element circuitry which is also based on the Ramey circuit consists of the cir- 
cuits developed by the French company SEA [5]. The functional schematic is shown in Fig. 6. The special fea - 
ture of this type of circuit is the common voltage supply for both parts of the circuit, i.e., the controlling and 
the operating parts. Due to this peculiarity it is impossible, without the additional winding w';, to obtain the 
logical function of signal repetition when the elements are connected in series, since one and the same source 
voltage is the magnetic polarity-reversing voltage for preceding and succeeding elements. Thus, without the ad- 
ditional winding, each element implements the “not” function, i.e., is an inverter. This element can also serve 
as the basis for shift registers and pulse counters, but with inversion from element to element. To obtain the 
function of repetition, one must connect winding w, via input 1 to the voltage source (the dashed line on Fig. 6) 
and, as the input winding, use winding w';. Then the ampere-turns of the signal will compensate the demagnet- 
izing ampere-turns of winding w, and, when a signal is applied to the input, the core's magnetic polarity will 
not be reversed, i.e,, a voltage appears at the output and, when the signal is taken off, the output voltage be- 
comes equal to zero. 


Thus, depending on the connection of windings w, and w';, the element u'r consideration can implement 
the function either of repetition or of “not.” 


Joining this element with an “or” circuit (Fig. 2) allows one to obtain the direct or the inverse (depending 
on the connections) “or” function for a theoretically unlimited number of inputs, just as with elements of the 
first type. The "and" function cannot be obtained directly with just one such element. 


Figure 7 gives a circuit which allows the inverse “and” function to be obtained for two inputs, Ug; and Uy. 
Thus, to obtain the direct "and" function of just two signals, it is necessary to connect one "not" element to the 
output of such an element. Consequently, in circuits requiring a large number of multi-input "and" elements, 
the use of logical elements of this type is less efficient than the use of elements of the first type, since a much 
larger number of elements is required here. For constructing a memory element in this system of logical ele- 
ments, at least two standard elements are also necessary (Fig. 8). 


Logical elements of this type permit a voltage source fluctuation of 410%. The advantage of this type of 
element as compared with the first type is the decrease in the number of diodes. However, the logical elements 
of the second type also have disadvantages, namely: they are very sensitive to source voltage fluctuations and 
have a less flexible structure, i.e., a large number of elements is required when the use of “and” elements with 
several inputs simplifies the circuit as a whole. 


B. Pulse Element Circuits of Elements in Parallel with Loads 





The basic difference between the operation of pulse elements with elements in parallel with loads from 
that of elements in series with loads is that a load voltage appears in elements of the parallel type only when the 
core's magnetic polarity is reversed, whereas a load voltage appears in elements of the series type only when the 
core is saturated and cannot reverse its magnetic polarity. 
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Fig. 9. Circuit for an “or” element. Fig. 10, Circuit for a “not” element. 


In the group of circuits of logical elements with loads connected in parallel, one should first mention the 
transformer element [6], which is widely used in computing technology and is being introduced into remote con- 
trol technology. Figure 9 shows the functional schematic of such an element, which is implementing the "or" 
function with three inputs. Single-polarity pulses from a current source pass through the core via special winding 
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Wp- If a signal appears on at least one of the inputs Ucgy—Ucs in an interval between these pulses and reverses the 
core’s magnetic polarity, then the succeeding pulse from the source winding reestablishes the previous saturated 
state of the core, wherein load voltage Uy, appears on load winding wy at the moment of magnetic polarity re- 
versal by the supply pulse. With the core's magnetic polarity reversed by the input signal, the voltage in the out- 
put winding is blocked by the diode. Using elements of this type, one can construct shift registers, pulse counters 
and other circuits. 
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Fig. 11. Circuit for an “and” element. Fig. 12. Circuit for a memory element. 


To obtain the “not” function in an element similar to the one just described, special winding ws (Fig. 10) 
is introduced, the ampere-turns of this winding being directed toward the ampere-turns of winding wy, blocking 
its action when a signal is applied to the input of winding ws. Thus, signals Uy and Ugg must be introduced in 
order to obtain the “not” function in this element. 


The element which implements the “and” operation in this system of logical elements has two cores (Fig. 
11). The element of the circuit of Fig. 11 has a nonzero output voltage only when the first core only has its 
magnetic polarity reversed, This occurs only when 
there are signals on both inputs, Uy and Ugg, which in- 
duce, in the second core, ampere-turns which are equal 
and oppositely directed. 


For the creation of an “and” element which per- 
mits the logical multiplication of a large number of 
signals, the circuit is further complicated, Because of 
this, circuits with logical elements of the transformer 
type are so constructed that the use of “and” elements 
is avoided, or so that the number of these elements be minimal. 








The memory element in this system of elements also requires the use of two cores, as do all the other cir- 
cuits of pulse logical elements. Its circuit is shown in Fig. 12. 


The basic advantage of logical elements of the transformer type is that they operate reliably with fluctua - 
tions of the supply voltage within wide limits (450%). This is explained by the fact that the integral ob- 


tained at the output of each element, is determined only by the magnitude of 2w,%s of the core and, starting 
with some necessary minimum quantity of magnetic polarity-reversing ampere-turns, does not depend on the 
magnitude of these ampere-turns. 


The disadvantage of the logical elements of the transformer type is the low flexibility of the circuits when 
the “and” element is not used or the increase in the number of elements when elements implementing this func- 
tion are used. 


We turn briefly to the elements where the load is connected directly in parallel with the working winding. 
The circuit of two such elements, one of which is connected to the output of the other, is shown in Fig. 13, This 
circuit can form the basis for shift registers or pulse counters. It possesses the following special features. 





The circuit of one such element is a repeater, in which an 
output voltage appears when its core has its magnetic polarity re- 
versed in a working half-period, since the load is connected in 
parallel with the working winding. When another such element is 
connected in as the load, the working half-period of the controlling 
element is the controlling half-period of the controlled element. 
In this half-period, the reversal of magnetic polarity of the cores 
of the first and second elements occurs simultaneously and, if the 

ty ts 
elements are identical, then =| u,dt = ~ \ u,dt. In this circuit, uy and u, are one and the same voltage 
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Fig. 14. Circuit for an “and” element. 


and, seemingly, w, must be equal to wz. However, this would lead to the circumstance that a succeeding ele- 
ment would have the same effect on a preceding element as the preceding element has on the succeeding one. 
Thus, for wy = Wg, there will not be unidirectionality of action in the elements. Unidirectionality is easily ob- 
tained in these elements, just as in elements implemented by other circuits, by making the number of turns of 
the input winding significantly smaller than the number of turns of the output winding, so that the back effect is 
effectively negligible. In order to dissipate the excess voltage on the input element, ballast resistor Rj, is con- 
nected in series with the input (Fig. 13, shown in dotted lines). 
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Fig. 15. Circuit for a shift register without inversion. 
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Fig. 16. Circuit for a shift register with inversion. 


It must be mentioned that the circuit of Fig. 13, which is being considered here as a circuit with the load 
connected in parallel (a circuit with a current source), can also be considered as a circuit with the load connected 
in series (circuit with a voltage source) if, as the input and output quantities, we take the voltages on the im- 
pedances® and not on the choke windings. Depending on what one considers as the input and output quantities 
(voltages on the impedances or on the chokes), the functions “and” and “or,” implemented by this circuit, are 
interchanged. Thus, if this circuit is considered as a parallel one, then the circuit with the diodes connected to 
produce the “and” function for any number of inputs will have the same form as the circuit with the diodes con- 
nected so as to give the “or” function with a series-connected load. Figure 14 shows such a circuit with three 
inputs. With a parallel connection of the load, the “or” function is obtained in a more complicated way, similar 
to the circuit which implements the “and” function using elements with series-connected loads (Fig. 3). 


In carrying out our purpose of making a comparative survey of the pulse magnetic logical element circuits, 
we formulate the following conclusions. 


1, In all the circuits considered for creating memory elements, at least two simple elements, forming the 
arms of a flip-flop, are required, and none of the circuits has any essential superiority over the others. 


2. The two-winding logical elements with biased diode circuits (the circuits of Figs. 1-5) have single-core 
modifications of the basic logical functions “and,” “or,” and “not” for any number of inputs, and operate reliably 
within a broad range of supply voltage fluctuation (430%), but have a larger number of diodes than the other cir- 
cuits, 


*It is easily shown that the circuit of Fig. 1 can be taken to this same form. 
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3. The three-winding logical elements (the circuits 


arn of Figs. 6-8) have fewer diodes, but permit fluctuations 
Y out of only +10% in the supplied voltage, and do not have 
Vin | Vout single-core modifications of the logical function “and” 
- for an arbitrary number of inputs. 
a b 


4, The transformer logical elements (the circuits 
of Figs. 9-12) operate reliably over a wide range of sup- 
ply voltage fluctuations (+50%) and have fewer diodes 
(in comparison with the circuits of Figs. 1-5), but have 
no single-core modifications of the logical “and" function. 


Fig. 17. Pulse supply: a)The circuit; b) The 
form of the output pulses, 


5. Logical elements with loads in parallel permit the easy implementation of the “and” and “or” functions 
with any number of inputs and with fewer diodes than the elements of Figs. 1-5, but have | wer gains and permit 
smaller fluctuations of the supply voltage (4+ 10%). 


6. All the types of logical elements considered are roughly equivalent from the point of view of power de- 
livered for equal physical dimensions. 


1. The repeater elements of each type can be used as the basis for shift registers (cf. for example, Figs. 
13, 15 and 16). To convert them to ring counters, it is necessary to supply counting pulses, rather than supply 
pulses, to the circuit elements and, at their inputs, to supply starting pulses. Registers constructed of elements 
of the first, third, and fourth types shift pulses without inverting them, while elements of the second type invert 
the shifted pulses. 


8. The dimensions of all the types of logical elements considered can be significantly lessened if short 
pulses of the same frequency are supplied to them (Fig. 17b). It is necessary to mention, at the same time, that 
the pulse supply (Fig. 17a) also executes simultaneously, the function of a stabilizer, since it gives a constant 


t 
magnitude of \ uy dt. This is particularly important for circuits of the second type, which are relatively sen- 


0 


sitive to fluctuations in the supply voltage. 


On the basis of all that has been said, one may conclude that, of all the systems of logical elements con- 
sidered, the most promising for industrial application is that whose elements have biased diode circuits, since 
the system of these elements is quite flexible, i.e., allows one to create any modification of “and,” “or,” and 
“not” for any number of inputs and with a minimum number of cores, and permits large fluctuations in voltage 
supply. The disadvantage of this system is the large number of diodes which, when silicon diodes are used for 
objects working in difficult temperature conditions, raises the cost of these elements. 


In industrial objects, where the dimensions of the elements are of significance, it makes sense to use sys- 
tems of small-size elements, drawing from a pulse source. 


C, Circuits of Logical Elements with DC Outputs 





Single-core circuits in which the feedback coefficient is brought to the necessary value by means of a con- 
denser are used as logical elements by the General Electric Company [3]. Figure 18 shows the circuit of such an 
element. Thanks to the presence of the condenser, the average value of the current in the feedback winding is 
much higher than in the working winding. 


By varying the number of turns in the feedback winding and the capacity of the condenser, one can change 
the feedback gain of this circuit and obtain different slopes of the element's input-output characteristic. With a 
sufficiently high feedback gain, the element enters the relay (switching) mode, and can be used as a memory 
element. This same circuit can be easily made the base of the logical “not” element (the characteristic of Fig. 
19) and an “or® element. The “and” element, in this system of logical elements also, is more complicated than 
the other elements. For reliable operation of the “and" element, it must have as many cores as there are inputs 
(Fig. 20), so that the full output voltage only appears when all the cores have their magnetic polarity reversed. 


A disadvantage common to a?" these circuits is the necessity of using condensers, which either significantly 
increase the size of the elements or are the cause of variations in circuit parameters, which can lead to the 
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Fig. 18. a) Functional schematic of a repeater; b) re- 
peater characteristics. 
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Fig. 19. Characteristic 
of the “not” element. 





























Fig. 21. a) Circuit of a memory element; b) characteris- 
tic of the memory element. 
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Fig. 22. Circuit for a universal element. 


necessity of controlling and adjusting the circuits during the period of usage. The “and" element, moreover, is 
relatively complex to use, since it contains several cores with individual and common windings. In our opinion, 
1 the memory element is the only one clearly superior to the pulse elements considered above, since only in this 
{ system of logical elements can it be implemented in a one-core modification. 


The memory element ¢ircuit and characteristic in Fig. 21a and b) operates reliably for significant fluctua- 
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Fig. 23. Characteristics of the universal element. 


tions in supply voltage, with the bias for the windings coming from a common source. It has two input windings, 
one for switching on, the other for switching off. 


In our opinion, it is sometimes advantageous to use this modification of the memory element in conjunc- 
tion with pulse “and,” “or” and “not” elements. 


D. Logical Elements Based on Ordinary Magnetic Amplifiers 





The General Electric Company [3], in addition to single-core elements which execute simple logical func- 
tions, also developed logical elements which implement more complex functions. These elements are ordinary 
magnetic amplifiers with several control windings and combinations of input impedances, One of the control 
windings is used for biasing. 


The circuit of a universal element, connected for the implementation of several logical functions, is shown 
on Fig, 22, This element has three inputs A, B, and C, and one output. Inputs A and B are applied in parallel to 
one of the two control windings. Input C is supplied to the second control winding, which can be connected in 
opposition relative to the first. 


The two impedances R, in the circuit of signal C allow, by switching of these resistors, the magnimde of 
the controlling ampere-turns of signal C to be varied by a factor of two for one and the same signa! level. 


The circuit also contains a step-by-step control of the magnitude of the bias, by means of which the ain- 
pere-turns of the bias can be varied from zero to three times the unit signal ampere -turns, 


Figure 23 gives the element's input-output characteristics for some of the logical functions which can be 
implemented by this element. The “or” function for three signals, A + B + C, is implemented by applying sig- 
nals to one of the C inputs and to the A and B inputs. The ampere-turns of the bias must equal the ampere-turns 
of any of the th cee signals. 


To obtain the logical "and" function, A- B- C, the magnitude of the bias ampere-turns must be triple that 
of each of the signals. The function of inversion (negation), “not,” is obtained by applying a signal to input A or 
input B and with no bias. 








With no bias, and with identical magnitudes of signals A, B, and C, the element implements the function 
AB + AC + BC, 





If the first and second inputs of signal C are connected in parallel, one can implement the function of trans- 
mitting C unless it is inhibited by either A or B, i.e., (A + B)C. 


The “memory” function is implemented in this element by introducing in it, in the relay mode, an addi- 
tional positive feedback path through resistor R,. The signal for switching off the memory is applied to inputs 3 
and 4, connected in parallel, which allows the memory to be switched off even when a signal to switch on is 
present, 


The logical elements obtained are quite universal, and can replace a number of monofunctiona! elements 
but, just because of this, they are more cumbersome and complicated than the monofunctional elements. More- 
over, they execute the necessary functions reliably only with stable signal magnitudes. Therefore, if the plug-in 
principle underlies the construction of logical circuits, the use of multifunctional elements is inefficient. 
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INFLUENCE OF EDDY CURRENTS ON THE CHARACTERISTICS 





OF MAGNETIC AMPLIFIERS WITH FEEDBACK 


AND LOW CONTROL CIRCUIT IMPEDANCE 


M.A. Rozenblat 


(Moscow) 


The theoretical influence of eddy currents in cores with rectangular hysteresis 
loops is considered as it bears on the characteristics of the basic types of single-cycle 
magnetic amplifiers with feedback, including fast-acting amplifiers, for low im- 
pedance in the control circuit for various forms of voltage supply curves. The basic 
conclusions are borne out by experiment. 


Magnetic alloys with rectangular hysteresis loops (50NP, 65NP, and others) are widely used in contemporary 
magnetic amplifiers. This allows a significant decrease in the weight and dimensions of an amplifier for a given 
power, a lowering of its inertia and an improvement in the linearity of its input-output characteristic. The 
characteristics of such amplifiers, particularly when high-quality germanium and silicon diodes are used, depend 




















Fig. 1. 


essentially on the form of the portions of the dynamic hysteresis loop by which the 
magnetic state of the amplifier’s cores varies [1]. Deviations of the form of the 
parts of the dynamic hysteresis loop from the form of the corresponding static 
hysteresis loop can be due to eddy currents, magnetic viscosity and inhomogeneity 
of the magnetic properties of the core material along its cross section. 


The present work is devoted to the theoretical consideration of the effect 
of one of these factors, namely, eddy currents, on the characteristics of the basic 
types of magnetic amplifiers with feedback, in the case when cores with ideal 
rectangular static hysteresis loops, such as shown on Fig. 1, are used. It is assumed 
that the ratio of the outer to the inner diameter of the cores is sufficiently small 
so that one can ignore the inhomogeneities of their magnetization due to the de- 
crease in magnetic field strength as the point of the core under consideration is 
removed further from the core's geometric center [2]. 


The mode of operation of a magnetic amplifier and the effect of eddy cur- 
rents depend essentially on the impedance of the control circuit. In this work, we 


consider the basic two-half-period magnetic amplifiers (Fig. 2) and also fast-acting magnetic amplifiers follow - 
ing the circuit of Fig. 3, all with very small impedances R, in the control circuit: 


we 
R<Z 4, 


where w_ and w,, are the number of turns in the working winding and the control winding, respectively, and R 

is the active impedance of the working circuit, equal to the sum of the load impedance R,, the forward im- 
pedance of the diodes Ry and the working winding impedance Rg. With this assumption, R, cannot limit the mag- 
nitude of the current at double frequency directed into the control circuit. We consider load impedance R; to be 


purely active, forward impedance Rg of the diodes to be constant and the back impedances to be infinitely large. 
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Fig. 3. 


As the basis of our accounting for the effect of eddy currents, we take the equations for the dynamic hys- 
teresis loop which were obtained by the author in (3) for various modes of magnetic polarity reversal of cores 
with rectangular static hysteresis loops. From Eqs. (2) and (3) of [3] we can obtain, in general form, the follow- 


ing expression for the dynamic hysteresis loop: 












lys~ 








B \ dB, n- 
H=+H,+ Smyd? (12+ p-) F 10 10. (1) 


where H,, is the static coercive force, By is the saturation induction, y is the specific electrical conductivity, and 
d is the material's thickness. 


The plus signs in formula (1) correspond to increasing induction, i.e., to the ascending arm of the hysteresis 
loop, and the minus signs correspond to decreasing induction, or to the descending arm of the hysteresis loop. 


Double-Half-Wave Amplifiers with Sinusoidal Voltage Supply 





With the assumptions made above, all the circuits of Fig. 2 have identical characteristics [1,4] and so we 
may limit our consideration to the first of them (Fig. 2a). 


We give the subscripts 1 and 2 to the parameters appertaining to cores 1 and 2, respectively. We now con- 
sider the half-period when the supply voltage e = E,, sinwt through diode 1 is applied to winding w_ of the first 
core (Fig. 2a). Under the action of this voltage, induc- 
tion By varies from some initial value By for wt = 0 to 
the saturation induction B, for wt = a, after which load 
current i increases by a jump from the low value of mag- 
netizing current i, to a value determined by the active 
impedance R of the load circuit. Simultaneously with the 
variation in By there is induced, in winding we of the first 
core, an emf e, which is applied to winding w, of the 
second core and forces its induction B, to vary from the 
initial value B, to value Bg at the end of the half-period. 
In the following half-period, the two cores interchange 
roles. In the steady state, By = Bag. 











D, The character of change of load current i and the 

+ magnetic state of the cores during one half-period, with 
the effect of eddy currents taken into account, is shown 
in Fig. 4. 


Field strengths H, and H, are induced, respectively, 
MH, by the component i, of total current i,, in winding w. of 
the first core and current i, in winding w, of the second 
core: 
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Fig. 4. 
where / is the mean core length. 


When inductions By and B, vary simultaneously, the following magnetizing current must flow through wind- 
ing w. of the first core: 


Vc. 


, . l 
ty = ly v_ is >= Vinw_ (H,—A,). (3) 


At the beginning of the process of induction change, H; = H, and H, =—H,. Therefore, the initial value 
of magnetizing current is 


2H, 
. oe ‘ 3’ 
‘no = Qanw (3") 


~ 


With no eddy currents present, this value of i, is maintained during the entire process of change of induc- 








tions By and By. If the voltage e = i,9% then current i in winding w. does not reach a value sufficient to change 
induction By, and the load circuit manifests itself as a purely active circuit, in which the current is 


i on E,,, sin wt 2 (4) 
R 
This occurs for 0 = wt =a», where 
i ok 
O@) = arcsin E— , (5) 
™m 


For ag = wt Sa, the load current i = i,,, and is defined by formula (3). After saturation of the core for 
wt = a, the load current is again defined by foemuta (4). The mean value of the load current is 


i 7 - (\ 4 dust +- indwt + ( z dos] ; 


aj= 


h, —— 
By integrating, we obtain 


E 
I= <a (2 — C08 a + cosa) + /,,, (6) 


where the mean value of the magnetizing current is 


oe =\ (As — H,) dwt. 


0 





(7) 


We note that the magninide of ag is ordinarily so small that we can, in practice, take ag = 0. In fact, the 
magnitude of i,, 9 in formula (5) is less than the minimum load current, |, min = Iy(a=n) By expressing E,, in 
terms of the maximum value of load current, I, .,3,: 


Ea = ¥ Ju eR, 
we find, from (5), that 
“ _ 2) ; 2 
@, < arc sin —Lmin_. — arcsin ——— 
* Lnax ®K ny 


where Kypy is the multiplicity of load current variation. 


Even for Ker = 10, we have ag < 0.064 and cosa, > 0.998. Therefore, in order to simplify the calcula- 
tions, we shall henceforth take a, = 0. Then, instead of (6), we obtain 


En 
Ly = op (1 + cosa) + Jy. (8) 


We now find the expression for the average current I, in the control circuit. In the intervala swt =1, 
the induction B,; = B, = const, and variation of B, can occur only under the action of a signal voltage E,. How- 
ever, with the assumptions we have made (R, -+ 0), the magnitude of E, is very small, and cannot give rise to 
any essential variation of B, in the time interval cited. We may, therefore, consider that By, just as By, varies 
only for 0 swt sa. Fora swt <1, current ig, in view of the very slow change of induction B,, cannot in- 
crease in absolute value beyond the value 

Hi 
4 = D4nwe * (9) 
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9) 








Thus, in taking into account that the change in the core's magnetic state occurs along the descending arm 
of the hysteresis loop, we obtain for the control current; 


{ F a bad 
fe ‘| Oma, \Hatot —f aoe | (10) 
0 a 


Equations (8) and (10) are the amplifier’s input-output characteristic, 


If we ignore the voltage drop Ri,, in the load circuit caused by the magnetizing current, then we can de- 
scribe the variations in inductions By and B, for 0 = wt <a by the equations 





t 
B, = By + ay \ Emsin wtdt = By + kB} (1 —cos wt), (11) 

sg 

108 ¢ 
B, = B,—-—; \ En sin widt = B, — kB, (4 —coswt), (12) 

where 
AB,, E,,,-108 

k= 2B, su ow, SB, — 


is the ratio of the maximum possible induction change AB,,, during the half-period (for a = 1) to twice the saturs- 
tion induction. In the sequel, we shall assume that k = 1. For wt = a, induction B, = By and, therefore, we have 
from (11) that 


By = B, — kB, (1 — cosa). (14) 


By substituting the values B = B, and B = B, in (1), we find the corresponding expressions for Hy and Hg, after 
which we can integrate Eqs. (7) and (10): 


a 
5 os Star = \ (2H + n®/yd*kB,- 107 (2 + k + kos a — 2k cos wt) x 
0 


x sin wt] dot = 57——| 2H. = + 2n/yd*B,k-10™* (1 — cosa)], (15) 


l xn i ¢ 
I.= iy | a | [He + nf qdh*B,-107 (1 — 008 wt) sin wt] x > 





l = 
* dat = View [H. + On/yd*k*B, - 10 lo (1 — cos a)*}. (16) 


If we neglect the component I, of the load current, then we can obtain, from (8) and (16), the following 
equation for the amplifier’s input-output characteristic in relative units: 


] 2nfyd*k*B I 
EE Re te ao 
fp =-1- 75 ( = ) (17) 
u L max. 


We note that the dynamic coercive farce of a core for a sinusoidal induction with crest value B,, = B, 
equals* [3] 


Hog= He + *fyd*B,-10°. (18) 
*This expression can be obtained from Eq. (1) by setting B = 0. 
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Equation (17) may, therefore, be given in the following form: 


\2 


Ie 2 H od I 19 
mo —teF() (tr), qi 


L.max 





Double-Half-Wave Amplifiers with Pulse Voltage Supply 





In many cases, magnetic amplifiers are supplied by virtually rectangular pulses (Fig. 5a) for the purpose 
of reducing the size of the amplifier for a given power. We denote by E the peak value of the supplied voltage 
pulses and by r the duration of the pulses, and we then consider 
e the operation of the circuit of Fig. 2a under the same assumptions 
previously made. 


During an operating half-period, the induction B, of the first 
core is varied, under the influence of supply voltage e = E, in cor- 
Ir . respondence with the following law (Fig. 5b): 





















| 7 
' i 
' i ; t 
’ a; 
, . B,=Byo+ ay | Eat = By + 
aij re 3 é 
3 
b + 50 t =Byy+2kB,* (20) 
art, re | Tt 
4 5, where By is the initial value of By for t = 0 and 
“8, ee we ww ww we me eee ow ewe = = = 
is Ex-10° (21) 


Pig. 5. =" 2w_SB, . 
Here, just as in the case of a sinusoidal voltage supply, k is 


the ratio of the maximum possible variation during one half-period of the induction of one of the cores under 
the influence of the supply voltage to twice the saturation induction. 


For induction B, we have 


B, = B,— 2kB, ~. (22) 


For t = t,, induction By becomes equal to B, and the first core saturates, while the load current reaches the 
value of 


i = E/R. (23) 


With this, B, reaches the value Byy = By and, due to the small value of the voltage signal, does not vary 
further until time t = T/2. Starting with t = T/2, the cores exchange roles. 


The average value of current in the active load is found from the formula 
t 


For the values of H; and H, which define the magnitudes of the magnetizing current I, (7), we find, from 
(1), (20), and (22), that 


2es—' 
ae ah 4) 


nyt 
a's 


aes i,dt +- 


0 


7 
a 


By 


H,= H.+ —|— 


qT —— 
pe mr or (1+ 5 


+ 2k- ‘), (25) 
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2nyd*k*B .t 





Hy= — He—-—aag nai 
By taking into account that 
Byy = By = By — 2kB, =, (27) 
we find that 
ta 
i= a 7 \ (H,—H,)dt= ties [He + a0" |, @8) 


0 
For a constant component of control current we find, by analogy with (10), that 
be T/2 
l 2 l 2nyd*k*B, tt 
le = Vanwy | \ Het = \ Hedt| = came (— He — ae a) - @9) 
0 a 


If we neglect component I, of the load current then, from (24), (29), and (9), we can obtain the amplifiers 
input-output characteristic equation 


I 2nryd*k?B ] S 
ha ~ ek Fi Trai " r—*-) , (30) 
u c . L max 
where 
2E 
= am 1 
Ty ynax Rk T° (31) 


By taking into account that T = 1/f, we find from (17) and (30) that the input-output characteristics of 
amplifiers with the circuits of Fig. 2 coincide (to within the magnetizing current 1, in the load circuit) when the 
amplifiers are supplied with sinusoidal voltages or with voltage pulses. 


Double Half-Wave Fast-Acting Magnetic Amplifiers 





We now consider fast-acting magnetic amplifiers, based on the Ramey circuit (Fig. 3), for the case when 
the signal voltage e, is a double-half-wave rectified (unsmoothed) sinusoidal voltage with crest value E,, either 
coinciding in phase or 180° out of phase with the supply voltage (Fig. 6). At the same time as the line voltage 
is being applied to working winding wi of one of the cores, via the corresponding diode, the second core is "de- 
magnetized” by the action of the total voltage 


é = E., sin wt — Ey sin of, (32) 


which is applied to its w,, winding. Here, E},, = E,jWc/w.. The change in induction in this core during one 
controlling half-period is defined by the formula 


t 
108 Ec 
B, = B, ——e\ edt = B,—kB (1-3) (1 — cos wt) (33) 
2 8 Ss 2 8 8 ’ 
x 


™m 


where formula (13) holds for k, At the end of the controlling half-period (wt = 7) we get 


, E 
Bix Buy B, — 2D; (1-3). (34) 


\ m 
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As before, the load current is defined by formula (8) with 


a 
IN the sole difference that the magnetizing current now equals 
wW 
tle 
“vA uN le = Sea al H,, dot, (35) 
0 

















, i.e., the supply of the working circuit in the circuit of Fig. 3, un- 
0 5 wt like the circuits of Fig. 2, attains, in the interval 0 = wt =a, the 
. x sx jx magnetizing current I, necessary for changing the induction of 
ae"e only one of the cores from By to Bs. Therefore, the magnitude of 
ten | I, in this circuit is less than for the circuits of Fig. 2. 
0 wt 
fi ax JX Since Eq. (14) is valid for the amplifier now under considera - 
Fig. 6, tion then, by taking (34) into account, we find that 


E 
cosa = 2—S — 1 
E 


“m 


and the load current (8) can be expressed in the following way: 


E I 
—< 4 7 ——- (36) 
Em L.max 


I) 
I 
L max 








From (1) and (33), we now find the expression for the field strength H, induced by current i, through the 
signal source: 


E \3 
H, = — H, — n*f/yd*k*B,-10-° (1 -- #) (1 — cos wf) sin wi. 


m 


By taking into consideration that, in the amplifier we are discussing, the induction B, varies during the en- 
tire controlling half-period, we find the following formula for the constant component of the control current: 


ae © l . E.\? 
vee Was | Hatt = ase | — He — BefnaP¥*By 80 *(1— +) | (37) 


If we ignore the second component of the load current in formula (36), which is due to the magnetizing 
current for one of the cores, we obtain, from (36) and (37), exactly the same Eq. (17) for the input-output charac- 
teristic of the fast-acting magnetic amplifier of Fig. 3 as we obtained earlier for the ordinary double-half-wave 
amplifiers of Fig. 2. 


Double-Half-Wave Amplifiers with an Arbitrary Form of Supply Voltage Curve 


We now show that the coincidence observed above of the input-output characteristics of amplifiers for dif- 
ferent laws of variation of induction B, during the controlling half-periods is not of a random character, but is 
also valid for any form of supply voltage curve, e = e(t). 


Assuming that, for t = 0 and t = T/2, voltage e = 0, and that e(T+t) = e(t), we denote, as before, 





T/2 
\ e(t)dt =k2wSB,-10% = (k <4). (38) 
9 


If the first core in the circuit of Fig. 2a is saturated for t = t, then, for the average value of load current, 
neglecting the magnetizing current which flows for 0 =t <t,, we obtain the expression 
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T/2 
\ e (t) dt. (39) 


ta 


i 
I= R 


4] i) 


For t, = 0, the load current attains its maximum value: 


4kw_ SB, 
I max = ~~ RT .408 ’ (40) 


For 0 st St,, inductions By and By vary in accordance with the following laws: 


t 


t 
ea wa Det ale (t) dt (41) 
0 


9 





B,= 





10 


For t = 0, we have By(0) = By and By (0) = B, and, for t = t,, By(ty) = Bs and By(t,) = Beg = By. Since 


2 ta 


T/2 
| (d= | ee) de— \ e(e) at, 


te 


a 


we then obtain, from (38)-(41), that 


2w_S-10-8 
Ty =1 1 wax —jr — (8, — By). (42) 


We now find the magnitude of the control current. By setting B = By in Eq. (1) and integrating, we find that 


RB? 
\ H,dt = — Ht + 5ryd?- 107% (x, rig ‘ 





2B, 
Therefore 
a T/2 
ly = ge re [| Het — \ Hat) = oie — H. FR 1B. — By). (43) 


By eliminating By — By from Eqs. (42) and (43), and by taking into account that T = 1/f, we obtain exactly 


the same Eas, (17) and (30) for the circuits of Fig. 2 as were obtained when the amplifier was supplied with sinus- 
oidal and pulte voltages, 


Experimental Data 





For the experitnental verification of the basic theoretical results obtained above, it was necessary to have 
cores with rectangular static hysteresis loops for which the relative influence of magnetic viscosity and inhomo- 
geneity of properties along the cross section would be negligibly small in comparison with the influence of eddy 
currents. For this purpose, vores made of alloy 65NP, 0.35 mm thick, were used for the experimental amplifier 
After annealing in a magnetic field, the cores had a coefficient of rectangularity (i.e., the ratio of the residual 
induction B, to the induction B,.{n a field B = 5 oe) exceeding 0.99 for a coercive force of Hg = 0.025 oe. The 
cores and pes amplifier had the following basic parameters: d = 0.35 mm, J = 78 mm, S = 31,5 mm, y = 4-104 
ohm cm™}, B, = 13000 gauss, f = 50.cps, w, = W., = 1200 turns, R = 150 ohm, and k = 0.9. Specially prepared 
plane germanium diodes were used in the amplifier in order to eliminate the effect of feedback currents on the 
amplifier's characteristics, The function l/h, max = f1(i,.) was computed by formula (17) without taking mag- 
netizing current I,, into account. This fuxction is shown by curve 1 on Fig, 7, Curve 2 in this figure represents 
the function 1, /l, max = f(I,), computed trom formulas (8), (15), and (17). By summing the ordinates of curves 
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i 1 and 2, we obtain the computed input-output characteris - 
> aay tic of the amplifier with magnetizing current being taken 
L max into account (curve 3 of Fig. 7). Curves 4 and 5 of Fig. 7 






















































































ad are experimental characteristics, obtained for the circuit 
of Fig. 2b when the amplifier was supplied by a sinusoidal 
, ag voltage (curve 4) and by rectangular pulses (curve 5). The 
[| 1 experimental characteristics for the other circuits of Fig. 2 
a a 2 | 08 virtually coincided with the characteristics obtained for 
/ the circuit of Fig. 2b. 
4 ) 
= Se ey a7 
Mi | 
“Tt taeet | 
See B r | 
seni L ets i 4 as 
4 ‘. : 
fh : 7 
—_—_———_———— —— = = sé Fi 
a 7 ; Af. H 
/ ! ' 
// | HU 
+—03 Hh 
4 | HH 
, il 
me td = 02 ' 
4 H ' 
/ Fig. 8. 
ee L. ] 
7 I pew 
JS TP oe Z 3 The curves given attest to a satisfactory coincidence 
TAs of theoretical and experimental characteristics, particularly 
wa - “08 -0b Ol 0 


in the region of small load currents. The discrepancies 
Fig. 7. observed in the region of large load currents between the 
theoretical and experimental characteristics, and also be- 
tween the experimental characteristics corresponding to 


sinusoidal and to pulse-supply voltages have a completely regular character, and are explained by the effect of 
magnetic viscosity. 


Figure 8 shows the oscillogram of the dynamic hysteresis loop of one of the cores for a sinusoidal variation 
of the induction with a frequency of 50 cps.. The dashed lines represent the static hysteresis loop, obtained by the 
ballistic method. The light lines, denoted by the letter T, correspond to the theoretical dynamic hysteresis loop 
computed from formula (1). It is obvious that, when the sign of the field strength changes, there is a noticeable 
lag in the variation of the induction in comparison with theoretical curve T, this lag not being caused by eddy 
currents, since the induction is practically invariant with this. Essential differences of the experimental dynamic 
hysteresis loop from the theoretical one are only observed at the initial portion of the descending (or ascending) 
arm of the hysteresis loop, which also explains the character of the deviations of the theoretical curves of Fig. 7 
from the experimental ones. 


The relative effect of viscosity, as compared to that of eddy currents, will be the greater, the slower the 
variation of magnetic flux. When an amplifier is supplied by a sinusoidal voltage, the speed of flux variation 
in the controlling half-period, with large load currents, is significantly less than for an amplifier supplied by 
rectangular voltage pulses. Therefore, the effect of viscosity is greater in the first case, and the upper portion 
of the amplifiers input-output characteristic (4) is shifted to the left in the region of large absolute values of 
the control current (Fig. 7). 


In the fast-acting amplifier of Fig. 3, the variation of induction in the controlling half-period occurs even 
more slowly than in the ordinary circuits of Fig. 2. This is explained by the fact that the induction varies from 
B, to Bag during the entire controlling half-period in the fast-acting circuit, whereas this occurs in the circuits of 
Fig. 2 only during the interval 0 = wt =a. We should, therefore, expect a more significant viscosity effect in 
the fast-acting circuit, which is attested to by the input-output characteristic of this circuit (curve 6 of Fig. 7). 
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In the region of small load currents, when the duration and speed of magnetic flux variation in the controlling 


half-period are practically identical for the circuits of Figs. 2 and 3, the characteristics of these circuits also co- 
incide., 


It is interesting to note that the presence of viscosity can lead to an improvement in the rectangularity of 
the dynamic hysteresis loop, an improvement in the linearity of the input-output characteristic and an increase 
in its steepness, As the frequency of the supply voltage is increased, the relative effect of the eddy currents in- 
creases, and a closer coincidence of theoretical and experimental characteristics is observed for all the circuits. 


SUMMARY 


In this work we considered theoretically the effect of eddy currents on the characteristics oi the basic types 
of magnetic amplifiers with feedback for very small impedances of the control circuit and with rectangular hys- 
teresis loops. Analytical expressions were obtained for the input-output characteristics of magnetic amplifiers 
with eddy currents present. It was established that their quantitative effect is the same for all types of amplifiers 
considered (to within the component of load current due to the magnetizing current of the unsaturated core), in- 
dependently of the form of the supply voltage curve. 


The experimental characteristics, obtained for a test magnetic amplifier, were found to be in sufficiently 
close agreement with the theoretical ones. With this, for 0.35 mm thick alloy 65NP, even for a supply frequency 


of 50 cps, there is a very significant effect of eddy currents on the characteristics of magnetic amplifiers with 
feedback. 


A significant effect of magnetic viscosity on magnetic amplifier input-output characteristics was observed, 
even for a marked development of a surface effect. The effect of viscosity, as distinguished from that of eddy 
currents, depends on the amplifier circuit and on the form of the supply voltage. It is the stronger, the more 
slowly the variation of magnetic flux occurs in the controlling half-period. 
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ON THE QUESTION OF DETERMINING THE OPTIMAL 
COEFFICIENTS OF IMPULSIVE RESPONSES 


A.M. Perel*man 


(Moscow) 


A formula is presented for determining the impulsive response of a system which 
will be optimal for “white noise.” 


In the case of an input signal consisting of a controlling action 


r 

: tt 

g(t) =>) 8) = (1) 
i=0 


and a disturbance in the form of “white noise," whose spectral density we denote by Sg, the system's optical im- 
pulsive response is defined by the expression [1, 2] 


K (t) = = nit. (2) 
i=0 


The coefficients y; of expression (2) are found by solving the system of r+ 1 equations 


4 pit-i+m 
+ no — (m=0,4,...,7), (3) 
® j—n 


i+itm = 
where the yyy are forces which we will consider as given. 


The solution of system (3) of equations can be expressed in terms of determinants 


r 
A 
“n= > Pm (i =0,41,...,r), (4) 


m=0 


where D is the determinant of the system and the A;)j are the co-factors of the corresponding elements. 


The first subscript of a co-factor corresponds to the number of the row in the determinant (the numeration 
starts with zero) and the second subscript corresponds to the column. 


The terms py Amji/D for up = (—1)™m! which enter into expression (4) are the coefficients % of the op- 
timal m~differentiator filter which is compensated for the first r—m derivatives. 


The solution of system (4) of equations is quite arduous for r > 3. We propose the following formula for 
determining the ratios A,,,;/D: 
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Am i 4yitm AN pm Arie So (5) 
p ae (ml)? (+m + 4) (il)? pms 


where a term of the form AQ denotes the number of ways of placing p things taken q at a time, i.e., the variation 
(permutations) of p things taken n at a time. 


Expression (5) allows the following formula to be obtained for determining the coefficients y;: 





r 2m+1 2i 
1= 50 >) (— Ay y,, Artima Aries (i= 0,4,..00) (6) 
han (ml)? (i + m + 1) (il)? TE 


Formula (5) was obtained inductively by considering the solution of system (3) of equations for a number 
of values of r and for m = 0,1,2, Table 1 gives the numerical values of the coefficients y; (Sy and T are omitted) 
for the optimal filters (ug = 1; u4=...= py = 0). Table 2 gives the values of y; for the optimal differentiator 
filters (uy =—1: pg = pg =... = py = 0). In both tables, between neighboring values of y, and also yj and 7j,3, 
there are given the ratios of the following coefficient to the preceding one. A consideration of these ratios al- 
lowed a number of regularities to be remarked, for which the numerators and the denominators of the ratios are 
written in the form of sums of regularly varying terms. The result of generalizing these regularities is formula (5). 


We consider two examples of the determination of the coefficients y; by formula (6). 


Example 1, To determine the coefficients of the impulsive response, optimal for “white noise,” of a double 
differentiator filter compensated by the first derivative. 


We have that r = 3, m = 2, jg = 2°, tg = a = fg = 0. Formula (6) takes the form: 
Ati+1 


4+i So 
+ 3) (@)® Ps 





6 
n= (— 1)? op 


From this we easily obtain the coefficients y» = 480 Sp/T®, 74 = —5400 Sp/T*, 7g = 12960 Sp/T®, and ys = 
= —8400 Sp/T*, which satisfy the given forces. 


Example 2. To determine the coefficients of the impulsive response, optimal for “white noise,” of the un- 
compensated triple differentiator filter. 


Here, r= 3, m = 3, pg = 32, wo = wa = wg = 0. Formula (6) takes the form: 


7 2i+1 
= —(— 1)'+3 Re. | Ay+i So 
a= 3. (i+ 4) (ly? pHs * 





From this we easily obtain the coefficients 7» = 840 Sp/T*, v1 = —10080 Sp/T®, yg = 25200 Sp/T®, ys = 
= —16800 S,/T", which satisfy the given forces. 
SUMMARY 


The formula presented allows one to determine easily the coefficients y; of the impulsive response, optimal 
for “white noise,” which satisfies given forces. This same formula can also be used for solving systems of linear 
algebraic equations whose numerical coefficients, reading from right to left and from top to bottom, comprise 
the number series ¥;, ¥%, Ys, ¥4, etc. 
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ON SETTING UP AN OPERATING PROGRAM 
FOR SEVERAL CONTROL SYSTEMS 


E.G. Vol’pert 


(Moscow) 


Frequently, a dispatcher or operator has to observe a significant number of quantities which characterize 
a process being controlled. For such observations there have recently been developed the so-called "run-around 
devices" which permit successive observation of the various parameters of the controlled process. However, in 
the observation of various quantities during equal intervals of time, there is the danger that some rapidly varying 
quantity will go beyond its permissible limits before “its turn” comes around. Conversely, a slowly varying quan- 
tity will be observed too frequently. 


We present below a method for improving the quality of operation of such systems by setting up a program 
based on probability relationships. 


Let n quantities be observed. We shall assume that, from the probabilistic point of view, these quantities 
are independent, or are weakly dependent. In fact, if there existed a functional dependence X,, = f (Xj), we 
could control quantity X; by controlling quantity X;. 


Let the signal from quantity X; be applied to the control or inspection devices during time interval t;. 


We shall assume that the system signals deviations of the observed quantities beyond definite limits. For 
the sake of generality, we shall assume that the deviations which give rise to signals are asymmetrically distrib- 
uted about the quantity's optimal value Xj9, and we shall denote them by 6X} and 5X". 


In addition, let the maximum admissible limits of deviations of quantity Xj be Aj and A‘. Then Xj will 
lie within the limits 
Nig — Ay < Xy < XyQ + AG. 
The quantities Aj are given by the conditions of satisfactory operation of the controlled object; in particular, 
they might be given by the conditions for accident-free or danger-free operation. 


If, at time t = 0, the quantity X; equaled Xj, then its possible values during time interval t are conditioned 
by the probability distribution function of tle random variable.X;. 


The probability that, with observation on regulation of quantity X; during time t, this quantity goes beyond 
the admissible limits is [ 1): 


PHP), ,1 pO, , pm . + 
(A;, A;) (BN; | ) (A;. A;) 


a,  ~. , Pe o£... @) 
(8X ;, 5X 5) (BN ;, 8X ;) (Aj. A;) 


Here, ve A%) is the probability that the magnitude of X; goes beyond the admissible limits before time 


© Rx, 5X") is the probability that the deviations of X; before time t will be less than the deviations giving 
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rise to signals, p(t) 


( at AN) is the probability that quantity X; goes out of bounds before time 2t, etc. 
io** jj 


We now elucidate formula (1). The circumstance that X; can go out of bounds before time t is expressed 
in formula (1) by the term ri As)’ However, during this same time, quantity X; can vary so slightly that no 
signal is forthcoming from it. In this case, X; will not be regulated, and will contimue to vary independently of 


the controller during time 2t. This case is expressed by the term MK” at) _, which is the probabili 
«ay é peril (8X},6X%) (A.A) seid 


that the quantity X; can go out of bounds during time 2t if it were not regulated at time t. 


Since probability P is less than one, the terms of the series will decrease rapidly. Therefore, for practical 
computations, we can limit ourselves to several of the first terms only. 


Now, let the maximum allowable probability that quantity X; go out of bounds be given. We denote it by 
Pmax: 
The magnitude of Pmyax depends on the nature of the quantity being controlled and its role in the process 
being controlled, The more heavily significant deviations of X; from Xj reflect themselves on the operation of 
the object being controlled, the smaller must P,,,, be. 
Thus, the determination of P,,,, depends on the properties 
o(%;) of the object and its individual parameters. 


It is obvious that P(") must be less than P.,,,. As the 
time interval t increases, the magnitde of p(t} increases. 
For some t, the equality ptt) = Prax Will hold. This t 
will also be the optimal control time tj. Indeed, increas- 
ing t makes it possible to devote a larger segment of time 
to the control process. On the other hand, the increase of 
ty is limited by Pmax- 


The method of determining t; is best illustrated by 
examples. 











For the quantity X;, let the probability distribution 
density be given in the form of curves, constructed for several times t (Fig. 1). 


Here, pit) is the curve of the probability distribution density for time interval ty, and pi") is the analogous 
curve for time t}. Then, 


Mgt Ai 
po, . =i- \ p? (X,) dX, (2) 
(Ay, Ay) 
Nig Ay 
Xigt 5X; 
PO, , = { (O(¥,) dX,, 
(®X;, BX)) Ri (4) ee, (3) 
Xig 5X; 
MigtA; 
pe) , =1-— \ (9 (X,) dX,. (4) 
(Ay A,) e # i) : 
X,,—A; 


10 1 


The integrals in formulas (2)-(4) are easily computed graphically from the areas bounded by the curves 
p(X;). 
For the determination of t;, it is necessary to find, on the basis of (1) and (2)-(4), the value of P* for sev- 


eral t. With Pyax known, the magnitude of t can be found by interpolation between the values of t,, to which 
the probabilities P\') correspond, one of which is greater than P,,,, and the other smaller. 























If the probability distribution densities are not given for the random variables X;, then one can quite simply 
determine the probabilities from the graphs of the distribution functions themselves (Fig. 2): 


ot) ‘ , . 2 . . 

Peaccny Xie t A) Fi Nig AD. (5) 
Mansangy = Ot N10 4 840. — Fe (Tuy 8X0 (6) 

pe > i { : F YX, ; 1 = ; ‘ Mts oat 

(Ay. AD [Pop (Xig + Ay) = Far (Nig i) (7) 


where F; and Fat are the ordinates of the probability distribution curves given, respectively, for t and 2t, 


It should be mentioned that it is completely nonobligatory to know the probability distribution density of 
the probability distribution function. In practice, it is necessary to construct, from experimental data, the graphs 
of several probabilistic characteristics of the quantities X; for various times t and, by interpolation, to find the 
optimal times t;. Expressing the probabilistic characteristics 
by functions is only convenient when these functions are 
known in the literature, or may be simply found, which 
rarely happens. 
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Thus, let the optimum time intervals for all the quan- 

: tities X4, Xq, ..., Xj, .-., Xp be found. We now show how 
clash f one sets up a program for signal transmissions, 
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We select a basic time interval ty, and express all the 
Fig. 2. other times, ty, ty, ..., tj, -.»» ta, im terms of it. It is neces- 
Sary to arrange matters so that the numbers expressing the tj 

in terms of tg will be mutually prime. With this, the first part of the program setup terminates. 


We Illustrate what has just been said by means of an example. Let there be three quantities to be regulated: 
X4, Xg, and Xs. For them we find the optimal intervals: t, = 3 sec, t, = 8 sec, and ty = 5 sec. Consequently, the 
first quantity must be sent every 3 sec, the second every 8 sec, and the third every 5 sec. 


Using this same example, we now consider the second part of the program setup. 


Thus, it is necessary to send quantity X, at times 3 sec, 6 sec, etc., quantity X, at times 8 sec, etc., and 
quantity Xg at times 5 sec, etc. If, at time t = 8 sec, we were to recommence the program at the beginning then, 
between transmissions of X;, there would be the time 2 + 3 = 5 sec (instead of 3 sec) and, between transmissions 
of Xz, time 5 + 3 = & sec (instead of 5 sec). 


If, now, the transmission of the quantities were programmed in accordance with the table 

















rf te | 3 | 5 | 6 | 8 | » | 40 | 12 | 15 | 16 | 17 
sec | 
Quantity sent | Vi Ns Ny Vs Ny Vs | Vy Vs | Y, Vs 
| 




















and, after this, the program were recommenced, then X3 would be sent after 7 sec, X, after 4 sec, and Xg after 
8 sec, 


As is clear from the example, it is desirable so to choose the moment of transition to the next cycle of 
the program that both the duration of the cycle, teng, and the duration of the neighboring steps, teng — te, 
tend + ts (and, in case of a large number of parameters, also t.,q — 2tg, etc.), be divided properly into cycles 
expressing the durations of the intervals tj. With this, the transition error is decreased, If t,,4 is made directly 
of severalt;'s, it isthen necessary to disperse the moments of transmission of the quantities. 


With this, of course, attention must be given to the physical nature of the X;. If an X; must be maintained 
accurately (for example, so as not to destroy the controlled object), it is then impossible to increase the corres- 
ponding tj, and some other t must be increased. 
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Technologically, this program can be implemented by means, for example, of a stepping switch with equal, 
or unequal, steps. 


Choosing optimal time intervals and introducing them into the regulation or observation program has sig- 
nificant advantages. 


If, as is ordinarily done, all quantities are checked successively with equal time intervals, then all the At's 
are equal among themselves. Here, At is the interval between signal transmissions from different quantities. Let 
T be the time of one observation cycle, and let us know the minimum admissible time intervals for the most 
rapidly varying i-th quantities. Then, T st, and At = T/n. 


If the optimum time intervals are found and are included in the program, then At = ty. 


Let us consider an example. Let five quantities be observed, where t; 4, = 2 sec and ty = 1 sec. Then, 
for the ordinary method of sending signals, At = 0.4 sec and, for the method proposed here, At = 1 sec. 


Thus, the incthod proposed here makes it possible to decrease the frequency of switching. It also eases the 
work of the operator and increases the length of service of the apparatus. 
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AN INSTRUMENT FOR DETERMINING THE FREQUENCY 
CHARACTERISTICS OF NONLINEAR SYSTEMS 


K.V. Zakharov and V.K. Svyatodukh 
(Moscow) 


The theory of operation, as well as the device itself, are described for a new 
electronic instrument which permits the aforementioned frequency characteristics to 
be determined with high accuracy in the frequency range f = 0.25 to 50 cps, directly 
in the experimental process. 


For the theoretical designs and dynamics simulations of automatic control systems, an ever increasing use 
is being made of the experimental determination of the frequency characteristics of the systems themselves and 
of their individual elements. However, whereas the methodology of experimental determination of the frequency 
characteristics of linear systems is developed in great detail, and many works devoted to this question have been 
published, these is still no satisfactory method for determining experimentally the frequency characteristics of 

















Fig. 1. 


nonlinear systems which, to a definite degree, inhibits the use of frequency characteristics for the analysis of non- 
linear systems. The difficulty which arises in the experimental determination of the frequency characteristics of 
nonlinear systems amounts to this, that the forced oscillations of their output signals due to a harmonic stimulus 
Up = AgSinwt at the input contains, not only the fundamental frequency, but also higher harmonics, and may be 
given in the form 
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Since, as the frequency characteristics of nonlinear systems, the amplitude and phase characteristics of the 
first harmonic, uy = Aysin(wt — 4), are used (cf. for example, [1,2]), the problem reduces to the determination 
of the amplitude A; and phase y of the first harmonic. 


In conjunction with this, the necessity arises for estimating, at least roughly, the higher harmonics of the 
output signal, so that one may judge the possibilities of using the frequency characteristics found for the analysis 
of the dynamics of a system which contains the nonlinear element investigated, and also so that one may gauge 
the possibilities of linearizing the equations of motion of these nonlinear elements [3, 4]. 


Although the estimates cited can, in theory, be given by means of the well-known methods of harmonic 
analysis [5,6], the practical use of these methods, which allow both the amplitude and the phase spectra to be 
determined, is arduous since, in addition to the large amount of computational work required by these methods, 
a preliminary recording (oscillographs) of the process under investigation is also required, which reduces the ac- 
curacy and speed of obtaining results. 


Recently, attempts have been made to create special instruments for determining the frequency character- 
istics of nonlinear systems. In particular, the electronic instrument described in [7] permits the experimental de- 
termination of the amplitude and phase of the first harmonic of a periodic signal in the frequency range f < 1 cps, 
by automatically calculating the first coefficients in the Fourier series, 


We give below a description of an electronic instrument, developed by the authors, by means of which one 
can, without resorting to oscillographs and complicated computations, carry out the harmonic analysis of the out- 
put signal of a nonlinear seryo-system whose input and output signals are in the form of electrical voltages in the 
frequency range f = 0.25 to 50 cps. The instrument described can be used, in particular, for taking off the fre- 


quency characteristics of servo-systems and also of power-drives of afitomatic control systems. The over-all view 
of the instrument is given in Fig. 1. 


Theory of Operation of the Instrument 








The theory of operation of the instrument is based on the compensation (of the harmonic being investigated) 
of the nonlinear system's output signal up, = A, sin(nwt ~ p,) by a specially formed signal u* = A*sin(nwt — $*) 
whose amplitude and phase are separately controlled, and whose frequency equals the frequency of the harmonic 
being investigated. The fact of coincidence of the amplitude and phase of the compensating signal u* with the 
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Fig. 2. 1) Sine wave generator; 2) nonlinear system under investigation; 
3) scale block; 4) forming device; 5) comparison block; 6) controlled 
selection system; 17) indicator. 





amplitude and phase of the investigated harmonic is established by the absence, in the difference Au = u ~u*, 
of the harmonic of frequency nw, which is accomplished by applying the signal Au to the input of a controlled 
selection system which is tuned to the frequency nw of the harmonic being investigated, and which has a suffi- 
ciently narrow pass band. 
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Let us consider the operation of the instrument in the case when one wishes to determine the amplitude 
and phase of the fundamental harmonic of a nonlinear system's output signal (the block schematic of the device 
is shown in Fig. 2). The electrical signal ug = Agsinwt from the sine wave generator is applied simultaneously 
to the input of the nonlinear system being investigated and to the input of the forming device. In the forming 
device, by means of a block for tuning the phase and a block for tuning the amplitude, there are implemented 
individual transformations of the phase and the amplitude, as well as variations in the sign of the signal ug. The 
output signal of the forming device, u* = A*sin(wt — ~*) and the investigated signal u, amplified by the scale 
block, are applied to the comparison block. 


After this, the output signal of the comparison block 


co 
Au=u—u* = A, sin (mt — v,) + > A,, sin (nw/ —*,), 
n=2 


where 


Ay = Ay 125 cos (hr — ¢°) + (F-)» 9a = are lg Gente — a can 
is applied to the controlled selection system, which is intended to separate the fundamental harmonic and to 
suppress the higher harmonics. From the amplitmde of the oscillations of the output of the controlled selection 
system, All ~ A,sin(wt — ~,), one can judge how well the amplitude and phase of signal u* coincide with the 
amplitude and phase of the fundamental harmonic of signal u: for Ay = A* and y = ~*, the amplimide of the 
selection system's output signal will be minimal. In particular, if the selection system does not pass the higher 
harmonics at all, or if the investigated system is linear, then the amplitude of the controlled selection system's 
output signal, for Ay = A® and y, = ~*, will equal zero. For observing the output signal of the selection system, 
we use an indicator whose natural frequency is significantly higher than the frequency being investigated. From 
the output signal af the comparison block, Au, with total compensation of the fundamental harmonic, one can 
make a summary estimate of the higher harmonics. 


For obtaining the amplitude and phase of higher harmonics, one must provide for the possibility of obtain- 
ing, together with the harmonic signal of the fundamental frequency, the signals at multiples of this frequency. 
If a sine potentiometer, turned by a drive motor, is used as the transducer for the fundamental frequency then, to 
obtain multiples of this frequency, one can also use a sine potentiometer, one whose shaft is geared to the shaft 
of the potentiometer for the fundamental frequency through a reducer which would provide a rotational speed 
which is an integral multiple of the speed of the drive shaft. 


If the system being investigated is linear, its frequency characteristics can be determined by means of our 
instrument, not only for harmonic input signals, but also for such periodic signals as rectangular impulses, saw - 
tooth signals, etc., which follows from the principle of isolating the first harmonic which underlies the operation 
of the instrument. 


Theory of Operation of the Instrument's Basic Blocks 





The forming device is a series connection of a phase-shifter and a voltage amplifier (Fig. 3). As is ob- 
vious from the forming device's transfer function: 


ay oe 


(p's = K ephtnemeny 
u (Pp) a( é re Be 


‘ 5 
>’ ee j= 


where T = RC is the time constant of the RC network and K, is the over-all gain, its amplitude characteristic 
A*/Ag = Kg does not depend on either the frequency w or the time constant T. The phase lag y*, defined by 
the formula tan(~*/2) = wT, depencis only on time constant T and frequency w. This allows the amplitude A* 
and the phase ~* of the compensating signal to be controlled separately, which significantly simplifies their 
tuning. 


The value of the amplitude characteristic A*/A, is determined directly by indicator divider Dy, with the 
gain K, of vacuum tube amplifier VTA, taken into account. This value equals Ky = aK,, where a is the scale 
coefficient for divider Dy. (For convenience, the circuit is provided with two possible gains, Ky = 1 and K, = 10.) 
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The value of time constant T, necessary for computing the phase shift ~* = 2arctanwT, is determined by 
the scale indication of variable impedance R and switch Ky for capacitors Cy and 10C;. Switching of the capa- 
citors provides a widening of the range of variation of time constant T. 


For variations of w from 0 to oo, the phase shifter provides phase shifts in the range 0 = p* < 180°. There 


is no difficulty, in principle, in creating still greater phase shifts, e:ther positive or negative, by varying the sign 
of the signal being investigated. 
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The scale block is the conjunction of a vacuum-tube amplifier with a voltage divider, designed to bring 
the scale of the nonlinear system's output signal to the scale of the input signal ug. The inclusion of the scale 
block in the scheme simplifies the operation, since the value of amplitude characteristic K,, after being thus 
brought to the proper scale, can be determined directly from the scale of divider Dy (with account taken of gain 
K, of vacuum-tube amplifier VTA). The amplifier in the scale block also provides the transformation of the 
voltage of the output transducer of the system under investigation to the normal working voltage, which is par- 
ticularly important in case one is measuring a mechanical output biasing by méans of tension transducers, or by 
transducers with low output voltages, In addition, a voltage from a controlled potentiometer is applied to the 
input of the scale block amplifier for the compensation of the de component of the investigated signal. 





Selection system. As the selection system in the instrument being described here, we use a vacuum tube 
amplifier with a twin-T bridge in the feedback path (Fig. 4). Such a circuit is a. band-pass filter whose passband 
is defined by the magnitudes of resistors Ry, Rg, and Rg, capacitors Cy, Cy, and Cg, and gain Ky [8]. For a suffi- 

ciently high degree of selectivity (i.e. for a sufficiently narrow 
passband) this system, tuned to the frequency nw of the har- 
l vr As | monic being investigated, does not pass, for all practical pur- 
poses, any higher or lower harmonics, which permits one to 
judge, from the absence of oscillations at its output, that the 
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ye! zy harmonic under investigation is absent from its input signal. 
? th rey Tuning of the selection system to the investigated fre- 
, } | quency is carried out by simultaneous variations of resistors R,, 
7 pp die Rg, and Rg in such manner that the ratios Ry: RgsRg = 121:0.5 
| —t 1G Y ’ j are kept. With these ratios of the resistances, and with capacitor 
oe——--—-— 2 


ratios of Cy:C_g:Cg = 1:1:2, the frequency passed by the sys- 
Fig. 4. tem, defined by the feedback cut-off frequency, equals ws, = 

= 1/RyCy. For the use of the selection system, a tenfold change 

of we was provided for by the capability of switching condensors 
Cy, Cy, and Cs. To ease the computation involved in determining the phase shift ~* = 2arctanwT, the scale of 
Ry was graduated in units of angular frequency. The selectivity of the system described depends on the amplifier 
gain without feedback, Kg, and will be proportional to Kg. However, an increase in Ky protracts the transient 
response when the amplitude and phase of the input signal change. Since the duration of the transient response 
turns out to be proportional to the quantity Ky + 1, and inversely proportional to the cut-off frequency ws of the 
feedback, it is not practical at low frequencies to use a large value of Ky. In the selection system being described, 
it was made possible to increase the selectivity by a transition from low frequencies, f < 2 to 5 cps, to 
higher frequencies, f > 2 to 5 cps, by varying gain Ky in the range 5 = Ky $25. The given values of Ky provide 
variations in the suppression coefficient of the third harmonic, A(w)/A(3w), in the range 3.4 = A(w)/A(3w) = 16,4. 
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For the determination of the frequency characteristics of a number of nonlinear systems (boosters, steering 
gear, and other power drives), in which the amplitude of the third harmonic does not exceed 10 to 20% of the 
amplitude of the fundamental, it is ordinarily sufficient to have the coefficient of suppression of the third har- 
monic of the order of 5 to 8. In this case, instead of an amplifier with a twin-T bridge in the feedback path, one 
can also use a second-order oscillatory system with a low damping coefficient, on a low-frequency filter in which 
the cut-off frequency lies between the fundamental frequency and its triple. Thus, if the oscillatory system has 





the transfer function ©(p) = re A ~— a1! it will suffice, for attaining the sensitivity mentioned, that the 
condition ¢ = 0.5 to 1 be met. For ¢ = 1, we have the low-frequency filter &(p) = K/(Tp + 1). 


Sequence of Operations in the Determination of Frequency Characteristics 





1, Tuning to the scale of the nonlinear system's output signal. For purposes of scale-tuning, a low-frequency 
signal is impressed upon the system's input, for which the system's amplitude characteristic may be taken equal 
to unity and then, by varying the amplifier gain and the division coefficient of the scale block, equality of the 
input and output signal amplitudes can be reached. For this, the gain K, of the forming device must be unity. 





2. Tuning of the selection system. The selection loop used in the instrument for controlling both frequency 
and degree of selectivity can also be simultaneously used for measuring the frequency w. In case the frequency 
w of the forced oscillations is not exactly known, it is necessary to tune to the frequency w by changing variable 
resistors Ry, Rg, and Rg, while observing the results of this tuning by the amplitude of the selection system's output 
Signal. The frequency measured by the scale of Ry can be used for an approximate calculation of the phase shift 
~* = 2arctanwT. For the accurate calculation of the phase shift, it is necessary to determine, with the requisite 
accuracy, the transducer's fundamental frequency by means of a special instrument, namely, a freqmeter. 





3. Tuning the amplitude A* and the phase ~* of the compensating signal. Tuning of compensating signal 
u® should first be implemented by phase, and, thereafter, by amplitude, since this significantly shortens the time 
necessary for the tuning. The necessity for this sequence derives from the following. 
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The result of tuning the compensating signal is observed by seeing a minimum amplitude of the selection 
system's output signal Aii which, in practice, is proportional to the amplitude A, of the difference 


uy —U* == Ay sin (wt -— ty) ~ A* sin (ml - ) : Ay sin (ot — Y,)- 


To provide the foundation for the tuning sequence given, it suffices to consider the effect of A® and p* on 
the amplitude of 


CS RT Fo 
A, = Ay Vi 2 mn cos Ad + (5) (Ay = +, --- $*). 


As is obvious from the expression for the partial derivative 
OA,g A* sin Ay 


(6%) V1 2 cos A ++ (5-) 
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the minimum of the amplitude A, as phase p® varies (for a 
fixed value of A*) is reached for Ap = 0, Le. for the coinci- 
dence of the phases ~* and #, independently of the magnitude 
of amplitude A*. At the same time, it follows from the ex- 
pression for the partial derivative 





A 
0A, (4 — cos Ay ) Ay 


(4) ‘ V { — 27 cos ay + (7) 


that the minimum amplitude of A, as A® varies is attained 
at the value A*, defined by the equation A*/A; — cos Ay = 0, 
i.e. depending on the phase difference Ay. 














Therefore, the tuning of the amplitude must be done 
only for a tuned phase, However, if the amplitmde of the 
compensating signal differs significantly from amplitude of 
the harmonic being investigated, then the curve A, = f(A) 
has a very shallow minimum (Fig. 5), which makes an exact 
tuning of ~* difficult. Due to this, the uming of ~* and A* 
must be carried out in two steps; first a rough tuning of phase 
and amplitude and then, in the same order, a refined (Vernier) 
tuning. 





As an example, Fig. 6a and 6b shows the oscillograms 
Fig. 6. of the signals up, u, u*, Au and Aji, obtained in the determina - 

tion of the frequency characteristics of a nonlinear system 

with an odd static characteristic, du/dt = f(u ug). Figure 6a 
shows the signals when there was complete coincidence of the compensating signal u® and the first harmonic sig- 
nal u, and Fig. 6b shows them when these signals were mismatched in phase. As is clear from Fig, 6a and 6b, the 
output signal Ali of the controlled selection system does not, in practice, contain higher harmonics so that, from 
the amplitude of its oscillations, one can easily observe the fact of a mismatch of signal u* with the first har- 
monic signal u. 


Brief Analysis of Instrument Accuracy 





The accuracy of the determination of amplitude and phase characteristics using the instrument herein de- 
scribed is affected by the following factors: 


1, The relationship between the amplitudes of the fundamental and the higher harmonics of the nonlinear 
system's output signal, characterized by the coefficients a, = Ay/Ap. 


2. The values of the suppression coefficients of the selection system 


A(o) 
A (nw) : 





B, = 


where A(w) is the amplitude of the selection system at its resonant frequency. 


3. The accuracy of the visual observation of the fundamental harmonic in the selection system's output 


n 


A 
signal Ai, characterized by the coefficients +,, — 7-=— 7: for which the presence of the fundamental 
A n 4 


harmonic in signal Au is still sufficiently readily observed (here, we denote by Ae and An the amplitudes of the 
fundamental and the greatest n-th harmonic at the selection system's output). 


4. The accuracy of determining the original frequency w. 


5. The accuracy of establishing the amplitude A* and the phase ~* of the compensating signal provided 
by the forming device. 
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u“ 
'\ [/\ The error arising from the first three factors is determined by the 
a ; accuracy of compensation of the fundamental harmonic, and will 
decrease proportionally to the ratio AA/Ay. Indeed, in case of coin- 
LA/ cidence of phase (y~* = y,), the relative error of amplitude determina - 
tion is | Ay — A*| /Aq = A,/Ay For amplitude coincidence (A* = Aj), 
Au the amplitude of the first harmonic of signal Au, for small values of 
Ay, can be put in the form A, = A| Ay|, from which it follows that 
b VU Av) AF t the error in the determination of phase is | Ay| = A,/Ay From the 
definition of coefficients ap, Bp, and yp it follows that if, in signal 
Aa Aili, the basic harmonic observed for y,, is not less than some value v. 
then, for the relative amplitude A,/A,, we have the estimate A,/A; « 
[\ [\ « ManBnYp- Analogous estimates will hold for the amplitude and 
: phase: 
V 
UV (As AJP Tl pape —_* 
Ay ae * Se 
Fig. 7. a) Investigated signal 
(An>g = 0) and u = Aysinwt + For the majority of nonlinearities arising in practice, apn min > 2 
+ 0.5 Aysin3wt, ag = 2; b) input (for example, for an element with a relay characteristic, ap min = &3 © 
signal to the selection system (dif- = 3). Therefore, for an estimate of accuracy, we can take ap min = 2: 
ference signal) Au = 0.007 The minimum value of the coefficient y, is determined, generally 
A,sinwt + 0.5 A,sinwt, Aa/Ay = speaking, by subjective factors, since it depends on the operator's train: 
= 0,007; c) output signal of the ing. Figure 7 shows the example of a signal Au for ys = 5 (An>s3 = 9), 
selection system Aii = k(0.2sin wt + from which it is obvious that the presence of the fundamental harmonic 
+ sin3wt), By = 14.4 and y3 = 5. in signal Au when yz = 5 is observed quite clearly. In the future, for 


estimating accuracy, we shall take yj, = 5. 


Suppression coefficient 6,, is defined by the amplitude characteristic of the selection system. For the 
chosen selection system [8] 


+ (#5) 


‘ 





t bo 2 
(Ke-+t* + (75) 


where w is the frequency relative to resonant frequency wp. 


It follows from this that 


For Kg = 25 and n = 3 we have that Bs = 15.4. By substituting the values agmin = 2, yg = 5, and Bg = 16.4 
in the expressions for the estimates of | Ay — A*| /A; and | A], we find that the maximum error due to the first 
three factors does not exceed 0.7% for amplitude and 0.4° for phase. 


The error in defining the initial frequency Aw affects only the accuracy of the phase determination. As 
follows from the expression for the phase shift, » = 2arctanwT, the error in phase determination is Ay = 
7 
~ Aw = 7 ene Aw or Av = Awsing/w. The maximum error occurs for » = 90° and equals Av max * 
= Aw/w; from whence it follows that the error does not exceed the relative error in frequency measurement Aw/w. 
Existing frequency meters can provide frequency measurements with an error of less than 1%. With this, the error 
Ap will not exceed 0.5°. 


The accuracy of establishing the phase ~* and amplitude A* provided by the forming device is determined 


1636 





or 





by the accuracy in reproducing the transfer function K,(Tp—1)/(Tp+1) by means of the chosen electrical system, 
and also by the deviations of the actual parameters from the design values. 


Detailed investigations, the presentation of which is not possible in this paper, show that the maximum er- 
ror due to the forming device cannot exceed 0.5" in phase and 1% in amplitude. 


From what has been presented, it follows that the total maximum error given by the instrument cannot, in 
practice, exceed 1 to 1.5° in phase and 1.5 to 2% in amplitude. 


Experimental data have backed up the data given here on the accuracy of the instrument's operation. 


Today, the instrument described is used for determining frequency characteristics of actual automatic con- 
trol systems and their elements in a number of undertakings. 


In conclusion, the authors wish to express their gratitude to G.V. Aleksandrov for his aid in building the in- 
strument, and also to A.F. Vasilevskaya, Yu.L Rudakov, and K.V. Pavlov for taking part in the development of 
the individual blocks of the instrument. 
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DIGITAL MODELS 


A.V. Shileiko 
(Moscow) 


A systematic survey is made of contemporary literature dealing with a class of 
new computing devices, namely, digital differential analyzers, digital computers op- 
erating by the method of summing increments, and with individual function generators 
whose input and output quantities are presented by the method of delta~modulation of 
pulsed electrical voltages. The theoretical bases of the operation of digital models is 
formulated, and comparisons are made of their characteristics with the characteristics 
of electronic digital computers and of analog computers. For this, a computing ma - 
chine is considered as an equivalent filter which transforms an input signal, and as the 
criterion for comparing different classes of machines, we take the product of the fil- 
ter’s resolving power by the width of the pass-band spectrum. 


Recently, there has been heightened interest in a new type of electronic computer, called digital differen- 
tial analyzers, or machines which operate by summing increments (incremental computers). The fundamental 
characteristic peculiarity of these machines is the possibility of isolating in their make-up individual functional 
blocks which execute a mmmber of basic mathematical operations, such as integration, summation, multiplying 
by a constant factor, etc. For the make-up, or programing, for a problem, the blocks are interconnected in ac~- 
cordance with a scheme which provides a given sequence of operation executions. From the user's point of view, 
the digital differential analyzers occupy a place intermediate between the universal (general-purpose) digital 
machines and analog computers, since they join the potential accuracy of digital computers and the simplicity 
of setting up for a problem characteristic of analog computers. 


Independently of digital differential analyzers, a number of devices have recently been developed, and 
described in the periodical literature [1-24], which employ the methods of presenting variables which have been 
adopted in digital differential analyzers,and which are designed for the execution of individual mathematical 
operations, such as multiplication, generating of sine functions, etc. The essential fact here is that such blocks, 
individually or in definite combinations, can be used for simulation in the true time scale of definite physical 
processes and, consequently, can operate as elements of automatic control systems, This allows us to unify them 
under the common denotation of “digital models," and to formulate certain theoretical bases for the operation of 
such devices, 


In digital models, the method of presenting quantities is that known in communication theory as delta- 
modulation.* In accordance with this method, only the increments of the original variables take part in all the 
transformations, and only at the final step of these transformations is the resulting quantity obtained in the form 
of a sum of the individual increments, The scales of all quantities taking part in the transformations are quan- 
tized, and the magnitude of an elementary increment is equated to the quantum step 6-. Elementary increment 


*It is specifically the delta-modulation method which is used, and not a number-pulse code which, in the gen- 
eral case, can be multipositional (cf. for example, A.A. Kharkevich, Outline of General Communication Theory, 
Gostekhizdat, 1955). 
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are presented, in digital models, by individual electrical pulses. The time scale is also quantized and, thus, the 
presentation of increment pulses can only occur at definite moments of time, to, ty, ... t, separated by the step 5, 
of the quantized time scale. The absence of a pulse at a given moment of time t; has, in one case, the meaning 
of a negative increment, ~65 (a binary incremental system) and, in other cases, the meaning of zero (the ternary 
incremental system). For presenting the sign of an increment in the ternary system, an additional pulse, trans- 
mitted along a second channel, is used. 





a,, 4,(¢,) »— 
Oy 4,(¢;) aera f a, 5 (t,) 
(4) >— 





fT 4, 





Fig. 1. 


Any information transformer which operates with quantities presented by the delta-modulation method re- 
ceives, at its inputs, successions of increments of the form 


@59 5x (t,), 4 By (by), oo Og 8. (ts. 
By By (to). yy By (ty)e-- +s Aig By (tis «+s (1) 
Bg Sy (t,)s Bay Bay (by )o - + = 9 Fagg Bay (Eg) s oe + 


where the subscript on a 6 indicates the quantity being presented, and the coefficients a,; can take the values 


of +1 and —1 in the binary system, and +1, 0, and —1 in the ternary system. At the information transformers 
output there is also formed a series of increments, of the form (Fig. 1) 


Brg bz (ty)s O48, (ty)a- - «+ %yi 8, (ty (2) 


The operation of an information transformer can be described, generally speaking, by the corresponding re- 
lationship between the original quantities, x, y, ..., w, and the resulting quantity z: 


=m f(z 9, .--, w). (3) 


Here, however, there is a definite inconvenience, related to the fact that the “instantaneous” yalues 
Xj, Yj» ---, Wj and 2; are obtained as the result of the accumulation of increments during the entire time of opera - 
tion of the transformer and, consequently, the relationships between these values are integral characteristics, and 
do not describe the operation of the transformer at each given moment of time. For describing transformer opera - 
tion, it is more convenient to use some linear function of the x, y, ..., Ww which does not possess the disadvantage 
just noted. Ordinarily, as such functions, we choose the derivatives with respect to time of the corresponding 
quantities, characterizing the speed of their change and, to describe the operation of information transformers, 
we use relationships of the form 


dz dx dy dw 
ai "ae dade): “ 


For quantities presented by the delta-modulation method, where quantized amplitude and time scales are 
used, approximate values of derivatives with respect to time are defined as 
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The step 5, of the quantized time scale is, as a rule, the same for all the transformers involved in the solu- 
tion of a given problem, The quantization step 6, is chosen for an established scale of the variables. With these 
conditions, the approximate normalized value of the speed of change of variable & can be defined by the expres- 
sion 
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Blocks Executing Individual Mathematical Operations 





The simplest form of digital devices which execute individual mathematical operations are the devices 
which realize relationships of the form 


dz dx 
—_ 2 7 
dt K di’ (7) 


where K is a constant. For the case when the increments are only positive (the variant of the binary system of 
increments in which the coefficients can take only the values of 0 and 1), but K = 0.5, the role of such a trans- 
former can be played by an ordinary trigger (flip-flop). Analogously, a series connection of two, three, etc. 
triggers can realize relationship (7). for the cases when K = 0.25, 0.125, etc. Using triggers, we can construct cir- 
cuits which realize relationship (7) for any K of the form 


K = 1/n, (8) 


where n is an integer. A circuit for n = 1, 2, ..., 16 is shown in Fig. 2. Here, we have used the convention that 
the 0 state in each of the four triggers (denoted by Tr on Fig. 2) corresponds to the conducting state of the right 
vacuum tube, and state 1 corresponds to the conducting 


7" M Ky wie Ne state of the left vacuum mbe. Excitation of each follow- 
| i,” | ! j ! 2 
as wa 3 ing cell is produced by negative pulses (i.e. at the mo 
























































in | bd ment of triggering a tube, whose anode is joined to the 
— feedback path). The output pulse, via the feedback path, 
beonk resets all cells to zero. The fact that any of the keys 
(contacts) K,~-K, is in its lower position is described by 
Fig. 2. the symbol 1, the symbol 0 being used when a key is in 


its upper position. 


The combination 1111 of key positions corresponds to the standard circuit of inter-stage connections for a 
binary counter. If we assume that a continuous sequence of electric pulses is applied to the circuit's input, then 
we shall have the following sequence of combinations of trigger states: 0000, 1000, 0100, 1100, ..., 1111. The 
negative output pulse will appear each time the circuit makes the transition from state combination 1111 to state 
combination 0000, ie., once for each 16 input pulses (n = 16). The combination 0000 of key positions would, 
with no feedback path, correspond to the sequence of combinations of trigger states: 0000, 1111, 0111, etc. How- 
ever, thanks to the feedback path, each input pulse gives rise to the appearance of an output pulse, and returns 
the circuit to state 0000 (n = 1). With the feedback path switched out, the combination 1100 of key positions 
would correspond to the sequence of trigger states: 0000, 1000, 0100, 1100, 0011, etc. With the feedback path 
closed, an output pulse returning the circuit to its initial state would appear at the moment when combination 
0011 appeared. This corresponds to the case n = 4. From a further consideration of the circuit of Fig. 2, we can 
immediately draw the conclusion that such a circuit, with an arbitrary number of triggers, will realize relation- 
ship (7) with the condition that 
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Fig. 4. 





m 
n==4 +4 >) b,-261, (9) 
{=1 


where the b; are coefficients assuming the values of 1 
or 0 and which characterize, in the notation adopted 
earlier, the states of the corresponding keys. Circuits 
of this type are described in more detail in [6]. There 
is a description, in [4], of a circuit for generating sinus- 
oidal functions which uses an analogous method. 


Another variant of the circuit for realizing rela - 
tionship (7) for any K is shown in Fig, 3. Here, there 
is again used a series connection of triggers, Try—Trp. 
If an arbitrary sequence of pulses, characterized by 
velocity v [the velocity is defined by formula (6)], is 
impressed on the circuit's input then, at the output of 
each trigger, there will appear a sequence of pulses 
with velocity 


“= (i me 4,2; << «itd, (10) 


where i is the ordinal number of the trigger. For this, 
the output pulses are connected to the anodes of the 
corresponding triggers in such fashion that no pair of 
pulses from any pair of triggers coincides in time (Fig. 
3b). If this latter coudition is met, then the velocity 
of the pulses at the circuit's output will be defined by 
the relationship 


Vv 

1 —j 
v = c,0-2 11 
out 24 at) 


where the coefficients c, characterize the states of the corresponding gates G; (a value of 1 for a coefficient c 
corresponds to an open gate G;, while the value of 0 corresponds to a closed gate), and n is the total number of 


triggers. Here, consequently, 


n 
K= >) oa. (12) 


In both circuits considered, the quantity K can, generally speaking, be some function of the quantity x: 


In this case, relationship (7) can be rewritten in the form 


K = K(x). (13) 


dz = K(x)dx (14) 


een ( K (2) de. (15) 


Thus, the circuits described are blocks for multiplication by a constant, when K is constant (for K = const, 


it follows immediately from (7) that z = Kx + 2), and integrators when K = K(x). The case is also possible when 




































K is a function of the independent variable x, but varies in a stepwise fashion in accordance with a definite pro- 
gram as x reaches the given values x4, Xg, ..., Xj. Under these conditions, the transformer's output z will take the 
form of a broken line, and the transformer itself, in conjunction with the corresponding program device, can be 
used for generating arbitrary functions of the variable x by the method of piecewise-linear approximation. 


Figure 4 shows the circuit of an information transformer which operates in the normal binary system (the 
system of +1 and ~—1) and which realizes relationships between the original variables of the form 
i 
¥; =, + \ oY at, 


te (16) 
il = Ayi de 

dt 
dt 


Here, systems Y and R are delay lines or shift registers (SR, and SR) whose constants are numbers, coded in binary 
(or any other number system), which are continuously circulating through the circuits closed through the one- 
place adders Ly and £,. At the circuit's input denoted by 5, there are impressed clock pulses, which mark the 
successive moments of time tp, ty, ..., tj. The fact of coincidence of the current 5, pulse with a pulse at input 
dy/dt (ay{ + 1 in the binary system) is fixed by gate G,, and the output signal of this gate is applied to the input 
of adder £,, incrementing the contents of register Y by a unit in the least significant place. The fact of nonco- 
incidence of a 5, pulse with a pulse at the dy/dt input is fixed by gate Gy, which gates the number 11...10 (the 
least significant order is on the right), the one’s complement of a unit in the least significant bit, from comple- 
ment former CF, to the input of adder Dy. The system just described operates in the mode of a reversible counter, 
accumulating the current value of y;, and realizing the first relationship in (16). Naturally, the frequency with 
which the orders follow each other in the register must not be less than the quantity 


f=N/6,, (17) 


where N is the capacity of the register. 


The current value of y; from the output of adder Ly is applied to the input of adder £,, either directly 
through gate G,, if the current increment of variable x is ay; = +1, or by its complemented code from comple- 
ment former CF2 if the current increment of variable x is a,j = —1. In the first case, y; is added to the quantity 
contained in register R and, in the second case, is subtracted from it. If the current value rj exceeds the capacity 
of register R, which is equal to the capacity of register Y, then at the output, dz/dt, of the circuit there appears 
a pulse which is simply the overflow pulse of register R. If the capacities of registers Y and R have the common 
value N, then the quantity y; can be given in the form 


y, = 2N/m, (18) 


where m = 1, since y; =2N, It is clear from (18) that, to obtain one overflow signal from register R, it is neces- 
Sary to transmit the contents of register Y m times to register R. By taking into account that each such tansmis- 
sion is executed under the stimulus of the current increment of the quantity x, we obtain the relationship 


which is the second relationship of (16).if we take A = 2"N. Relationship (18), the validity of which for positive 
increments and positive values of y; and rj is obvious, will also be valid for any rj, yj, and a,j if negative values 
of rj and yj are presented in the complemented code, and if the most significant bits of rj and yj; (this position is 
made binary even when the numbers 7; and y; are coded in some other number system) play the role of sign bits, 
where a 1 in this bit corresponds to plu:. This statement can be verified by a consideration of the corresponding 
examples (cf., [1]). In an analogous way, one can construct a circuit operating in the ternary incremental system. 


The circuit of Fig. 4, just as the previous ones, can be used for multiplying the speed of variation of quan- 
tity x by a constant factor, for integrating and fot generating functions by the method of piecewise-linear approx- 
imation. Moreover, the circuit can be used for one more special operation. Since, in the given case, the over- 
flow signal from register Y is not used, but the presence of such an overflow might, naturally, lead to an erroneous 
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result, measures are usually taken to stop the computer if the Y register overflows. However, if such a step does 

not, for some reason or other, occur, then an overflow of the register by one unit in the least significant bit leads 

to the replacement of the maximum possible value of y, which was stored in this register prior to the reception of 
the unit in the least significant bit which engendered the overflow 
(11.,.1 in the binary system), by its minimum possible value 


P tnsgnesiaapsingilteisin (00...0 in the binary system). Conversely, subtraction of a unit 
x 4 in the least significant bit from the minimum possible value of y 
al leads to the formation of its maximum possible value. Due to 
4 (zy) ; this property, the circuit of Fig. 4, with the condition that the 
at minimum or maximum value of y is initially placed in register 








Y, can operate as an independent relay amplifier, in which a vari- 
iv! ation in the contents of register Y by one unit in the least signifi- 
y ay cant place leads to a change in the velocity of the output signal 
Fig, 5. from plus the maximum to minus the maximum. Ordinarily, in 
this case, a signal with the maximum velocity (signal 5,) is ap- 
plied to input dx/dt of the circuit. Such an amplifier is used for 


forming discontinuous functions, for setting up adder circuits and, in a number of other'special cases (for more 
details cf. [1]). 








For the multiplication of two variables, the circuit shown on Fig, 5 is used. Multiplication is implemented 
in accordance with the relationship 


A (zy) = why + yAa + AyAs = wdy + (y + Ay) Az. (19) 


The letters X and Y denote adding registers, analogous to adding register Y on Fig. 4, in which the current 
values of x and y are formed by the formulas 


tj 


x= Zo+ \ rie dt. 
te 
Fi (20) 


dy 
¥s= Yor \ Ty dt. 

Letter R denotes a summing register which is also analogous to register R of Fig. 4. In each cycle of circuit 
operation, quantity x is applied to register R (by the direct or complemented code, depending on the value of ay,) 
until the current increment (denoted by x; on the circuit) has been added to it, and the quantity y is applied to 
the Y register after the current increment has been added (yj,; = yj + Ay). The overflow signals from regis- 
ter R are the increments of the product, in correspondence with (19). 


The Construction of Digital Differential Analyzers 





For the solving of problems, the individual functional blocks of a digital model are interconnected in cor- 
respondence with the initial equations in exactly the same way as in electronic analog computers. There are two 
basic types of digital differential analyzers today: parallel and serial. A parallel digital differential analyzer 
is a set of an arbitrary number of functional blocks, each of which constitutes an independent physical package. 
Such a machine contains a common power supply and a common source of basic synchronizing pulses. For a 
problem setup, the blocks are interconnected by cables. 


In serial digital differential analyzers, the contents of registers Y and R of each integrator are written suc- 
cessively in common delay lines with capacity 


Q = MN, (21) 


where N is the capacity of one integrator register and M is the number of integrators in the machine. The circuit 
of a series digital differential analyzer is fundamentally analogous to the circuit of Fig. 4, with the difference 
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that the capacity of the registers is defined by (21). Here, at each moment of time, the operations in some one 
integrator are carried out, so that, for the execution of the operations in all the integrators, a time interval of 


8, = Mr, (22) 


is required, where r is the time of execution of one elementary operation in one integrator. The overflow signals 
obtained at the end of each elementary operation in each integrator are written in a special register, and are di- 
rected from there to the inputs of other integrators in accordance with a program, introduced into the machine 
before initiation of problem solution and also stored in special registers, Ordinarily, the tracks on a magnetic 
drum are used as the registers of series digital differential analyzers. 


The basic data on the existing types of digital differential analyzers are given in the table. For compari- 
son purposes, we have included in this table the data on two general-purpose digital computers of the correspond- 
ing classes, 


Accuracy of Digital Models 


The question as to the accuracy of digital models is still insufficiently studied. The basic error sources in 
solving a problem on a digital model are, naturally, the quantization errors. Some decrease in these errors is at- 
tained by operating in the ternary incremental system, and also by the replacement of approximate integration 
by the rectangular formula (16) by integration by the trapezoidal formula 





ds = A (y+ > Ay) de. (23) 


Moreover, the accumulation of quantization errors in the execution of operations on two variables in a digi- 
tal model depends on the phase relationships between successive pulses representing each of these variables. This 
phenomenon is called phase errors. The process of phase error formation can be elucidated with the following 
example. In the register of some integrator, let there be stored zero which, in the corresponding mumber system, 
will have the form 100... At inputs dx/dt and dy/dt of this integrator, let there be applied signals of the form 
0, 1, 0, 1,..., 0, 1, with zero velocity. The contents of register Y will change in the following sequence: 100, 
++, O11, ..., 100, ..., etc. With this, the sequence of numbers 011..., 011..., 011... etc. will be applied to regis- 
ter R. The resulting error will obviously equal 


cas dt. (24) 


This same error can be made equal to zero if the sequence of dy/dt signals has, as before, the form 0,1, 
0,1, ..., While the sequence of dx/dt signals is replaced by 1,0, 1,0, ... The phase error is particularly marked 
in cases when the functional blocks are connected in a circuit shunted by a feedback path. 


Rapidity of Digital Models 


An estimate of the rapidity of digital models can '... made on the basis of the following considerations. 
Let it be necessary, as the result of solving a given problem, to obtain some function of time x(t) (periodic or 
admitting an artificial periodic extension) which, with sufficient accuracy, can be presented in the frequency in- 
terval 0 = w = Wax by the Fourier series 





x(t) =2+ >) A, sin (wt + %) = >) %,, (25) 
i=} img) 


each of whose harmonics will be presented in the machine by its derivative with respect to time 


dr 
a = A; @, 008 (w, t + 9) (26) 
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The quantity x; can be obtained in the digital model with some relative error ¢ only in case the selected 
system of scales satisfies the relationship 


A, 1 
hel cr 


In accordance with (26), the product Ayw,/5x; determines the greatest velocity of increment application 
necessary for forming the current value x;. In digital models it must obey the obvious condition 


A; < 1 (28) 
— wo —_ 
85, ™ Ly) 


From (27) and (28) we obtain 


baa | | 
“= . 4 ; (29) 


In the literature, the quantity 1/6, is usually measured by the number of “iterations” per second where, by 
iteration, one understands the time interval 5, during which one elementary operation is executed in all blocks 
of the machine. For given e, the quantity 


© max ec 
Imax = Qn S in’, (3°) 


characterizes the passband of a low-frequency filter equivalent to the given machine. In serial digital differen- 
tial analyzers, the quantity 1/6, is of the order of 50 to 100 and, in parallel machines, goes as high today as 10°, 


With an accuracy of 0.01% for fray one obtains, respectively, 1.5- 107° cps for series machines, and 1.5 cps for 
parallel machines. 


Today, electronic computers are conventionally divided into two basic classes ~ digital and analog. If we 
accept the definition suggested by Yu. Ya. Basilevskii* for a mathematical machine as a “device which operates 
in accordance with a system of noncontradictory logical rules, accepts initial information, transforms it in ac- 
cordance with a given algorithm and then puts out the resulting information.” then a natural criterion for placing 
a given machine in one of the two classes will be the method used in it for coding information. The numerical 
codes used in digital machines are uniform codes with relatively low radices (almost always 2, 8, or 10) and high 


positionality. In analog machines (contimuous or discrete), codes are used with very high radices and low posi- 
tionality (usually 1). 


The deltamodulation method falls in the second group of codes and, on this basis, digital models must be 
assigned to the class of analog computers. The relatively high (as compared with the majority of other machines 
of this class) accuracy of digital models is explained by the high noise stability characteristic of the deltamodula - 
tion method, and the simplicity of exact reciprocal transformations of the deltamodulation code and numerical 
codes, which guarantee accuracy in introducing initial data and in outputing results. 


The correctness of assigning digital models to the class of analog devices (and not isolating them, for ex- 
ample, in an independent class) is borne out by a number of additional considerations. The method of preparing 
and setting up for a problem solution on a digital model is almost completely analogous with the corresponding 
methods used in operations with mechanical differential analyzers, and very close to the methods applied in opera - 
tions with electronic models (analog computers). One of the most characteristic features of general-purpose digi - 
tal computers is the possibility of executing individual operations, or small cycles of operations, whose results can 
depend as little as desired on the previous and following results. This capability is used for solving such logical 
problems as, for example, the translation of text from one language to another. In contradistinction to this, the 
solution of a problem on a digital model is obtained only as the sum of the results of a large number of elemen- 
tary operations. Finally, the presentation of the program in the form of codes sorted in the device's operational 
memory, and transformed together with the rest of the information, gives the general-purpose digital machine 
the ability to modify its programs which is completely lacking in digital models. 


* Questions In the Theory of Mathematical Machines. Vol. 1 (edited by Yu Ya. Basilevskii) [in Russian) (Fizmat- 
giz, 1958). 
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The logical operations cited in the literature [1,2] which are executed by digital models do not go beyond 
conditional transfers from one set of apparatus to another, or variations of initial conditions, which are also com- 
pletely realizable in ordinary analog machines. 


If one takes into account the heterogeneity of the functional blocks of digital models, and the compara - 
tively broad class of problems solvable by them [1], one will apparently be correct in classing them as general- 
purpose analog machines. A strict distinction should also be drawn between digital models and special - purpose 
digital computers, Thus, for example, the “digital differential analyzer” which realizes the iterative process of 
solving differential equations, as described in [1], is a special-purpose digital machine, since all the information 
in it is coded in a numerical code. 

SUMMARY 


1, Digital models are special forms of analog computers in which the quantities which present the original 
quantities of the problem are deltamodulation voltage pulses. 


2. In comparison with electronic analog computers, digital models have a much higher potential accuracy 
with a lower speed of action, The quantity w,,,,,/e has the same order of magnitude for digital models as for 
analog computers. 


3. AS concerns amount of apparatus and power requirements, digital models are roughly analogous to elec- 
tronic digital computers of the corresponding classes, although in parallel digital models, with their building - 
block design, the amount of apparatus depends on the complexity of the problem to be solved. Although differ- 
ing advantageously from general-purpose computers in the greater simplicity of programing problems for them, 
digital models have significantly less flexibility than general-purpose computers. 


4. Slowly flowing processes can advantageously be simulated on digital models. 
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PROCEEDINGS OF THE MOSCOW SEMINAR ON THE MAGNETIC 
ELEMENTS OF AUTOMATION DURING THE YEAR 1959 


G.V. Subbotina 


During the first half of 1959, the Moscow Seminar on magnetic amplifiers and contactless magnetic ele- 
ments of automation, remote control and computing technology for the Institute of Automation and Remote Con- 
trol of the AN SSSR continued its work. 


The session of January 28th consisted of the paper of M.A. Boyarchenkov, “Reversible dc instruments with 
magnetic amplifiers.” In considering the possible variants of control schemes for dc motors, the speaker based 
his work on the developed circuit with combined control of the magnitude of the armature current, by means of 
saturated reactors, and of the magnitude and direction of the exciter current. For this purpose, a two-cycle fast- 
acting magnetic amplifier was provided in the exciter circuit. 


On the 11th of February, N.P. Vasil'eva gave a paper on "Magnetic logical elements." The paper presented 
the state of the work being conducted abroad on magnetic logical elements, gave a classification of the known 
systems of logical elements and presented the results of the investigations and developments at the IAT AN SSSR 
in systems of logical elements of the types “and,” “or.” “not,” and “memory.” 


On February 25th, R.A. Lipman and I.B. Negnevitskii gave the communication, “Questions in the design of 
single-cycle magnetic amplifiers with internal feedback and dc loads.” The speakers suggested carrying out the 
computations for the gain of this type of amplifier by using experimental demagnetization curves. 


On March 11th, M.S. Libkind spoke on the development of a new type of three-phase ferromagnetic device 
for the control of ac power circuits without distortion of the form of the current curves. 


On March 25th, A.L. Pisarev, in the paper “Investigation of three-phase magnetic amplifier circuits,” spoke 
of the results of the analysis of the possible types of six-diode, three-phase magnetic amplifier circuits, and con- 
sidered their special features. The paper contained recommendations as to the use ofspecific circuits. 


The session of April 8th consisted of the presentation of the paper of K.D. Mart'yanova, “Temperature 
characteristics of industrial magnetically soft alloys." 


Two papers were given on April 22nd. In the paper, "The use of magnetic amplifiers for transforming 
single-phase current to three-phase current and vice versa,” A.M. Bamdas and V.A. Kulinich spoke of the develop- 
ment and analysis of the operation of original static electromagnetic transformers. In the paper of A.M. Bamdas 
and S,V. Shapiro, "New power actuating stabilizers with transformers controlled by shunt magnetization,” results 
in the development of actuator stabilizers were discussed. 


On May 13th, V.P. Molomin spoke of the work "Controlling three-phase asynchronous motors by means of 
magnetic amplifiers from a single-phase network.” 


On May 27th, V.S. Volodin and M.A. Rozenblat spoke on the development of sources of stable de voltages 
based on magnetic and semiconducting elements. 
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LIST OF LITERATURE ON QUESTIONS OF MATHEMATICAL 
SIMULATION (BY ANALOG COMPUTERS) FOR 1957 (continued) 


Compiled by E.O. Vil'td and R.S. Landsberg 
Edited by Cand. Tech. Sci. B. Ya. Kogan 


VI. Electromechanical Simulation Devices 





"Balanced-beam computing device.” Electronic and Radio Engr. 34, 9, 351 (Sept., 1957) 2 illus. 
"Electro-magnetic computer." Proc. Control and Automat. 4, 8, 311-312 (August, 1957) 4 illus. 
M. Fogiel, “Absolute value via a differential.” Instrum. and Automat. 30, 4, 707 (April, 1957). 


R. Gendreu, "Servomechanisms in analog computers.” [in French] Ann. Radioelectricité 12, 47, 45-61 
(Jan., 1957) 8 illus, 


P.J. Gravell and Th.F. Kavanagh, “Computers ~ the key to modern manufacturing scheduling.” Trans. IRE. 
PGIE — 4, 90-93 (March, 1957) 1 illus. 


S.L. Hess, "A simple analog computer for determination of the Laplacian of a mapped quantity.” Bull. 
Amer. Meteor. Soc. 38, 2, 67-73 (February, 1957). 


N. Hogben, "A new harmonic integrator." Trans. Instn Naval Archit. 99, 2, 171-172 (1957). 


W.E. Howland, E.A, Trabant and G.A. Hawkins, "A mechanical computing device for heat-conduction 
problems.” Trans. ASME 79, 3, 675-680 (April, 1957) 5 illus., 5 refs. 


E. Jucker, “A Swiss analog computer.” (in German] Bull. Assoc. Suisse Electriciens 48, 23, 1017-1020 (1957) 
15 illus. 


N. Rouche, "Electromechanical analog computers and the Lagrange-Maxwell equations.” [in French) Bull. 
Soc. Franc. Electric. 75, 7, 181-186 (March, 1957). 


J. Valat, "An electromagnetic simulator for a Hill's equation.” Compt. rend. Acad. Sci. Paris 244, 20, 
2462-2465 (May 13, 1957). 


VII. Special-Purpose Analog Computers 





1, Devices for Solving Systems of Algebraic Equations. Root Indicators. 





A.N. Lebedev, “Stability of inertia -less mathematical models.” Priborostroenie 7, 6-8 (1957) 4 illus., 4 
refs, 


H. Adler, “An electrical instrument for solving polynomial equations.” [in German} Nachrichtentechnik 7, 
8, 335-342 (August, 1957) 7 illus., 18 refs. 


V. Horvath, “Practical methods for extracting roots by successive integrations.” [in Czechoslovakian} 
Strojnoelektrotechn. Casor 8, 3, 189-197 (1957). 


S.K. Ip, "An electrolytic tank as an analog-computing machine for factorizing high-degree polynomia!s.” 
Quarterly J. Mech. and Appl. Math. 10, 3, 368-384 (August, 1957) 7 illus., 3 tables, 12 refs. 
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L. Levine and H.F. Mefesinger, “An automatic analog computer method for solving polynomials and find- 
ing root loci.” Convent. Rec. IRE 5, 4, 164-172 (March 18-21, 1957) 7 illus., 8 refs. 


2. Devices for Solving Integral Equations 





M.E, Fisher, “On the continuous solution of integral equations by an electronic analog. Pt. 1." Proc. Cambr, 
Philos. Soc. 53, 162-174 (Jan., 1957). 


3. Correlators 





Osima Sintaro, Takasima Kenenié, Nisilati Kokki and Vatanabé Khirosi, "An analog-type correlator.” (in 
Japanese) Electr. Commun. Lab. Techn. J. 6, 3-4, 203-230 (1957); Elektrotekhnika 7, 159 (1959). 


C. van Schooneveld, “An electronic analog correlator." [in Dutch] Tijdschr. Nederl. Radiogen 22, 4, 205- 
237 (1957). 


4. Various Devices 





S.P. Kapitsa, “The mechanical computation of harmonic conjugate functions,” in the book: Computing 
Technology, Vol. 1 [in Russian] (AN SSSR, Moscow, 1957) pp. 167-169, 1 illus., 2 refs. 


V.A. Frank, “Fourier analog-computer of the Hagg-Laurent type.” J. Scient. Instrum. 34, 5, 210-211 (May, 
1957) 2 illus.,2 refs. 


R.F. Scott, “Hydraulic analog computer for studying diffusion problems in soil” Geotechnique 7, 2, 55-72 
(June, 1957). 


VIII. Analog-Digital and Digital-Analog Converters 





V.M. Kuntsevich, “Analog-to-digital converters (a survey).” Avtomatika 2, 70-87 (1957) 10 illus., 24 refs. 


C.R, Borley, C,H. Braybrook and L. Coates, "Automatic positioning systems for machine tools.” J. Brit. 
Instn Radio Engrs 17, 9, 513-521 (Sept., 1957) 8 illus., 17 refs. 


G.G, Bower, “An analog~-digital converter.” [in German] Regelungstechnik 5, 11, 418-421 (1957) 4 illus.; 
5, 12, 459-466 (1957) 19 illus., 22 refs.; Control Engng 4, 4, 107-118 (1957). 


Discussion on “High-speed analog -to-digital converters" by G.I. Herring and D. Lamb. J. Brit. Instn Radio 
Engrs 17, 11, 631 (Nov., 1957). 


S.K. Feingold, “The logic of V-brush analog -to-digital converters.” J. Amer. Instrum. Soc. 4, 2, 66-68 
(February, 1957) 4 illus., 1 ref. 


H.J. Finden and B.A. Hotlock, “The inductosyn and its application.” J. Brit. Instn Radio Engrs 17, 7, 369- 
383 (July, 1957) 21 illus. 


H.W, Gettings, “Analog-digital analog recording." Automat. Control 4, 4, 20-21 (1957). 


A.D. Glick, “High-speed digital-to-analog conversion by integration of a variable-rate pulse train." Proc. 
Western Joint Computer Conference, Los Angeles, Calif., February 26-28, 1957. (New York, IRE, 1957) pp. 128- 
133, 9 illus, 


G.J. Herring and D. Lamb, “High-speed analog -to-digital converters.” J. Brit. Instn Radio Engrs 17, 8, 407- 
420 (August, 1957) illus., 16 refs. 


M.L. Klein, F.K, Williams and H,C, Morgan, “Digital transducers.” Instrum. and Automat. 30, 7, 1299- 
1300 (July, 1957) 2 illus. 


G.H, Myers, “A cyclic digital-to-analog decoder.” Convent. Rec. IRE 5, 4, 156-163 (March 18-21, 1957) 
9 illus., 2 refs. 


H.N. Putschi, J.A. Raper and J.J. Suran, “Digital-analog converter provides storage." Electronics I, 30, 12, 
148-151 (Dec., 1957) 4 illus. 









































H. Stern, “General-purpose integrator for date reduction.” Instrum. and Automat. 30, 2, 254-255 (February, 
1957) 3 illus, 


A. Tiffany, “A simple shaft digitizer and store.” Electronic Engng 29, 358, 568-574 (Dec., 1957) 10 illus., 
4 tables, 18 refs. 


IX. Applications of Analog Computers 





1, Applications of Analog Computers for the Solution of Automatic Control Problems 





A.M. Batkov, The Solution of One Class of Nonstationary Problems of the Statistical Dynamics of Auto- 
matic Control Systems Using Electronic Analog Computers. Dissertation [in Russian] (Mosk. Inzh, -Fiz. In-t., Mos- 
cow, 1957). 


Yu.E. Efroimovich, “Analog and control computers for steel arc furnaces,” Stal’ 7, 602-608 (1957) 7 illus, 


G.A, Tatsitov, “Obtaining Mikhailov curves graphically,” in the book: Collected Papers on Electric Instru- 
ment Construction [in Russian] (Leningrad, 1957) pp. 94-110, 11 illus., 10 tables, 10 refs. (Leningrad In-t. Techn. 
Mekh. i Optiki, 28.) 


G.A, Bekey, “Analog computers in process control.” Canad. Chem. Process 41, 3, 90-92, 94 (1957). 


H.G. Doll and T.M. Stout, "Design and analog-computer analysis of an optimum third-order nonlinear 
servomechanism.” Trans. ASME 79, 3, 513-525 (April, 1957) 24 illus., 25 refs. 


G, Gandusio and F. Guerrini, “A simple analog computer for industrial control.” [in Italian] Strumentazione 
e Automazione 6, 108-111 (March, 1959) 6 illus. 


E.L. Harder, “Analog computers in industrial control systems,” in: Automation in Business and Industry, 
edited by E.M. Grabbe (J. Wiley and Sons, Inc., London, 1957) pp. 456-493, 23 illus., 6 refs. 


T.J. Higgins, “Electroanalogic methods,” Appl. Mech. Revs 10, 2, 49-54 (February, 1957); 8, 331-335 
(August, 1957) 712 refs. 


H.H, Idzerda and L. Ensing, “An electronic simulator for solving control problems." [in German] Regelungs- 
technik 5, 6, 199-205 (1957) 22 illus. 


J.M.L. Janssen and R.P. Offereins, "Use of an electronic simulator for finding the optimal tuning of regulators.” 
[in German] Regelungstechnik 5, 8, 264-270 (1957) 16 illus., 9 refs. 


J.M.L, Janssen, "Use of an analog computer in solving the control problems of the Shell Laboratory.” [in 
Dutch] Ingenieur 69, 32, 083-087 (August 9, 1957) 14 illus., 6 refs. 


G. Kourim, “Superheaters and regulators in models of control systems.” [in German] Regelungstechnik 5, 
12, 466-469 (1957) 8 illus., 6 refs. 


K.J. Lesemann, “Applications of analog computers in the design of control systems.” (in German] Regelungs~- 
technik. Report on the Proceedings in Heidelberg, Sept., 1956 [in German] (Oldenbourg Verlag, Munich, 1957) 
pp. 333-342, 16 illus., 8 refs. 


E.H. Ludwig, "The use of analog devices as control devices enlarges the capabilities of control technology.” 
[in German] Regelungstechnik. Report on the Proceedings in Heidelberg, Sept., 1956 [in German] (Oldenbourg 
Verlag, Munich, 1957) pp. 352-357, 8 illus,, 1 ref. 


S.W. McEihenny, “Design and analysis of the motor-alternator speed regulator for the New Haven speed- 
merchant locomotive.” Applicat. and Ind. 29, 10-16 (March, 1957) 16 illus. 


B. Mazelsky and H.B. Amey, “On the simulation of random excitations for airplane response investigations 
on analog computers.” J. Aeronaut. Sci. 24, 9, 633-649 (Sept., 1957). 


].G. Miles and W.G, Proctor, *The use of an analog computer in control system analysis.” Brit. Commun, 
and Electronics 4, 2, 86-90 (1957) 13 illus. 


D.R. Miller, *Problem solving.” Petrol. Engr 29, 3, E12, 14-6 (March, 1957). 
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A, Nathan and J. Mohler, *Demonstration of conditional stability on an analog computer.” Trans. IRE EC-6, 
4, 287 (Dec., 1957) 4 illus, 


D.W. Russell, W.K. McGregor and L.F. Burns, “The analog computer as a process controller.” Control Engng 
4, 9, 160-165 (Sept., 1957). - 


J.W. Schwartzenberg, “An analog study of a high-speed recording servomechanism.” Trans. ASME 80, 2, 
490-496 (February, 1958) 14 illus., 6 refs. 


J.\.. Shearer, “Nonlinear analog study of a high-pressure pneumatic servomechanism.” Trans. ASME 79, 3, 
465-472 (April, 1957) 12 illus., 9 refs. 


"Synchronous -machines analogs.” J. Inst. Electr. Engrs 3, 27, 164-165 (March, 1957). 


W.T. Thomson, “Analog computer for nonlinear system with hysteresis.” J. Appl. Mech. 24, 2, 245-247 
(June, 1957) 5 illus., 3 refs. 


Ch, Vazaca and G, Janculescu, “An electronic analog automatic control system.” [in Rumanian] Automat. 
i electron. 1, 5, 200-208 (1957). 


M.F, Versini, "The use of analog computers in solving problems of industrial control systems.” [in German] 
Dtsch. Elektrotechn, 11, 11, 499-500 (1957) 10 illus. 


Ch.W. Worley, “Simulation of a liquid level control system with an analog computer." Computers and 
Automat, 6, 2, 30-31 (1957). 


2. Use of Analog Computers and Elements in Aviation 





A.S. Lock, “Simulation, computers and telemetry," Missile Control [Russian translation} (Gostekhteorizdat, 
Moscow, 1957) pp. 692-759, 39 illus., 2 tables. 


“Admiralty tactical simulator.” Electr. J. 158, 11, 766 (1957). 
* Analog computers simulate aircraft take-off." Brit. Communs, and Electronics 4, 1, 40 (1957) 1 illus. 


W.R. Blunden, F.H. Hooke and H.K. Messerle, “A study of aircraft response to sudden and graded wind gusts 
by means of a differential analyzer.” Brit. J. Appl. Phys. 8, 4, 167-173 (April, 1957) 9 illus., 8 refs. 


H, Constant, “Simulated flight testing of aeroengines.” Engineering 183, 4762, 750-753 (June 14, 1957) 
6 illus. 


D.L. DeMyer, “Flight simulator tests fire-control radars.” Electronics 30, 6, 168-170 (June, 1957) 2 illus. 
R.K. Frick, “Analog computers in propeller control design.” Electr. Engng 76, 10, 858 (Oct., 1957). 


R.K. Frick, “The role of analog computers in propeller control design.” Applic. and Ind. 33, 257-263 (Nov., 
1957) 9 illus., 6 refs, 


J. Hudson and K.K. Smell, “Simulators in Navy missile training.” Missiles and Rockets 2, 1, 57-58 (Jan., 1957). 


G. Isakson, “Simple model study of transient temperature and thermal! stress distribution due to aerodynamic 
heating.” J. Aeronaut. Sci. 24, 8, 616-619 (August, 1957). 


A.D. Jeffery, LM. Dempsey, H. Saville and K. Gill, "The application of analog computer technique to the 
design of aero engine control systems.” J, Brit. Inst. Radio Engrs 17, 11, 633-647 (Nov., 1957) 19 illus. 


N.F. Kfoury, “Analog computer control compound aircraft machining.” Electr. Engng 76, 10, 915-919 
(Oct., 1957) 14 illus. 


E. Lyon and D.S. Peck, “AC or DC computers for flight simulators?” Aviat. Age 28, 6, 70-75 (1957). 


W. Makinson and G.M, Hellings, “Synthetic aids to flying.” J. Roy. Aeronaut. Soc. *1, 560, 509-528 
(August, 1957). 


“A radar trainer incorporating an analog computer.” Brit. Commun. and Electronics 4, 2, 101-102 (Feb- 
ruary, 1957) refs. 














J.V. Roberts, “Short's helicopter flight simulator,” Brit. Communs and Electronics 4, 9, 542-545 (Sept., 1957) 


5 illus, 


"Simulating helicopter flight,” Engineering 184, 4780, 500-501 (Oct. 18, 1957) illus. 
"Simuiation-key to systems development,” SAE Journal 65, 12, 30-31, 36, 37 (Nov., 1957) 5 illus., 1 table. 
"Simulation kinetic heating in high-speed aircraft,” Engineering 184, 4768, 111 (July 26, 1957). 


F.R. Spearman, "Analog computing applied to guided weapons,” Discovery 18, 5, 192-197 (May, 1957) 7 
illus, 


"A three~-axis simulator for controlled flight test of airborne systems,” Computers and Automat. 6, 5, 18 
(May, 1957). 


B.H. Venning, “Analog computer development with reference to helicopter applications,” J. Helicopter 
Assn 11, 2, 77-103 (1957). 


V.IL. Wethe, “Airborne computers will reduce cockpit workload,” Aero-Engng Rev. 16, 75-78 (May, 1957). 
R.W. Williams, “The applications of analog computers in the development of guided weapons,” Engl. Elec- 
tric J. 15, 1, 34-45 (March, 1957) 12 illus. 


3. Various Applications 





V.V. Aleksandrov, “The use of electronic mathematical machines for the automation of chemical produc- 
tion,” Khimicheskaya Prom-st' 2, 52-56 (1957) 7 illus., 5 refs. 


S.A. Aleskerov, "On designing electromagnetic systems on electrical models," Automation and Remote 
Control (USSR) 18, 764-772 (1957) 4 illus., 8 refs. 


LS. Bruk, “On problems in the domain of automating the control of energy system modes of operation, and 
the prospective applications of computing technology,” Elektrichestvo 12, 31-34 (1957). 


A.L. Vainer and S.M. Fertik, "The use of electrical simulation in the design of lightening-protection ground- 
ing devices,” Trudy Khar'k Politekhn. In-ta, Ser. Elektroradiotekhn. 12, 2, 179-203 (1957) 26 refs. 


Yu.E£. Efroimovich, “Analog and control computers for steel arc furnaces,” Stal’ 7, 602-608 (1957) 7 illus. 


V.P. Kononov, “Electrical simulation of steady-state modes in complex heat-fixing systems by means of 
electronic simulators,” in the book: Trudy 4-i Konfer. Molodykh Uchenykh (Moscow, 1957) pp. 172-190, illus., 
6 refs. (Energet. In-t. AN SSSR.) 


D.P. Morozov and Yu.A. Bortsov, “Simulating transient responses in complex electric drive circuits,” in the 
book: Trudy Mosk. Energet. In-ta. Vyp. 29 (Gosénergoizdat, Moscow-Leningrad, 1957) pp. 126-143, illus. 


"Radar trainer with analog computers,” Vopr. Radiolokats. Tekhniki 3, 88-91 (1957) 2 illus. 
“Type SL21 radar trainer,” Vopr. Radiolokats. Tekhniki 1, 82-88 (1957) 5 illus., 2 refs. 


A.S. Aldred and P.A. Doyle, “Electronic-analog-computer study of synchronous-machine transient stability,” 
Proc, Instn Electr. Engrs 104, A, 14, 152-160 (April, 1957) 15 illus.,18 refs. 


"Analog computer for admiralty research at Imperial College,” Brit. Communs and Electronics 4, 5, 302 
(May, 1957); Engineer 203, 5276, 376 (1957). 


C.P. Atkinson, "Superharmonic oscillations as solutions to Duffing*s equation as solved by an electronic 
differential analyzer,” J. Appl. Mech, 24, 4, 520-528 (Dec., 1957) 8 illus., 9 refs. 


L.V. Boffi and V.B. Haas, “Analog computer representation of alternators for parallel operations,” Commun. 
and Electronics 30, 153-157 (May, 1957) 6 illus., 2 refs, 


P.E.A. Cowley, “The application of an analog computer to the measurement of process dynamics,” Trans. 
ASME 179, 4, 823-832 (1957) 25 illus., 13 refs. 


R. DeCote and W.J. Horvath, “An electronic computer for vector electrocardiography,” Trans, IRE PGME-8, 
31-37 (July, 1957) 5 illus., 1 ref. 






























N.D. Diamantides, “Analog computer solves statistical problems,” Automat. Control 4, 5, 60-65 (1957). 


N.D. Diamantides and M. Horowitz, “Autocorrelation of the earth's crust with analog computers,” Rev. 
Scient. Instrum. 28, 5, 353-360 (May, 1957) 7 illus. 


F.D. Ezekiel and H.M. Paynter, “Computer representation of engineering systems involving fluid transients,” 
Trans. ASME 79, 8, 1840-1850 (Nov., 1957) 15 illus., 19 refs. 


V. Frank, “A Fourier analog computer of the Hagg-Laurent type,” J. Scient. Instrum. 34, 5, 210-211 (1957) 
2 Illus,, 2 refs. 


W. Fritzsche, “Computers as automatic drive elements” [in German] AEG Mitt. 47, 11/12, 428-432 (1957) 
5 illus., 8 refs. 


"The G.E.C, electronic computer service,” Ind. Progr. and Developm. 18, 2, 51 (1957). 
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